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Since the subject of traffic dynamics has captured the interest of physicists, many surprising effects 
have been revealed and explained. Some of the questions now understood are the following: Why 
are vehicles sometimes stopped by so-called "phantom traffic jams", although they all like to drive 
fast? What are the mechanisms behind stop-and-go traffic? Why are there several different kinds of 
congestion, and how are they related? Why do most traffic jams occur considerably before the road 
capacity is reached? Can a temporary reduction of the traffic volume cause a lasting traffic jam? 
Under which conditions can speed limits speed up traffic? Why do pedestrians moving in opposite 
directions normally organize in lanes, while similar systems are "freezing by heating"? Why do 
panicking pedestrians produce dangerous deadlocks? All these questions have been answered by 
applying and extending methods from statistical physics and non-linear dynamics to self-driven 
many-particle systems. This review article on traffic introduces (i) empirically data, facts, and 
observations, (ii) the main approaches to pedestrian and vehicle traffic, (iii) microscopic (particle- 
based), mesoscopic (gas-kinetic), and macroscopic (fluid-dynamic) models. Attention is also paid 
to the formulation of a micro-macro link, to aspects of universality, and to other unifying concepts 
like a general modelling framework for self-driven many-particle systems, including spin systems. 
Subjects like the optimization of traffic flows and relations to biological or socio-economic systems 
such as bacterial colonies, flocks of birds, panics, and stock market dynamics are touched as well. 
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I. INTRODUCTION 

A. Motivation: History of traffic modelling and traffic 
impact on society 

The interest in the subject of traffic dynamics is 
surprisingly old. In 1935, Greenshields has carried out 
early studies of vehicular traffic, and already in the 
1950ies, a considerable publication activity is docu- 
mented by journals on operations research and engi- 
neering. These papers have already addressed topics 
like the fundamental diagram between traffic flow and 
vehicle density or the instability of traffic flow, which 
are still relevant. The reason becomes understandable 
through the following quotation from H. Greenberg in 
1959: 

"The volume of vehicular traffic in the past several 
years has rapidly outstripped the capacities of the na- 
tion's highways. It has become increasingly necessary 
to understand the dynamics of traffic flow and obtain 
a mathematical description of the process. " 

More than fourty years later, the situation has de- 
teriorated a lot. Some cities like Los Angeles and San 
Francisco suffer from heavy traffic congestion around 



the clock. In Europe, as well, the time that drivers 
spend standing in traffic jams amounts to several days 
each year. During holiday seasons, jams may grow up 
to more than 100 kilometers in size. Vehicle emissions 
like SO2, NOa;, CO, CO2, dust particles, smog, and 
noise have reached or even exceeded a level compa- 
rable to those by industrial production and those by 
private households, being harmful to the environment 
and human health. On average, every driver is in- 
volved in one accident during his lifetime. In Germany 
alone, the financial damage by traffic due to accidents 
and environmental impacts is estimated 100 billion 
dollars each year. The economic loss due to congested 
traffic is of the same order of magnitude. However, 
the traffic volume is still growing because of increased 
demands for mobility and modern logistics, including 
e-commerce. 

Without any doubt, an efficient transportation sys- 
tem is essential for the functioning and success of mod- 
ern, industrialized societies. But the days, when free- 
ways were free ways, are over. Facing the increasing 
problems of individual traffic, the following questions 
came up: Is it still affordable and publicly acceptable 
to expand the infrastructure? Will drivers still buy 
cars in view of streets that effectively turn into park- 
ing lots? It is for these reasons that automobile com- 
panies started to worry about their own future and 
decided to spend considerable amounts of money for 
research on traffic dynamics and the question how the 
available infrastructure could be used more efficiently 
by new technologies (telematics). 

At this time, physicists were motivated to think 
about what physics could contribute to the field 
of traffic dynamics. Although, among mathemati- 
cians, physicists, and chemists, there were already 
some early pioneers like Whitham, Prigogine, Mon- 
troll, and Kiihne, the main activities started in 
1992 and 1993 with papers by Biham et al. (1992), 
Nagel and Schreckenberg (1992), as well as Kerner 
and Konhauser (1993). These papers initiated an 
avalanche of publications in various international 
physics journals. Meanwhile, it is difficult to keep 
track of the scientific development and literature. 
Therefore, it is high time to write a review on traffic, 
which tries to bring together the different approaches 
and to show up their interrelations. Doing so, I will 
take into account many significant contributions by 
traffic engineers as well. Hopefully, this will stimulate 
more intense discussions and cooperations among the 
different disciplines involved. 

In the following sections, I will try to give an 
overview of the field of traffic dynamics from the per- 
spective of many-particle physics. 

B. Driven many-particle systems 

A big challenge for physicists arc the dynamics and 
pattern formation in systems far from equilibrium, es- 



pecially living systems (Nicolis and Prigogine, 1977; 
Hakcn, 1977, 1983, 1988; Pasteels and Deneubourg, 
1987; Weidlich, 1991; Feistel and Ebeling, 1989; Vic- 
sek, 1992; Kai, 1992; DeAngelis and Gross, 1992; Val- 
lacher and Nowak, 1994; Cladis and Palffy-Muhoray, 
1995; Helbing, 1995a). Non-equilibrium systems are 
characterized by not being closed, i.e., by having an 
exchange of energy, particles, and/or information with 
their environment. Consequently, they often show 
complex behavior and, normally, there are no general 
results such as the laws of thermodynamics and statis- 
tical mechanics for closed systems of gases, fluids, or 
solids in equilibrium. Nevertheless, they can be cast 
in a general mathematical framework (Helbing, 2001), 
as I will show here. 

1. Classical mechanics, fluids, and granular media 

Let us start with Newton's equation of motion from 
classical mechanics, which describes the acceleration 
x a (t) of a body a of mass m a subject to pair interac- 
tions with other bodies 0: 

m a x a (t) = ^2 F ap(t) • (1) 

/3(#«) 

The interaction forces F a p(t) are mostly dependent 
on the locations x a (t) and x@{t) of the interacting 
bodies a and at time t. Often, they depend only 
on the distance vector d a p = (xp — x a ), but in spe- 
cial cases, they are also functions of the velocities 
v a(t) = x a (t) and vp(t) — x@(t). For potential forces, 
the above many-body system can be characterized by 
a Hamilton function. A typical example is the descrip- 
tion of the motion of celestial bodies. 

In driven many-body systems such as fluids un- 
der the influence of pressure gradients and boundary 
forces or vibrated granular media like sand, we have 
to consider additional interactions with the environ- 
ment. Therefore, we need to consider additional terms. 
This concerns (external) driving forces F (x,t) due 
to boundary interactions and gravitational or electri- 
cal fields, (sliding) friction forces Ff T (t) — —j a v u(t) 
with friction coefficient j a , and individual fluctuations 
Ca (t) reflecting thermal interactions with the environ- 
ment (boundaries, air, etc.) or a variation of the sur- 
face structure of the particles: 

m a x a (t) = F (x a (t), t) - j a v a (t) + ^2 F <xf){t) 

+ Ca(*). (2) 

From these so-called "microscopic", molecular- 
dynamics type of equations, one can systematically 
derive "macroscopic" , fluid-dynamic equations for the 
spatio-temporal evolution of the particle density, the 
momentum density or average velocity, and the en- 
ergy density or velocity variance (related to the tem- 
perature). The methods for the construction of this 
micro-macro link will be sketched in Sec. III.E. 



In driven systems, the ongoing competition between 
the driving forces and the dissipative friction forces 
leads to a spatio-temporal redistribution of energy, 
which produces a great variety of self-organization 
phenomena. This results from non-linearities in the 
equations of motion, which allow small initial pertur- 
bations to be enhanced and non-equilibrium patterns 
to be dynamically stabilized. In fluids one can find 
the formation of waves or vortices, bifurcation scenar- 
ios like period doubling behavior, the Ruelle-Takens- 
Newhouse route to chaos, intermittency, or turbu- 
lence, depending on the particular boundary condi- 
tions (Joseph, 1976; Landau and Lifshits, 1987; Drazin 
and Reid, 1981; Swinney and Gollub, 1985; Schuster, 
1988; Grofimann, 2000). In this review, it is important 
to know that turbulence normally requires three- or 
higher-dimensional systems, so it is not expected to 
appear in one- or two-dimensional vehicle or pedes- 
trian traffic. However, the instability mechanism for 
the explanation of the subcritical transition to turbu- 
lence in the Hagen-Poiseuille experiment (Gebhardt 
and Grofimann, 1994; Grofimann, 2000) may also be 
relevant to traffic systems, as pointed out by Krug (cf. 
the discussion of metastability in Sec. IV. A. 4). 

In vibrated granular media, one can find emer- 
gent convection patterns (Ehrichs et al, 1995; 
Bourzutschky and Miller, 1995; Poschel and Herr- 
mann, 1995), collective oscillating states (so-called os- 
cillons, see Umbanhowar et al, 1996), spontaneous 
segregation of different granular materials (Poschel 
and Herrmann, 1995; Santra et al, 1996; Makse et al, 
1997), or self-organized criticality with power-law dis- 
tributed avalance sizes in growing sand heaps (Bak 
et al, 1987, 1988; Bak, 1996) or in the outflow from 
hoppers (Schick and Verveen, 1974; Peng and Herr- 
mann, 1995). 

In spite of many differences between flows of flu- 
ids, granular media, and vehicles or pedestrians due 
to different conservation laws and driving terms, one 
can apply similar methodological approaches, e.g. 

(i) microscopic, molecular dynamic models (see, e.g., 
Hoover, 1986; Goldhirsch et al, 1993; Buchholtz and 
Poschel, 1993; Hirshfcld et al, 1997; Sec. III. A), 

(ii) lattice gas automata (Frisch et al, 1986; Chen 
et al, 1991; Peng and Herrmann, 1994; Tan et al, 
1995) or cellular automata (see Sec. III.B), 

(iii) gas-kinetic (Boltzmann- and Enskog-like) models 
(Boltzmann, 1964; Enskog, 1917; Chapman and Cowl- 
ing, 1939; Cohen, 1968, 1969; Lun et al, 1984; Jenk- 
ins and Richman, 1985; Kobryn et al, 1996; Dufty 
et al, 1996; Lutsko, 1997; Cercignani and Lampis, 
1988; McNamara and Young, 1993; Sela et al, 1996; 
see Sec. III.E), and 

(iv) fluid-dynamic models (Haff, 1983; Goldhirsch, 
1995; Sela and Goldhirsch, 1995; Hayakawa et al, 
1995; Du et al, 1995; see Sec. III.D). 



C. Self-driven many-particle systems 

In self-driven many-particle systems, the driving 
force is not of external origin (exerted from out- 
side), but associated with each single particle and self- 
produced. This requires each particle to have some 
kind of internal energy reservoir (Schweitzer et ai, 
1998b; Ebeling et ai, 1999). 

Self-driven "particles" are a paradigm for many ac- 
tive or living systems, where they are a simplified 
and abstract representation of the most important dy- 
namic behavior of cells, animals, or even humans. In 
order to reflect this, I will generalize Eq. (2) a little, 
replacing the external driving force F (x a ,t) by an 
individual driving force F^(t). Moreover, Newton's 
third law Fp a (t) — —F a p(t) (actio = reactio) does 
not necessarily apply anymore to the "self-driven" , 
"self-propelled", "motorized", or "active" "particles" 
we have in mind. I will show that these minor changes 
will imply various interesting phenomena observed in 
nature, for example in biological, traffic, or socio- 
economic systems. In this context, the masses m a are 
sometimes not well-defined, and it is better to rewrite 
the resulting equation by means of the scaled quan- 
titles F° a (t) = 7a«2(*)e£(*), 7a = m Q /r Q , F aP {t) = 
Tn a f a fj(t), and Ca(t) = J a £u(t), where the accelera- 
tions f a p{t) are often loosely called forces as well: 



dv a (t) _ v° a {t)e° a (t) + £ a (t) - v a {t) 



dt 



E 



(3) 



From this equation we can see that, with a relaxation 
time of T a , the driving term w° (t)e° a (t)/r a and fric- 
tion term —v a (t)/r a together lead to an exponcntial- 
in-time adaptation of the velocity v a (t) to the desired 
speed w° (t) and the desired direction e a (t) of motion. 
This is, however, disturbed by fluctuations £ Q (i) and 
interactions f a /3(t) with other particles (3. It is clear 
that attractive forces f a p(t) will lead to agglomer- 
ation effects. Therefore, we will usually investigate 
systems with vanishing or repulsive interactions, for 
which one can find various surprising effects. 

A further simplification of the model equations can 
be reached in the overdamped limit r Q of fast 
(adiabatic) relaxation. In this case and with the ab- 
breviations V afj (t) = T a f al3 (t), £ a (t) = v a a (t)Xa{t), 

we obtain 

v a (t) = v a (t)e a (t) = v° a (t)e° a (t) + v a p(t) + £ a (t) 
= v Q a {t)[e° a {t)+Xa(t)] + Y l v a p(t). (4) 



1. The concept of social (behavioral) forces 

Human behavior often seems to be "chaotic" , irreg- 
ular, and unpredictable. So, why and under which con- 



ditions can we apply the above force equations? First 
of all, we need to be confronted with a phenomenon of 
motion in some (quasi-)continuous space, which may 
be also an abstract behavioral space or opinion scale 
(Helbing, 1992a, 1993b, 1994, 1995a). It is favourable 
to have a system where the fluctuations due to un- 
known influences are not large compared to the sys- 
tematic, deterministic part of motion. This is usually 
the case in pedestrian and vehicle traffic, where people 
are confronted with standard situations and react au- 
tomatically rather than making complicated decisions 
between various possible alternatives. For example, an 
experienced driver would not have to think about the 
detailled actions to be taken when turning, accelerat- 
ing, or changing lanes. 

This automatic behavior can be interpreted as the 
result of a learning process based on trial and error 
(Helbing et ai, 2001c), which can be simulated with 
evolutionary algorithms (Klockgether and Schwefel, 
1970; Rechenberg, 1973; Schwefel, 1977; Baeck, 1996). 
For example, pedestrians have a preferred side of walk- 
ing (Oeding, 1963; Older, 1968; Weidmann, 1993), 
since an asymmetrical avoidance behavior turns out to 
be profitable (Bolay, 1998). The related formation of 
a behavioral convention can be described by means of 
evolutionary game theory (Helbing, 1990, 1991, 1992a, 
c, 1993a, 1995a, 1996c). 

Another requirement is the vectorial additivity of 
the separate force terms reflecting different environ- 
mental influences and psychological factors. This is 
probably an approximation, but there is some exper- 
imental evidence for it. Based on quantitative mea- 
surements for animals and test persons subject to sep- 
arately or simultaneously applied stimuli of different 
nature and strength, one could show that the behavior 
in conflict situations can be described by a superposi- 
tion of forces (Miller, 1944, 1959; Herkner, 1975). This 
fits well into a concept by Lewin (1951), according to 
which behavioral changes are guided by so-called so- 
cial fields or social forces, which has been put into 
mathematical terms by Helbing (1991, 1992a, 1993b, 
1994, 1995a; see also Helbing and Molnar, 1995). The 
third Newtonian law, however, is usually not valid. 



D. What this review is about 

In the following sections, 

(i) I will focus on the phenomena actually observed in 
traffic and their characteristic properties, while I will 
discuss models only to the extend, to which they are 
helpful and necessary for an interpretation and a bet- 
ter understanding of the observations, 

(ii) I will discuss the main methods from statistical 
physics relevant for modelling and analyzing traffic 
dynamics (see the table of contents), 

(iii) I will discuss how far one can get with a physical, 
many-particle description of traffic, neglecting socio- 
psychological factors and human behavior. The limi- 



tations of this approach will be shortly discussed in 
Sec. IV.F. 

This review is intended to serve both, experts and 
newcomers in the field, so some matters will be sim- 
plified or explained in more detail for didactical rea- 
sons. Moreover, I will try to identify open problems 
and to shed new light on some controversial topics 
presently under discussion. The main focus will be 
on the various kinds of phenomena occuring in self- 
driven many-particle systems and the conditions, un- 
der which they appear. I will start with the subject 
of one-dimensional vehicle traffic, continue with two- 
dimensional pedestrian traffic and three-dimensional 
air traffic by birds. 

Those physicists who don't feel comfortable with 
non-physical systems may instead imagine special 
kinds of granular, colloidal, or spin systems driven by 
gravitational or electrical forces, or imagine particu- 
lar systems with Brownian motors (Hanggi and Bar- 
tussek, 1996; Astumian, 1997; Julicher et al, Prost, 
1997; Reimann, 2000). Moreover, I would like to en- 
courage everyone to perform analogous experiments 
in these related physical systems. 

Despite the complexity of traffic and complications 
by unknown, latent, or hardly measurable human fac- 
tors, physical traffic theory is meanwhile a prime ex- 
ample of a highly advanced quantitative description 
of a living system. An agreement with empirical data 
has been reached not only on a qualitative, but even 
on a semi-quantitative level (see Sec. IV.B.f). More- 
over, there are even "natural constants" of traffic, 
emerging from the non-linear vehicle interactions, see 
Sec. IV. A. 5. But, before getting into all this, I would 
like to mention some other interesting (self-)driven 
many-particle systems loosely related to traffic. 

E. Particle hopping models, power-law scaling, and SOC 

Many-particle systems in equilibrium can be very 
well understood with methods from thermodynam- 
ics (Keizer, 1987) and statistical physics (Uhlen- 
beck and Ford, 1963; Landau and Lifshits, 1980; 
Ma, 1985; Klimontovich, 1986; Huang, 1987). Exam- 
ples are phase transitions between different aggregate 
states of matter like vapor, water, and ice, or the mag- 
netization of spin systems composed of many "elemen- 
tary magnets" . It is known that the phase transition 
of disordered many-particle systems with short-range 
interactions into states with long-range order is easier 
in higher dimensions, because of the greater number 
of neighbors to interact with. Often, there is a cer- 
tain upper dimension, above which the system can be 
described by a mean field approach so that the fluctu- 
ations can be neglected. For lower-dimensional spaces, 
one can usually develop approximate theories for the 
influence of noise by means of suitable expansions and 
renormalization-group treatments based on scaling ar- 
guments. These give universal power-law scaling expo- 



nents for the ( "critical" ) behaviour of the system close 
to the phase transition point (critical point) (Stanley, 
1971; Domb and Green, 1972-1976; Ma, 1976; Hohen- 
berg and Halperin, 1977; Stinchcombe, 1983; Domb 
and Lebowitz, 1983-2000; Voss, 1985; Schuster, 1988). 
In non-equilibrium systems, one can frequently find 
self-organized criticality (Bak et al, 1987, 1988; Bak, 
1996) and fractal properties (Mandelbrot, 1983; Fam- 
ily and Vicsek, 1991; Vicsek, 1992; Barabasi and Stan- 
ley, 1995). Self-organized criticality (SOC) means that 
the respective system drives itself to the critical point, 
which is normally characterized by long-range corre- 
lations and scale-free power laws in analogy to ther- 
modynamics. 

In the considered many-particle systems, there is 
often also a certain lower dimension, below which 
the system is always disordered, because of the small 
number of neighbors. For example, it is known that 
neither ferromagnetic nor antiferromagnetic order is 
possible in one-dimensional equilibrium spin systems 
(Mcrmin and Wagner, 1966). The discovery that 
the situation can be very different in driven, non- 
equilibrium spin systems, has initiated an intense re- 
search activity regarding the phase transitions in one- 
dimensional driven diffusive systems far from equilib- 
rium. It turned out that many properties of equilib- 
rium systems can be generalized to non-equilibrium 
systems, but others not. 

By means of particle hopping models, it was possi- 
ble to gain a better understanding of directed percola- 
tion (Domany and Kinzel, 1984), of spontaneous struc- 
ture formation (Vilfan et al., 1994), of spontaneous 
symmetry breaking (Evans et al, 1995), of the rough- 
ening transition in certain growth processes (Alon 
et al., 1996), of the non-equilibrium wetting transi- 
tion (Hinrichscn et al, 1997), and of phase separa- 
tion (Evans et al, 1998; Hclbing et al, 1999b). Nev- 
ertheless, for these non-equilibrium transitions there 
is still no general theory available which would be of 
comparable generality as equilibrium thermodynam- 
ics or statistical physics. For further reading, I recom- 
mend the books by Spohn (1991), Schmittmann and 
Zia (1995, 1998), Derrida and Evans (1997), Liggett 
(1999), and Schutz (2000). 

Note that the above mentioned models usually as- 
sume a random sequential (asynchronous) update, i.e. 
the state of each particle is randomly updated after an 
exponentially distributed waiting time At. However, 
it is known that this is not realistic for traffic systems, 
which rather require a parallel (synchronous) update 
of the vehicle locations (Schreckenberg et al, 1995), as 
in other flow problems. To reflect the driver reaction 
to the corresponding change of the traffic situation, it 
is common to update the vehicle speeds synchronously 
as well. However, this is debatable, since some mod- 
els of spatio-temporal interactions in social systems 
show artefacts, if updated in parallel (Hubcrman and 
Glance, 1993). For example, a parallel update excludes 
the dynamical reachability of certain states (Schad- 



Schneider and Schreckenberg, 1998), which are called 
paradisical or Garden-of-Eden states (Moore, 1962). 



1. The asymmetric simple exclusion process (ASEP) 

The only particle hopping model I will discuss here 
is the asymmetric simple exclusion process. Particu- 
larly relevant contributions to this model go back to 
Spitzer (1970), Liggett (1975), Domany and Kinzel 
(1984), Katz et al. (1984), Liggett (1985), Krug 

(1991) , Dcrrida et al. (1992), Janowski and Lebowitz 

(1992) , Dcrrida et al. (1993), Ferrari (1993), Schutz 
and Domany (1993), Stinchcombe and Schutz (1995a, 
b), Kolomeisky et aZ.(1998), and Schutz (1998). For 
overviews see Spohn (1991), Schmittmann and Zia 
(1995, 1998), Liggett (1999), and Schutz (2000). In 
the following, I will restrict to the totally asymmetric 
simple exclusion process (TASEP). This model is de- 
fined by L sites j of a one-dimensional lattice, which 
can be either empty (corresponding to the occupation 
number rij — 0) or occupied by one particle (rij = 1). 
The particle locations and occupation numbers are up- 
dated every time step At. When updated, a particle 
in cell j hops to the right neighboring cell (j + 1) with 
probability q, if this is empty, otherwise it stays in 
cell j. The total rate of moving to the right is, there- 
fore, given by qrij(l — n J+ i)/At. The boundaries are 
characterized as follows: A particle enters the system 
at the left-most cell j = 1 with probability qo, if this 
is empty. Moreover, a particle in the right-most cell 
j = L leaves the system with probability q^. This 
corresponds to particle reservoirs at the boundaries 
which can be described by constant occupation prob- 
abilities no = qo/q and n^+i = (1 — Ql/q)- Although 
it takes enourmous efforts, the stationary states and 
even the dynamics of the TASEP can be analytically 
determined. For this, one has to solve the correspond- 
ing master equation (see Sees. III.C.3 and III.D.l). 



F. Active Brownian particles 

Like usual Brownian particles, active Brownian 
particles perform a random walk as well. However, 
they arc not only reactive to an external poten- 
tial U(x,t), but they are driven by an internal en- 
ergy reservoir ("pumped particles") or can actively 
change the potential U(x,t) while moving ("active 
walkers"). Therefore, I have recently suggested to call 
them Brownian agents. For an overview see Schweitzer 
(2001). 



1. Pumped Brownian particles 

Schweitzer et al. (1998b) suggest a model describing 
how self-driven particles may take up and consume 
the internal energy behind their driving force. Their 
model corresponds to Eq. (2) with the specifications 



F a f) — (i.e. no direct interactions), but F (x a (t)) is 
replaced by the expression — VJ7(x a ) + jE a (t)v a (t), 
where the first term is an external potential force and 
the last term an internal driving force (7 being the 
friction coefficient). The dimensionless (scaled) energy 
reservoir E a {t) is assumed to follow the equation 



dE a (t) 
alt 



Qo(x a (t)) -{c + d [v a (t)] 2 }E a (t) . (5) 



Herein, Qo(x) reflects the exploitation rate of energy 
resources, while the last term describes the consump- 
tion rate of energy, which grows quadratically with 
the speed (c and d being suitable parameter values). 
If the relaxation of this energy equation is fast, we can 
approximate the driving term by 



7 Qo(x a (t))v a (t) 
c + d[v a (t)} 2 



so that the internal driving and the dissipative fric- 
tion together can be represented by an active fric- 
tion coefficient j'(v) = j[l — Qo/(c + dv 2 )}. More- 
over, the driving direction e° a (t) is given by the nor- 
malized actual velocity e a (t) = v a (t)/\\v a (t)\\, and 
the desired velocity i>°(t) depends on the speed as 
well: w°(t) = Q \\v a (t)\\/{c + d[v a (t)} 2 } > 0. Notice 
that particle a takes up energy only when it moves 
with some finite speed v a (t) = \\v a (t)\\ 7^ (i.e. ex- 
ploits its environment). Therefore, in the absence of 
a potential U(x,t), we find the stationary solution 
v a — v a = 0- However, this is only stable under the 
condition jQ < cm a . Otherwise, particle a will spon- 
taneously start to move with average speed 



r a jQo 
m a d 



1/2 



(6) 



Depending on Qo(x) and U(x), this spontaneous mo- 
tion displays interesting dynamics such as limit cy- 
cles, deterministic chaos, or intermittency (Schweitzer 
et al., 1998b; see also Chen, 1997). It also allows the 
particles to climb potential hills, e.g. in a periodic 
ratchet potential U(x) (Schweitzer et al, 2000). 



2. Dissipative Toda and Morse chains 

Toda and Morse chains are particles coupled to a 
heat bath and moving on a ring with particular asym- 
metrical springs among neighbors, which are described 
by non-linear Toda or Morse potentials U'(xp — x a ) 
(Bolterauer and Opper, 1981; Toda, 1983; Jenssen, 
1991; Ebeling and Jenssen, 1992; Dunkel et al, 2001). 
Assuming pumped particles with an active friction co- 
efficient Y(v) similar to the previous paragraph, Ebel- 
ing et al. (2000) have found interesting dynamical pat- 
terns for overcritical pumping. This includes uniform 
rotations, one- and multiple soliton-like excitations, 
and relative oscillations. For Morse potentials, one 
also observes clustering effects which partly remind 
of jamming (Dunkel et al, 2001). 



3. Active walker models 

By modifying their environment locally, active walk- 
ers have indirect interactions with each other, which 
may lead to the formation of global structures. The 
versatility of this concept for the description of phys- 
ical, chemical, biological, and socio-economic systems 
is discussed by Schweitzer (2001). For example, a sim- 
ple model is given by m a dv a /dt = {—VU(x a ,t) — 
^v a {t) + Ca(t)] and changes of the environmental po- 
tential U(x,t) according to 

9U ^ ^ = -J^B S(x-x a (t))-CU{x, t)+D AU(x, t). 

a 

Herein, C and D are constants, and S(x — x a ) de- 
notes Dirac's delta function, giving contributions only 
at x — x a . Therefore, the first term on the right hand 
side reflects particles a leaving attractive markings at 
their respective locations x a . The last term describes 
a diffusion of the field U(x,t) (e.g., a chemical one), 
and the previous term its decay. 

The result of the dynamics is a spatial agglomer- 
ation process. While in the first stage, one finds the 
emergence of localized clusters at random locations, 
the clusters eventually merge with each other, result- 
ing in one big cluster. The growth and competition of 
these clusters can be described by an Eig en- Fisher-like 
selection equation (Schweitzer and Schimansky-Geier, 
1994; Fisher, 1930; Eigen, 1971; Eigen and Schuster, 
1979). 

4. Pattern formation of bacterial colonies 

Already in 1994, Eshel Ben- Jacob et al. have pro- 
posed a so-called communicating walker model for 
the beautiful pattern formation in certain bacterial 
colonies (for reviews see Ben- Jacob, 1997; Ben- Jacob 
et al., 2000). This model takes into account the con- 
sumption, diffusion, and depletion of food, the multi- 
plication under good growth conditions, and the tran- 
sition to an immobile spore state under starvation, 
as well as the effect of chemotaxis, i.e. the attrac- 
tive or repulsive reaction to self-produced chemical 
fields. The model allows one to reproduce the observed 
growth patterns of bacterial colonies as a function of 
the nutrient concentration and the agar concentration 
detcrminining how hard the bacteria can move. The 
same bacterial patterns have recently been reproduced 
by means of a macroscopic reaction-diffusion model 
(Golding et al., 1998). The formation of the observed 
dendritic structures is due to a so-called diffusion in- 
stability (Ben- Jacob, 1993, 1997). 

5. Trail formation by animals and pedestrians 

Another example is the observed formation of trunk 
trail systems by certain ant species, which is also based 



on chemotaxis. It can be modelled by two kinds of 
chemical markings, one of which is dropped during 
the search for food, while the other is dropped on the 
way back to the nest of the ant colony (Schweitzer 
et al., 1997). Other models for ant-like swarm forma- 
tion have been developed by Deneubourg et al. (1989) 
and Bonabeau (1996). 

Further kinds of active walker models have been 
proposed for the self-organization of trail systems by 
pedestrians or animals, where the markings corre- 
spond to "footprints" or other kinds of modifications 
of the ground which make it more comfortable to 
move (Helbing et al, 1997a, b; Helbing, 1998c). Inter- 
estingly enough, the model yields fair solutions and 
optimal compromises between short and comfortable 
ways, see Fig. 1 (Helbing, 1998c; Helbing and Vic- 
sek, 1999). The corresponding computer simulations 
are, therefore, a valuable tool for developing optimized 
pedestrian facilities and way systems. 

FIG. 1. Schematic representation of a human trail sys- 
tem (black solid lines) evolving between four entry points 
and destinations (full circles) on an initially homogeneous 
ground (after Helbing et al, 1997a, b, 2001c; Helbing, 
1998c; Helbing and Vicsek, 1999). When the frequency 
of trail usage is small, the direct way system (consisting 
of the four ways along the edges and the two diagonal 
connections) is too long to be maintained in competition 
with the regeneration of the vegetation. Here, by bundling 
of trails, the frequency of usage becomes large enough to 
support the depicted trail system. It corresponds to the 
optimal compromise between the diagonal ways and the 
ways along the edges, supplying maximum walking com- 
fort at a minimal detour of 22% for everyone, which is a 
fair solution. 



G. Vehicle and pedestrian traffic 

The next sections will mainly focus on unidirec- 
tional freeway traffic, but bidirectional traffic and 
two-dimensional city traffic are shortly sketched as 
well. I will start with an overview over the most 
important empirical findings. Afterwards, I will con- 
centrate on the different modelling approaches and 
discuss the properties of these models, starting with 
one-lane traffic of identical vehicles and then adding 
more and more details, including heterogeneous traf- 
fic on multi-lane roads. Finally, I will discuss two- 
dimensional pedestrian traffic and three-dimensional 
flocks of birds. 

II. EMPIRICAL FINDINGS FOR FREEWAY TRAFFIC 

As physics deals with phenomena in our physical 
world, a physical theory must face the comparison 
with empirical data. A good theory should not only 
reproduce all the empirically known properties of a 
system, but it should also make predictions allowing 
us to verify or falsify the theory. To gain on overview of 



data analyses, I recommend to read the review articles 
and books by Gcrlough and Hubcr (1975), Koshi et al. 
(1983), the Transportation Research Board (1985), 
May (1990), Daganzo (1997a, 1999a, b), Bovy (1998), 
and Kerner (1998b, c, 1999a, b, 2000a, b). 

A. Measurement techniques 

Probably the most refined technique to gain em- 
pirical data is based on aerial photography, allowing 
us to track the trajectories of many interacting vehi- 
cles and even their lane changing maneuvers (Treit- 
erer and Taylor, 1966; Treiterer and Myers, 1974). A 
method suitable for experimental investigations are 
car-following data. Depending on the equipment of 
the cars, it is possible to determine the location and 
speed, possibly the acceleration and clearance (netto 
distance), and sometimes even lane changing maneu- 
vers of the equipped car or of the respective following 
vehicle (see, for example, Koshi et al., 1983; Bleile, 
1997a, b, 1999). 

However, most data are obtained by detectors lo- 
cated at certain cross sections x of the freeway. For 
example, single induction-loop detectors measure the 
numbers AiV of crossing vehicles a during a certain 
sampling interval AT as well as the times t° a and 
t\ when a vehicle a reaches and leaves the detec- 
tor. This facilitates to determine the time headways 
(brutto time separations) 

and the time clearances (netto time separations or 
time gaps) t° a — t x a _ x including their respective dis- 
tributions, as well as the vehicle flow 

Q(x,t) = ^ (7) 

and the time occupancy 0(x,t) — J2 a (^a ~~ t^)/AT 
(where ao < a < ao + AN and ao denotes the last ve- 
hicle before the sampling interval begins). The newer 
double induction-loop detectors additionally measure 
the vehicle velocities v a and the vehicle lengths l ai al- 
lowing us to estimate the headways (brutto distances) 
d a = v a At a (assuming constant vehicle speeds) and 
the clearances (netto distances) s a = (d a — l a -i). The 
velocity measurements are normally used to obtain 
the (arithmetic) average velocity 

ao+AN 

V(x,t) = (v a ) = — v « ( 8 ) 

but the velocity variance 

0(x, t) = ([v a - (v a )} 2 ) = (K) 2 ) - (v a ) 2 (9) 

and the local velocity distribution may be deter- 
mined as well. The vehicle density p(x, t) is often 



calculated via the fluid-dynamic formula (46), i.e. 
p{x,t) = Q(x,t)/V(x,t). Another method is to re- 
late the density with the time occupancy, as in the 
formula p(x, t) = 0(x, t)/[L(x, t) +Ljj], where L(x, t) 
is the average vehicle length during the measurement 
interval and Ld the detector length (i.e. the loop ex- 
tension) in driving direction (May, 1990). 

The problem of the above measurement methods is 
that the velocity distribution measured over a time 
interval AT differs from the one measured on a free- 
way section of length AX. In other words, temporal 
and spatial averaging yield different results, since fast 
vehicles cross a section of the freeway more frequently 
than slow ones, which is intuitive for a circular multi- 
lane freeway with fast and slow lanes. The problem 
of determining the empirical density via the formula 
p = Q/V is that it mixes a temporal average (the flow) 
with a spatial one (the density). This can be com- 
pensated for by the harmonic average velocity V(x, t) 
defined by 

vh - (i) • (10) 

giving a greater weigth to small velocities (Gcrlough 
and Huber, 1975; Leutzbach, 1988). Using formula 
(10) instead of the commonly applied formula (8) to- 
gether with the relation p = Q/V results in simi- 
larly shaped velocity-density relations, but the density 
range is much higher (see Fig. 2). Moreover, the ve- 
locity and flow values at high densities are somewhat 
lower. The disadvantage of Eq. (10) is its great sensi- 
tivity to errors in the measurement of small velocities 

FIG. 2. Empirical velocity-density relations for different 
definitions of the average velocity. Symbols represent av- 
erages of 1-minute data determined via the harmonic ve- 
locity formula V = l/(l/v a ), while the solid line is a fit 
function to the velocity averages determined via the con- 
ventionally applied arithmetic formula V = (v a ). (After 
Helbing, 1997a.) 

B. Fundamental diagram and hysteresis 

Functional relations between the vehicle flow 
Q(x,t), the average velocity V(x, t), and the vehicle 
density p(x,t) or occupancy 0(x,t) have been mea- 
sured for decades, beginning with Greenshields (1935), 
who found a linear relationship. The name "funda- 
mental diagram" is mostly used for some fit function 

Q e (p) = pV e (p) (11) 

of the empirical flow-density relation, where V e (p) 
stands for the fitted empirical velocity-density re- 
lation (see Fig. 2), which is monotonically falling: 
dV e (p)/dp < 0. Most measurements confirm the fol- 
lowing features: 

(i) At low densities p, there is a clear and quasi one- 
dimensional relationship between the vehicle flow and 



the density. It starts almost linearly and is bent down- 
wards. The slope at low densities corresponds to the 
average free velocity Vq, which can be sustained at fi- 
nite densities as long as there are sufficient possibilities 
for overtaking. That is, the velocity-density relation of 
a multi-lane road starts horizontally while it tends to 
have a negative slope on a one-lane road. 

(ii) With growing density, the velocity decreases 
monotonically, and it vanishes together with the flow 
at some jam density Pjam, which is hard to determine 
because of the above mentioned measurement prob- 
lems. Estimates by different researchers from various 
countries reach from 120 to 200 vehicles per kilometer 
and lane, but values between 140 and 160 vehicles per 
kilometer are probably most realistic. 

(iii) The vehicle flow has one maximum Q max at 
medium densities. 

(iv) The empirical flow-density relation is disconti- 
nous and looks comparable to a mirror image of the 
Greek letter lambda. The two branches of this reverse 
lambda are used to define free low-density and con- 
gested high-density traffic (see, e.g., Koshi et al, 1983; 
Hall et al, 1986; Neubert et al., 1999b; Kerncr, 2000a; 
cf. also Edie and Foote, 1958). In congested traffic, the 
average vehicle velocity is significantly lower than in 
free traffic. Therefore, the free and congested part of 
the flow-density relation can be approximately sep- 
arated by a linear flow-density relation pV scp , where 
Vsep is the average velocity below which traffic is char- 
acterized as congested. Around this line, the density 
of data points is reduced (see Fig. 3a). 

FIG. 3. Empirical time series of the hysteretic breakdown 
of traffic flow plotted (a) as a function of the density and 
(b) as a function of the time. High-flow states are only ob- 
served shortly before the breakdown of traffic. Note that, 
immediately after the breakdown, the flow drops tem- 
porarily below the typical bottleneck flow (see Persaud 
et al, 1998; Cassidy and Bertini, 1999). While free traffic 
is characterized by a quasi one-dimensional flow-density 
relationship, in the congested regime the flow-density data 
are widely scattered. The free and congested traffic regimes 
can be separated by a line pV Be p, where we have chosen 
Veep = 70 km/h, here. 

The reverse-lambda-like structure of the flow- 
density data implies several things: First, we have a 
discontinuity at some critical density p cr (Edie, 1961; 
Treiterer and Myers, 1974; Cedcr and May 1976; 
Payne, 1984), which has animated researchers (Hall, 
1987; Dillon and Hall, 1987; Persaud and Hall, 1989; 
Gilchrist and Hall, 1989) to relate traffic dynamics 
with catastrophe theory (Thorn, 1975; Zeeman, 1977). 
Second, there is a certain density region p c \ < p < p C 2 
with p C 2 = Pa, in which we can have either free or 
congested traffic, which points to hysteresis (Treit- 
erer and Myers, 1974). The tip of the lambda corre- 
sponds to high-flow states of free traffic (in the pass- 
ing lanes), which can last for many minutes (see, e.g., 
Cassidy and Bertini, 1999). However, these are not 



stable, since it is only a matter of time until a tran- 
sition to the lower, congested branch of the lambda 
takes place (Persaud et al, 1998; see also Elefteriadou 
et al., 1995). The transition probability from a free 
to a congested state within a long time interval is 
below the flow Q c \ = Q e (p c i) and 1 above the flow 
Qc2 = Qe{pc2)- Between Q c \ and Q C 2, it is monotoni- 
cally increasing, see Fig. 4 (Persaud et al., 1998). The 
transition from congested to free traffic does not go 
back to high-flow states, but to flows Q e (p) belonging 
to densities p < p C 2, typically around p cl . The high- 
flow states in the density range p c \ < p < p C 2 are 
only metastable (Kcrner and Rchborn, 1996b, 1997; 
Kerner, 1998b, 1999a, b, c, 2000a, b). 

FIG. 4. Probability of breakdown of free traffic as a func- 
tion of the flow for different waiting times (after Persaud 
et al, 1998). 

(v) The flow-density data in the congested part are 
widely scattered within a two-dimensional area (see, 
e.g., Koshi et al, 1983; Hall et al, 1986; Kiihnc, 
1991b). This complexity in traffic flow is usually inter- 
preted as effect of fluctuations, of jam formation, or 
of an instability of vehicle dynamics. The scattering 
is reduced by increasing the sampling interval AT of 
the data averaging (Leutzbach, 1988). 

(vi) By removing the data belonging to wide moving 
jams (see Fig. 20), Kerner and Rchborn (1996b) could 
demonstrate that the remaining data of congested 
traffic data still display a wide and two-dimensional 
scattering (see Fig. 3a), thereby questioning the ap- 
plicability of a fundamental diagram and defining the 
state of "synchronized flow" ("synchronized" because 
of the typical synchronization between lanes in con- 
gested traffic, see Fig. 15a, and "flow" because of flow- 
ing in contrast to standing traffic in fully developed 
jams). Therefore, Kerner and Rehborn classify three 
traffic phases: 

- free flow, 

- synchronized flow (see Sees. II. C, II. D, and II. E. 2 
for details), and 

- wide moving jams (i.e. moving jams whose width in 
longitudinal direction is considerably higher than the 
width of the jam fronts). 

In contrast to synchronized flow, a wide moving jam 
propagates either through free or any kind of sychro- 
nized flow and through bottlenecks (e.g. on- and off- 
ramps), keeping the propagation velocity of the jam's 
downstream front (Kerner, 2000a, b). In contrast to 
wide moving jams, after synchronized flow has oc- 
curred upstream of an on-ramp, the downstream front 
of synchronized flow is fixed at the on-ramp. In an ini- 
tially free flow, two types of transitions are observed: 
either to synchronized flow or to (a) wide moving 
jam(s). Both of them appear to be first-order phase 
transitions accompanied by different hysteresis and 
nuclcation effects (Kerner and Rehborn, 1997; Kerner, 



1998a, 1999a, 2000c). 

(vii) The slopes S£ T (t) = [Q(x,t + AT) - 
Q(x,t)]/[p(x,t + AT) — p(x, t)] of the lines connect- 
ing successive data points are always positive in free 
traffic. In synchronized flow, however, they erratically 
take on positive and negative values, characterizing a 
complex spatio-temporal dynamics (Kerner and Reh- 
born, 1996b). This erratic behavior is quantitatively 
characterized by a weak cross-correlation between flow 
and density (Neubert et ai, 1999b). Banks (1999) 
showed that it could be interpreted as a result of ran- 
dom variations in the time clearances (partly due to 
acceleration and deceleration maneuvers in unstable 
traffic flow). These variations are, in fact, very large 
(cf. Fig. 6). Banks points out that, if drivers would 
keep a save clearance s*(v) = s' + Tv in congested 
traffic, the flow would grow with decreasing (effective) 
safe time clearance T according to 

Qe(p) = ± (l ~ -?-) , (12) 

where s' denotes the minimum front-bumper-to-rear- 
bumper distance and Pj am = l/l' = 1/(1 + s') the 
jam distance with the (average) effective vehicle length 
I' = (l+s 1 ). Simplifying his argument, in the congested 
regime positive slopes can result from a reduction in 
the safe time clearance T with growing density, while 
negative slopes normally correspond to a reduction in 
the speed. Hence, the slopes S^ T (t) do not necessarily 
have the meaning of a speed of wave propagation. 

(viii) The flow-density data depend on the measure- 
ment cross section. While the congested branch is very 
pronounced upstream of a bottleneck, it is virtually 
not present far downstream of it. However, immedi- 
ately downstream of a bottleneck, one finds a pos- 
itively sloped congested area slightly below the free 
branch ( "homogeneous-in- speed states"). It looks sim- 
ilar to the upper part of the free branch, but with 
a somewhat lower desired velocity. I suggest that this 
may indicate a transition to free traffic in the course of 
the road (see Figs. 5 and 43b) and that this "recover- 
ing traffic" bears already some signatures of free traffic 
(cf. also Persaud and Hurdle, 1988; Hall et al, 1992). 
For example, on-ramp flows just add to the flows on 
the freeway, which may produce states with high flows 
(Kerner, 2000b). Some observations, however, ques- 
tion the interpretation of homogeneous- in-speed states 
as recovering traffic, as they can extend over stretches 
of more than 3 kilometers (Kerner, 1999a). 

FIG. 5. Homogeneous-in-speed states scatter around a 
line through the origin, indicating constant speeds. The 
data are typical for "recovering traffic" at cross sections 
located downstream of a bottleneck, where congested traf- 
fic relaxes to free traffic. 



C. Time headways, headways, and velocities 

Time headways show an astonishing individual 
variation, supporting Banks' theory discussed above. 
When distinguished for different density ranges, time- 
headway distributions have an interesting property. 
It turns out that the distribution becomes sharply 
peaked around approximately 1 s for densities close to 
congestion (Smulders, 1990; Hoogendoorn and Bovy, 
1998c; Neubert et ai, 1999b). For lower and higher 
densities, the distribution is considerably broader, 
see Fig. 6 (Tilch and Helbing, 2000). This indicates 
that congestion is an overcritical phenomenon occur- 
ing when the maximum efficiency in terms of time 
headways (related to high-flow states) cannot be sup- 
ported any longer. One may conclude that the effec- 
tive time clearance has, then, reached the safe time 
clearance T a . A further increase in the density forces 
to reduce the speed, which automatically increases the 
time headway At a = (T a +l' a /v a ), even if the effective 
clearance remains T a (cf. Banks, 1991b). 

The distributions of vehicle headways d a = v a At a 
can be determined as well, but it is always surprisingly 
broad (see Fig. 7). Therefore, Bleile (1999) suggests to 
consider this by an individual distance scaling of the 
velocity-clearance relation. 

FIG. 6. Empirical time-headway distributions for differ- 
ent vehicle densities and traffic regimes, (a) Due to the 
larger proportion of long vehicles (see Sec. II. F), the time- 
headway distribution is broader on the right lane com- 
pared to the left one. (b) Before the breakdown of free 
traffic flow, the time-headway distribution is particularly 
narrow, afterwards it becomes wide. (After Tilch and Hel- 
bing, 2000.) 



FIG. 7. The distribution of vehicle distances d is rather 
broad at all densities, although the average vehicle dis- 
tance decreases with the inverse density 1/p (after Tilch 
and Helbing, 2000). 

What are the reasons and consequences of the wide 
scattering of time headways and clearances, apart 
from driver-dependent preferences? By detailed in- 
vestigations of clearance-over-speed relations, Koshi 
et al. (1983) have observed that vehicles keep larger 
distances in real congested traffic than in free traffic 
or undisturbed congested traffic under experimental 
conditions (cf. Fig. 8). 

FIG. 8. Average vehicle distances d as a function of the 
individual vehicle velocity v, separately for time intervals 
with free and congested traffic (according to the classifi- 
cation introduced in Fig. 3a). Vehicles keep increased dis- 
tances in congested traffic flow ( ) upstream of a bot- 
tleneck (a), but not in an undisturbed section far enough 
away from it (b). The broken and solid lines are fit curves 
in order to guide the eyes. (After Tilch and Helbing, 2000; 
note the related studies by Koshi et al, 1983; Dijker et al, 
1998.) 



Plotting the speed over the vehicle distance, one can 
find a density-dependent reduction of the speed (see 
Fig. 9), which has been called "frustration effect". 

FIG. 9. Average vehicle velocities v as a function of the 
headways d = v At for various density regimes, separately 
for free traffic (f) and congested traffic (c) on (a) the left 
and (b) the right lane for a cross section upstream of a bot- 
tleneck. (From Tilch and Helbing, 2000; note the related 
study by Neubert et al, 1999b.) Headways of less than 5 m 
are due to measurement errors. The dependence of the v(d) 
relation on the vehicle density and traffic regime is called 
frustration effect (see also Brilon and Ponzlet, 1996). 

In my opinion, this is partly an effect of a reduc- 
tion in the vehicle speed with decreasing distance due 
to safety requirements, combined with the wide dis- 
tance scattering illustrated in Fig. 7 (Tilch and Hel- 
bing, 2000). Another relevant factor is the influence 
of the relative velocity Av on driver behavior (Blcilc, 
1997a, 1999). The average clearances of vehicles which 
are driving faster or slower than the respective leading 
vehicles are naturally increased compared to vehicles 
driving at the same speed, see Fig. 10 (Neubert et al, 
1999b; Tilch, 2001). In summary, driver behavior is 
not only influenced by the clearance to the next car, 
but also by the relative velocity and the driver veloc- 
ity (Koshi et al, 1983; Bleile, 1997, 1999; Dijker et al, 
1998; Neubert et al, 1999b). 

FIG. 10. Empirical clearance s = (d — V) as a function of 
the relative velocity Av, for free traffic (f) and congested 
traffic (c) in different density regimes. The clearance is 
minimal for identical vehicle velocities. (From Tilch 2001, 
after Neubert et al., 1999b.) 

Note that the relative velocity in congested traffic is 
oscillating because of the instability of traffic flow, see 
Fig. 11 (Hoefs, 1972; Helbing and Tilch, 1998). As a 
consequence, the relative velocity variance 9, as com- 
pared to the average velocity V, is considerably higher 
in congested than in free traffic (see Figs. 12 and 17). I 
suggest that this, together with the factors mentioned 
in the previous paragraph, may explain the empirical 
relations depicted in Fig. 9. If this is confirmed, one 
can drop the assumption of a frustration of drivers in 
congested traffic. 

FIG. 11. Measured oscillations of the clearance s = (d—l) 
and of the relative velocity Av around Av — (cf. Hoefs, 
1972) indicate an instability in car-following behavior (af- 
ter Helbing and Tilch, 1998). 



FIG. 12. Empirical standard deviation U 8(t) of vehi- 
cle velocities divided by the average velocity V(t) as a 
function of time (from Tilch, 2001). The relative varia- 
tion yj 0{t) /V {t) is particularly large during the rush hour, 
when traffic is congested. This reflects the oscillations in 
the following behavior due to unstable traffic flow (cf. 
Fig. 11) and implies that measurements of average ve- 
locities V have a small reliability in the congested traffic 
regime. 

Measurements of velocity distributions for traffic 
with small truck fractions are in good agreement with 
a Gaussian distribution, see Fig. 13 (Pampel, 1955; 
Leong, 1968; May, 1990). However, if the sampling 
intervals AT are taken too large, one may also ob- 
serve bimodal distributions (Phillips, 1977; Kiihne, 
1984a, b), reflecting a transition from free to congested 
traffic. As expected for a Gaussian distribution, 1- 
minute averages of single- vehicle data for the skewness 
(K - V) 3 }/9 3 / 2 and the kurtosis [({v a - V) 4 }/9 2 - 3] 
are scattering around the zero line (Helbing, 1997e, a, 
c, 1998a). 

FIG. 13. Comparison of empirical velocity distributions 
at different vehicle densities p ( — ) with frequency poly- 
gons of grouped normal distributions having the same 

mean value and variance ( ). A significant deviation is 

only observed for p = 40 veh. /km/lane, probably because 
of the instability of traffic flow. (After Helbing, 1997a, c, 
e, 1998a.) 

Like the velocity distribution itself, the higher ve- 
locity moments are sensitive to the choice of the sam- 
pling interval AT. Large values of AT may lead to 
peaks in the variance 9 when the average velocity V 
changes abruptly. The additional contribution can be 
estimated as [dV{x, t)/dt] 2 (AT) 2 /4, where the over- 
bar indicates a time average over the sampling interval 
(Helbing, 1997a, b, c). 1-minute data of the variance 
9 are more or less monotonically decreasing with the 
vehicle density (see Fig. 14), while 5-minute data have 
maxima in the density range between 30 and 40 ve- 
hicles per kilometer and lane, indicating particulary 
unstable traffic (Helbing, 1997a). 

FIG. 14. Density-dependent standard deviation \/6 of 
vehicle velocities in a single lane of the freeway (•: 1-minute 
data, O: average values for a given density p). The veloc- 
ity variance 9 is more or less monotonically falling with 
increasing density. The greater scattering above 30 vehi- 
cles per kilometer indicates unstable traffic flows. (After 
Helbing, 1997a, b, e.) 

Note that the average 9 = J2i P$i I (Hi Pi) 01 the 
velocity variances 9i in the different lanes i is by an 
amount of ((V* — V) 2 ) lower than the velocity vari- 
ance evaluated over all lanes, in particular at low den- 
sities pi. This is due to the different average velocities 
Vi in the neighboring lanes. For a two-lane freeway, 
the difference in average speeds decreases almost lin- 



early up to a density of 35 to 40 vehicles per kilo- 
meter, while it fluctuates around zero at higher den- 
sities, see Fig. 15a (Helbing, 1997a, b; Hclbing et al, 
2001b). This reflects a synchronization of the velocities 
Vi{x,t) in neighboring lanes in congested traffic, both 
in wide moving jams and synchronized flow (Koshi 
et al, 1983; Kerner and Rchborn, 1996b; for reduced 
speed differences in congested traffic see also Edie and 
Foote, 1958; Forbes et al, 1967; Mika et al, 1969). 
However, according to Figs. 15b, c, there is a differ- 
ence in the speeds of cars and long vehicles ( "trucks" ) , 
which indicates that vehicles can still sometimes over- 
take, apart from a small density range around 25 vehi- 
cles per kilometer, where cars move as slow as trucks 
(Helbing and Hubcrman, 1998). However, direct mea- 
surements of the number of lane changes or of overtak- 
ing maneuvers as a function of the macroscopic vari- 
ables are rare (Sparman, 1978; Hall and Lam, 1988; 
Chang and Kao, 1991). It would be particularly in- 
teresting to look at the dependence of lane-changing 
rates on the velocity-difference among lanes. 

FIG. 15. (a) Difference between the average velocity 
Vi(p) in the left lane and Vi(p) in the right lane as a func- 
tion of the lane-averaged vehicle density p (from Helbing 
et al, 2001b; see also Helbing, 1997a, b; Helbing et al, 
2001a). (b), (c) Average velocities of cars and trucks (long 
vehicles) in the left and the right lane as a function of the 
density. The difference of these empirical curves shows a 
minimum around 25 vehicles per kilometer, where cars are 
almost as slow as trucks (after Helbing and Huberman, 
1998; Helbing et al, 2001b; Helbing, 2001). However, at 
higher densities, cars are faster again, which shows that 
there must be overtaking maneuvers in the congested den- 
sity regime. 



D. Correlations 

tn congested traffic, the average velocities in neigh- 
boring lanes are synchronized, see Fig. 16 (Mika et al., 
1969; Koshi et al, 1983; cf. also Kerner and Rch- 
born, 1996b, for the non-linear features of synchro- 
nized flow). The synchronization is probably a result 
of lane changes to the faster lane, until the speed 
difference is equilibrated (for recent results see, e.g., 
Lee et al, 1998; Shvetsov and Helbing, 1999; Nelson, 
2000). This equilibration process leads to a higher ve- 
hicle density in the so-called "fast" lane, which is used 
by less trucks (see, for example, Hall and Lam, 1988; 
Helbing, 1997a, b). Nevertheless, density changes in 
congested neighboring lanes are correlated as well. 
Less understood and, therefore, even more surpris- 
ing is the approximate synchronization of the lane- 
specific velocity variances 8i over the whole density 
range (Helbing, 1997a, b, e). It may be a sign of adap- 
tive driver behavior. 



FIG. 16. Comparison of the temporal evolution of differ- 
ent aggregate (macroscopic) traffic variables on neighbor- 
ing lanes i upstream of a bottleneck (from Helbing, 1997a, 
b, e). (a) The average velocities Vi(t) synchronize in the 
congested traffic regime during the rush hours, and the 
densities pi(t) are varying in a correlated manner. While in 
the left lane, the speed is higher in free traffic, in congested 
traffic it is the density (both of which is a consequence of 
the smaller truck fraction), (b) The standard deviations 
\f9~i(t) of individual vehicle velocities show always a ten- 
dency to be synchronized. If (and only if) the time interval 
of data averaging is large (here it is 5 min), the variance 
shows peaks due to sudden changes in the average velocity 
Vi(t). 

Correlations arc also found between the density 
and flow or average velocity. At small densities, there 
is a strong positive correlation with the flow, which 
is reflected by the almost linear free branch of the 
fundamental diagram. In contrast, at large densities, 
the velocity and density are strongly anticorrelated 
(see Fig. 16) because of the monotonically decreas- 
ing velocity-density relation (see Fig. 2). The average 
velocity and variance are positively correlated, since 
both quantities are related via a positive, but density- 
dependent prefactor, see Fig. 17 (Helbing, 1996b, 
1997a, c, Trcibcr et al, 1999; Shvetsov and Helbing, 
1999; Helbing et al, 2001a, b). 

FIG. 17. The density- dependent variance prefactors 
Ai(pi) = 9i(pi)/[Vi(pi)] 2 (reflecting something like 
squared relative individual velocity variations) have sig- 
nificantly increased values in the congested traffic regime, 
which may be a result of unstable traffic flow (see Figs. 11 
and 12). The empirical data can be approximated by fit 
functions of the form Ai(pi) = Aoi + AAi{l +cxp[— (p; — 
Pd) / Api]}~ , where Acn and Aoi + AAi are the variance 
prefactors for free and congested traffic, respectively, p c i is 
of the order of the critical density for the transition from 
free to congested traffic, and Apt denotes the width of the 
transition. (From Shvetsov and Helbing, 1999.) 

Finally, it is interesting to look at the velocity corre- 
lations among successive cars. Neubert et al (1999b) 
have found a long-range velocity correlation in syn- 
chronized flow, while in free traffic vehicle velocities 
are almost statistically independent. This finding is 
complemented by results from Hclbing et al (2001b), 
see Fig. 18. The observed velocity correlations may 
be interpreted in terms of vehicle platoons (Wagner 
and Peinke, 1997). Such vehicle platoons have been 
tracked by Edie and Foote (1958) and Treiterer and 
Myers (1974). 



FIG. 18. (a) Velocity correlation in synchronized flow as 
a function of the number of consecutive vehicles. The data 
indicate long-range correlations, which are not observed in 
free traffic or for traffic variables such as the time headway 
or headway. (After Neubert et al., 1999b.) (b) Correlation 
between the velocities of successive vehicles in the right 
lane as a function of the density (thin lines), determined 
from single- vehicle data (from Helbing et al, 2001b). As 
expected, the correlation coefficient is about zero at small 
densities. It reaches a maximum at around 20 to 30 vehi- 
cles per kilometer, where the transition from free to con- 
gested traffic occurs. Afterwards, it stays on a high level 
(around 0.65). In the left lane (not displayed), the velocity 
correlation is a little bit higher, probably because of the 
smaller fraction of trucks (long vehicles). The correlation 
coefficient is not sensitive to the type of the leading ve- 
hicle (car or truck), but to that of the following vehicle, 
which is reasonable. The thick line (- • -) indicates that, 
even at small densities of 10 veh. /km/lane, the velocity 
correlation between cars depends strongly on the head- 
way to the next vehicle ahead and becomes almost one for 
very small headways (after Helbing et al, 2001b). This is, 
because different velocities would imply a high danger of 
accidents, then. 

E. Congested traffic 

1. Jams, stop-and-go waves, and power laws 

The phenomenon of stop-and-go waves (start-stop 
waves) has been empirically studied by a lot of au- 
thors, including Edie and Foote (1958), Mika et al. 
(1969), and Koshi et al. (1983). The latter have found 
that the parts of the velocity-profile, which belong 
to the fluent stages of stop-and-go waves, do not sig- 
nificantly depend on the flow (regarding their height 
and length), while the frequency does. Correspond- 
ingly, there is no characteristic frequency of stop-and- 
go traffic, which indicates that we are confronted with 
non-linear waves. The average duration of one wave 
period is normally between 4 and 20 minutes for wide 
traffic jams (see, e.g., Mika et al., 1969; Kiihne, 1987; 
Helbing, 1997a, c, e), and the average wave length be- 
tween 2.5 and 5 km (see, e.g., Kerner, 1998a). Analyz- 
ing the power spectrum (Mika et al, 1969; Koshi et al, 
1983) points to white noise at high wave frequencies 
u> (Helbing, 1997a, c, e), while a power law with 
exponent <f> = 1.4 is found at low frequencies (Musha 
and Higuchi, 1976, 1978). The latter has been inter- 
preted as sign of self-organized criticality (SOC) in the 
formation of traffic jams (Nagel and Herrmann, 1993; 
Nagel and Paczuski, 1995). That is, congested traffic 
would drive itself towards the critical density p cr , re- 
flecting that it tries to re-establish the largest vehicle 
density associated with free flow (see below regarding 
the segregation between free and congested traffic and 
Sec. IV. A. 5 regarding the constants of traffic flow). 

Based on evaluations of aerial photographs, Treit- 
erer (1966, 1974) has shown the existence of "phan- 
tom traffic jams", i.e. the spontaneous formation of 



traffic jams with no obvious reason such as an acci- 
dent or a bottleneck (see Fig. 19). According to Da- 
ganzo (1999a), the breakdown of free traffic "can be 
traced back to a lane change in front of a highly com- 
pressed set of cars" , which shows that there is actu- 
ally a reason for jam formation, but its origin can be 
a rather small disturbance. However, disturbances do 
not always lead to traffic jams, as del Castillo (1996a) 
points out. Even under comparable conditions, some 
perturbations grow and others fade away, which is in 
accordance with the metastability of traffic mentioned 
above (Kcrncr and Konhauscr, 1994; Kerner and Rch- 
born, 1996b, 1997; Kerner, 1999c). 

FIG. 19. Vehicle trajectories obtained by Treiterer and 
Myers (1974) by aerial photography techniques (reproduc- 
tion from Leutzbach, 1988). The location (in units of 100 
ft) is displayed over the time (in seconds). While the slopes 
of the trajectories reflect the individual vehicle velocities, 
their density represents the spatial vehicle density. Corre- 
spondingly, the figure shows the formation of a "phantom 
traffic jam", which stops vehicles temporarily. Note that 
the downstream jam front propagates upstream with con- 
stant velocity. 

While small perturbations in free traffic travel 
downstream with a density-dependent velocity C (p) < 
V (Hillegas et al., 1974), large perturbations propa- 
gate upstream, i.e., against the direction of the vehicle 
flow (Edie and Foote, 1958; Mika et al., 1969). The 
disturbances have been found to propagate without 
spreading (Cassidy and Windover, 1995; Windower, 
1998; Muhoz and Daganzo, 1999; see also Kerner and 
Rchborn, 1996a, and the flow and speed data re- 
ported by Foster, 1962; Cassidy and Bertini, 1999). 
The propagation velocity C in congested traffic seems 
to be roughly comparable with a "natural constant" . 
In each country, it has a typical value in the range 
Co = 15 ± 5 km/h, depending on the accepted safe 
time clearance and average vehicle length (see, e.g., 
Mika et al., 1969; Kerner and Rehborn, 1996a; Cassidy 
and Mauch, 2000). Therefore, fully developed traffic 
jams can move in parallel over a long time periods and 
road sections. Their propagation speed is not even in- 
fluenced by ramps, intersections, or synchronized flow 
upstream of bottlenecks, see Fig. 20 (Kerner and Rch- 
born, 1996a; Kcrncr, 2000a, b). 

FIG. 20. Example of two wide moving traffic jams prop- 
agating in parallel with constant speed through free and 
congested traffic and across three intersections II, 12, and 
13 (from Kerner, 2000a, b; see also Kerner and Rehborn, 
1996a; Kerner, 1998b). 

Wide moving jams are characterized by stable wave 
profiles and by kind of "universal" parameters (Kerner 
and Rehborn, 1996a). Apart from 

(i) the propagation velocity Co , these are 

(ii) the density pj am inside of jams, 

(iii) the average velocity and flow inside of traffic jams, 



both of which are approximately zero, 

(iv) the outflow Q ou t from jams (amounting to about 
2/3 of the maximum flow reached in free traffic, which 
probably depends on the sampling interval AT), 

(v) the density p out downstream of jams, if these are 
propagating through free traffic ( "segregation effect" 
between free and congested traffic). When propagat- 
ing through synchronized flow, the outflow Q ut° = 
Qe(Pout°) °f wide moving jams is given by the den- 
sity Pout° °f the surrounding traffic (Kerner, 1998b). 
The concrete values of the characteristic parameters 
slightly depend on the accepted safe time clearances, 
the average vehicle length, truck fraction, and weather 
conditions (Kerner and Rehborn, 1998a). 

2. Extended congested traffic 

The most common form of congestion is not local- 
ized like a wide moving jam, but spatially extended 
and often persisting over several hours. It is related 
with a capacity drop to a flow, which is (at least in the 
United States) typically 10% or less below the "break- 
down flow" of the previous high- flow states defin- 
ing the theoretically possible capacity (Banks, 1991a; 
Kerner and Rehborn, 1996b, 1998b; Persaud et al, 
1998; Westland, 1998; Cassidy and Bertini, 1999; see 
also May, 1964, for an idea how to exploit this phe- 
nomenon with ramp metering; Persaud, 1986; Banks, 
1989, 1990; Agyemang-Duah and Hall, 1991; Daganzo, 
1996). For statistical reasons, it is not fully satisfac- 
tory to determine capacity drops from maximum flow 
values, as these depend on the sampling interval AT. 
Anyway, more interesting are probably the bottleneck 
flows Qhot after the breakdown of traffic. Because of 
relation (113), these may be more than 30% below 
the maximum free flow Q max (see Fig. 3b). Note that 
bottleneck flows depend, for example, on ramp flows 
Q rmp and may, therefore, vary with the location. The 
congested flow immediately downstream of a bottle- 
neck defines the discharge flow Q ou t > Qbot, which 
appears to be larger than (or equal to) the character- 
istic outflow Q ou t from wide jams: Qout > Qout- 

In extended congested traffic, the velocity drops 
much more than the flow, but it stays also finite. The 
velocity profiles V(x, t) can differ considerably from 
one cross section to another (see Fig. 23). In con- 
trast, the temporal profiles of congested traffic flow 
Q(x,t), when measured at subsequent cross sections 
x of the road, are often just shifted by some time in- 
terval which is varying (see Fig. 21). This may be ex- 
plained with a linear flow-density relation of the form 
(12) with dQ e {p)/dp = Co = -l/(7> jam ) (Hall et al, 
1986; Ozaki, 1993; Hall et al, 1993; Dijkcr et al, 1998; 
Westland, 1998), together with the continuum equa- 
tion for the conservation of the vehicle number (see 
the kinematic wave theory in Sec. III.D.l; Cassidy and 
Mauch, 2000; Smilowitz and Daganzo, 1999; Daganzo, 
1999a). However, a linear flow-density relation in the 



congested regime is questioned by the "pinch effect" 
(see Sec. II.E.3 and Fig. 23). 

FIG. 21. Example of time-dependent traffic flows at 
three subsequent cross sections of a 3 kilometer long free- 
way section without ramps (including a hardly notice- 
able breakdown from free to "synchronized" congested flow 
during the rush hours between 7:30 and 9:30 am, cf. the 
corresponding Fig. 16.) The flow profiles at subsequent 
cross sections are mainly shifted by some time interval, al- 
though fluctuations appear in the congested regime (after 
Helbing, 1997a). 

The extended form of congestion is classical and 
mainly found upstream of bottlenecks, so that it nor- 
mally has a spatially fixed ("pinned") downstream 
front. In contrast, the upstream front is moving 
against the flow direction, if the (dynamic) capacity 
Qbot of the bottleneck is exceeded, but downstream, 
if the traffic volume, i.e. the inflow to the freeway 
section is lower than the capacity. Hence, this spa- 
tially extended form of congestion occurs regularly 
and reproducible during rush hours. Note, however, 
that bottlenecks may have many different origins: on- 
ramps, reductions in the number of lanes, accidents 
(even in opposite lanes because of rubbernecks) , speed 
limits, road works, gradients, curves, bad road condi- 
tions (possibly due to rain, fog, or ice), bad visibility 
(e.g., because of blinding sun), diverges (due to "neg- 
ative" perturbations, see Sec. IV. B, weaving flows by 
vehicles trying to switch to the slow exit lane, or con- 
gestion on the off-ramp, see Daganzo et al., 1999; Law- 
son et al., 1999; Muhoz and Daganzo, 1999; Daganzo, 
1999a). Moving bottlenecks due to slow vehicles are 
possible as well (Gazis and Herman, 1992), leading 
to a forward movement of the respective downstream 
congestion front. Finally, in cases of two subsequent 
inhomogeneities of the road, there are forms of con- 
gested traffic in which both, the upstream and down- 
stream fronts are locally fixed (Trcibcr et al., 2000; 
Lee et al, 2000; Kerner, 2000a). 

While Kerner calls the above extended forms of con- 
gested traffic "synchronized flow" because of the syn- 
chronization among lanes (see Sees. II. C, II. D), Da- 
ganzo (1999b) speaks of 1-pipe flow. Kerner and Reh- 
born (1997) have pointed out that the transition from 
free to extended congested traffic is of hysteretic na- 
ture and looks similar to the first-order phase tran- 
sition of supersaturated vapor to water ("nucleation 
effect"). It is often triggered by a small, but over- 
critical peak in the traffic flow (playing the role of 
a nucleation germ). This perturbation travels down- 
stream in the beginning, but it grows eventually and 
changes its propagation direction and speed, until it 
travels upstream with velocity Co- When the pertur- 
bation reaches the bottleneck, it triggers a breakdown 
of traffic flow to the bottleneck flow Qbot- In summary, 
synchronized flow typically starts to form downstream 
of the bottleneck (cf. Fig. 39b) . 

Kerner (1997, 1998a) has attributed the capacity 



drop to increased time headways due to delays in ac- 
celeration. This has recently been supported by Neu- 
bert et al. (1999b). At a measurement section located 
directly at the downstream front of a congested region, 
they found the appearance of a peak at about 1.8 sec- 
onds in the time-headway distribution, when conges- 
tion set in, see Fig. 22. The reachable saturation flows 
of accelerating traffic are, therefore, Q ut ~ 2000 vehi- 
cles per hour. This is compatible with measurements 
of stopped vehicles accelerating at a traffic light turn- 
ing green (Androsch, 1978) or after an incident (Raub 
and Pfcfcr, 1998). 

FIG. 22. Time-headway distributions for different den- 
sity regimes at a cross section where vehicles accelerate 
out of congested traffic. While freely flowing vehicles use 
to have short time headways around 1 s, the time headway 
of accclarating vehicles is typically around 1.8 s. (After 
Neubert et al, 1999b.) 

Kerner (1998b) points out that the above transition 
to "synchronized" congested flow is, in principle, also 
found on one-lane roads, but then it is not anymore 
of hysteretic nature and not connected with synchro- 
nization. Moreover, he and Rchborn (1996b) distin- 
guish three kinds of synchronized flow: 

(i) Stationary and homogeneous states where both 
the average speed and the flow rate are stationary 
during a relatively long time interval (see, e.g., also 
Hall and Agyemang-Duah, 1991; Persaud et al., 1998; 
Westland, 1998). I will later call these "homogeneous 
congested traffic" (HCT). 

(ii) States where only the average vehicle speed is 
stationary, named "homogeneous-in- speed states" (see 
also Kerner, 1998b; Lee et al, 2000). I interpret this 
state as "recovering traffic", since it bears several sig- 
natures of free traffic (cf. Sees. II. B and IV. C. 2). 

(iii) Non-stationary and non-homogeneous states (see 
also Kerner, 1998b; Cassidy and Bertini, 1999; Treiber 
et al., 2000). For these, I will also use the term "os- 
cillating congested traffic" (OCT). 

Not all of these states are long-lived, since there 
are often spatio-temporal sequences of these different 
types of synchronized flow, which indicates continuous 
(second- order) transitions. However, at least states (i) 
and (iii) are characterized by a wide, two-dimensional 
scattering of the flow-density data, i.e. an increase in 
the flow can be either related with an increase or with 
a decrease in the density, in contrast to free flow (see 
Sec. II.B). 

3. Pinch effect 

Recently, the spontaneous appearance of stop-and- 
go traffic has been questioned by Kerner (1998a) and 
Daganzo et al. (1999). In his empirical investigations, 
Kerner (1998a) finds that jams can be born from ex- 
tended congested traffic, which presupposes the pre- 
vious transition from free to synchronized flow. The 



alternative mechanism for jam formation is as follows 
(see Fig. 23): Upstream of a section with homoge- 
neous congested traffic close to a bottleneck, there is 
a so-called "pinch region" characterized by the spon- 
taneous birth of small narrow density clusters, which 
are growing while they travel further upstream. Wide 
moving jams are eventually formed by the merging 
or disappearance of narrow jams, which are said to 
move faster upstream than wide jams. Once formed, 
the wide jams seem to suppress the occurence of new 
narrow jams in between. Similar findings were re- 
ported by Koshi et al. (1983), who observed that "rip- 
ples of speed grow larger in terms of both height and 
length of the waves as they propagate upstream" . Da- 
ganzo (1999b) notes as well that "oscillations exist 
in the 1-pipe regime and that these oscillations may 
grow in amplitude as one moves upstream from an 
active bottleneck" , which he interprets as a pumping 
e/Ject based on ramp flows. The original interpretation 
by Koshi et al. suggests a non-linear self-organization 
phenomenon, assuming a concave, non-linear shape of 
the congested flow-density branch (cf. Sec. IV. B. 2). 
Instead of forming wide jams, narrow jams may co- 
exist when their distance is larger than about 2.5 km 
(Kerner, 1998a; Treiber et al, 2000). 

FIG. 23. Example for the time-dependent average veloc- 
ities in the three lanes of a freeway at subsequent cross 
sections (after Kerner, 1998a). The figures show a transi- 
tion from free to synchronized flow at cross section D5, the 
emergence of narrow jams in the "pinch region" at the up- 
stream cross section D4, and the formation of a few stable 
wide jams (see further upstream cross section Dl). The 
wide moving jams arise by the merging or disappearance 
of narrow jams between cross sections D4 and Dl. 



F. Cars and trucks 

Distinguishing vehicles of different lengths ("cars" 
and "trucks"), one finds surprisingly strong variations 
of the truck fraction, see Fig. 24 (Treiber and Hel- 
bing, 1999a). This point may be quite relevant for the 
explanation of some observed phenomena, as quanti- 
ties characterizing the behavior of cars and trucks are 
considerably different. For example, this concerns the 
distribution of desired velocities as well as the distri- 
bution of time headways (see Fig. 25). 

FIG. 24. The truck fraction is varying considerably in 
the course of time and should, therefore, be taken into 
account in empirical analyses of traffic flow. (From Treiber 
and Helbing, 1999a.) 



FIG. 25. Separate time-headway distributions for cars 
and trucks on the right lane (after Tilch and Helbing, 
2000). The time headways in front of trucks (long vehicles) 
are considerably longer than those of cars, in particular in 
congested traffic. 



G. Some critical remarks 

The collection and evaluation of empirical data is 
a subject with often underestimated problems. To 
make reliable conclusions, in original investigations 
one should specify 

(i) the measurement site and conditions (including ap- 
plied control measures), 

(ii) the sampling interval, 
(hi) the aggregation method, 

(iv) the statistical properties (variances, frequency 
distributions, correlations, survival times of traffic 
states, etc.), 

(v) data transformations, 

(vi) smoothing procedures, 

and the respective dependencies on them. 

The measurement conditions include ramps and 
road sections with their respective in- and outflows, 
speed limits, gradients, and curves with the respec- 
tively related capacities, furthermore weather condi- 
tions, presence of incidents, and other irregularities. 
Moreover, one should study the dynamics of the long 
vehicles separately, which may also have a significant 
effect (see Sec. II.F). 

A particular attention has also to be paid to the 
fluctuations of the data, which requires a statistical in- 
vestigation. For example, it is not obvious whether the 
scattering of flow-density data in synchronized flow 
reflects a complex dynamics due to non-linear inter- 
actions or whether it is just because of random fluc- 
tuations in the system. In this connection, I remind of 
the power laws found in the high-frequency variations 
of macroscopic quantities (see Sec. II.E.l). 

The statistical variations of traffic flows imply that 
all measurements of macroscopic quantities should be 
complemented by error bars (see, e.g., Hall et at, 
1986). Due to the relatively small "particle" num- 
bers behind the determination of macroscopic quanti- 
ties, the error bars are actually quite large. Hence, 
many temporal variations are within one error bar 
and, therefore, not significant. As a consequence, the 
empirical determination of the dynamical properties 
of traffic flows is not a simple task. Fortunately, many 
hard-to-see effects are in agreement with predictions 
of plausible traffic models, in particular with deter- 
ministic ones (see Sees. IV. B and IV. C. 2). 

Nevertheless, I would like to call for more refined 
measurement techniques, which are required for more 
reliable data. These must take into account correla- 
tions between different quantities, as is pointed out by 
Banks (1995). Tilch and Hclbing (2000) have, there- 
fore, used the following measurement procedures: In 
order to have comparable sampling sizes, they have av- 
eraged over a fixed number AV of cars, as suggested 
by Helbing (1997e). Otherwise the statistical error at 
small traffic flows (i.e. at small and large densities) 
would be quite large. This is compensated for by a 
flexible measurement interval AT. It is favourable that 
AT becomes particularly small in the (medium) den- 



sity range of unstable traffic, so that the method yields 
a good representation of traffic dynamics. However, 
choosing small values of AiV does not make sense, 
since the temporal variation of the aggregate values 
will mainly reflect statistical variations, then. In order 
to have a time resolution of about 2 minutes on each 
lane, one should select AV = 50, while AV = 100 can 
be chosen when averaging over both lanes. Aggregate 
values over both lanes for AN = 50 are comparable 
with 1-minute averages, but show a smaller statisti- 
cal scattering at low densities (compare the results in 
Helbing 1997a, c with those in Helbing, 1997e). 

Based on the passing times t° a of successive vehicles 
a in the same lane, we are able to calculate the time 
headways At a . The (measurement) time interval 

a +AN 

AT= ( 13 ) 

a=a +l 

for the passing of A7V vehicles defines the (inverse of 
the) traffic flow Q via 1/Q = AT /AN = (At a ), which 
is attributed to the time t = (t° a ) = J2 a t° a / AN - 

FIG. 26. The covariance between headways d a and in- 
verse velocities l/v a shows significant deviations from zero 
in congested traffic, while it approximately vanishes in free 
flow. Even after traffic has recovered, there seem to remain 
weak correlations between headways and vehicle speeds for 
a considerable time. These are probably a reminiscence of 
congestion due to platoons which have not fully dissolved. 
(After Tilch, 2001.) 

Approximating the vehicle headways by d a = 
v a At a , one obtains 

g-<^-(£)-<«® + «"(*■•;=)■ 

where cov(d a ,l/v a ) is the covariance between the 
headways d a and the inverse velocities l/v a . We ex- 
pect that this covariance is negative and particularly 
relevant at large vehicle densities, which is confirmed 
by the empirical data (see Fig. 26). Defining the den- 
sity p by 

1/p = (d a ) (14) 

and the average velocity V via Eq. (10), we obtain the 
fluid-dynamic flow relation (46) by the conventional 
assumption cov(d a ,l/v a ) = 0. This, however, over- 
estimates the density systematically, since the covari- 
ance tends to be negative due to the speed-dependent 
safety distance of vehicles. In contrast, the common 
method of determining the density via Q/(v a ) sys- 
tematically underestimates the density (see Fig. 27). 
Consequently, errors in the measurement of the flow 
and the density due to a neglection of correlations 
partly account for the observed scattering of flow- 
density data in the congested regime. However, a con- 
siderable amount of scattering is also observed in flow- 
occupancy data, which avoid the above measurement 
problem. 



FIG. 27. 50- vehicle averages of (a) the density Q(l/v a ) 
and (b) the conventionally determined density Q/{v a ), 
both times as a function of the density p according to 
the proposed definition (14). The usually assumed rela- 
tions arc indicated by solid lines. (From Tilch and Helbing, 
2000.) 



III. MODELLING APPROACHES FOR VEHICLE 
TRAFFIC 

In zcroth order approximation, I would tend to say 
that each decade was dominated by a certain mod- 
elling approach: 

(i) In the 50ies, the propagation of shock waves, i.e. 
of density jumps in traffic, was described by a fluid- 
dynamic model for kinematic waves. 

(ii) The activities in the 60ies concentrated on "mi- 
croscopic" car-following models, which are often called 
follow-the-leader models. 

(iii) During the 70ies, gas-kinetic, so-called Boltz- 
mann-like models for the spatio-temporal change of 
the velocity distribution were florishing. 

(iv) The simulation of "macroscopic", fluid-dynamic 
models was common in the 80ies. 

(v) The 90ies were dominated by discrete cellular au- 
tomata models of vehicle traffic and by systematic in- 
vestigations of the dynamical solutions of the models 
developed in the previous 50 years. 

(vi) Recently the availability of better data brings up 
more and more experimental studies and their com- 
parison with traffic models. 

Altogether, researchers from engineering, mathe- 
matics, operations research, and physics have prob- 
ably suggested more than 100 different traffic models, 
which unfortunately cannot all be covered by this re- 
view. For further reading, I therefore recommend the 
books and proceedings by Buckley (1974), Whitham 
(1974), Gerlough and Huber (1975), Gibson (1981), 
May (1981), Vumbaco (1981), Hurdle et a/. (1983), 
Volmullcr and Hamerslag (1984), Gartner and Wil- 
son (1987), Leutzbach (1988), Brannolte (1991), Da- 
ganzo (1993), Pave (1993), Snorek et al. (1995), Wolf 
et al. (1996), Lesort (1996), the Transportation Re- 
search Board (1996), Gartner et al. (1997), Helbing 
(1997a), Daganzo (1997a), Rysgaard (1998), Schreck- 
enberg and Wolf (1998), Bovy (1998), Brilon et al. 
(1999), Ceder (1999), Hall (1999), and Helbing et al. 
(2000c, 2001a, b). Readers interested in an exhaustive 
discussion of cellular automata should consult the de- 
tailed review by Chowdhury et al. (2000b). 

In the following subsections, I can only introduce a 
few representatives for each modelling approach (se- 
lected according to didactical reasons) and try to show 
up the relations among them regarding their instabil- 
ity and other properties (see Sees. IV. A to IV. D). 

Before, I would like to formulate some criteria for 
good traffic models: For reasons of robustness and cal- 
ibration, such models should only contain a few vari- 
ables and parameters which have an intuitive mean- 



ing. Moreover, these should be easy to measure, and 
the corresponding values should be realistic. In ad- 
dition, it is not satisfactory to selectively reproduce 
subsets of phenomena by different models. Instead, 
a good traffic model should at least qualitatively re- 
produce all known features of traffic flows, including 
the localized and extended forms of congestion. Fur- 
thermore, the observed hysteresis effects, complex dy- 
namics, and the existence of the various self-organized 
constants like the propagation velocity of stop-and-go 
waves or the outflow from traffic jams should all be 
reproduced. A good model should also make new pre- 
dictions allowing us to verify or falsify it. Apart from 
that, its dynamics should not lead to vehicle collisions 
or exceed the maximum vehicle density. Finally, the 
model should allow for a fast numerical simulation. 



A. Microscopic follow-the-leader models 

The first microscopic traffic models were proposed 
by Reuschel (1950a, b) and the physicist Pipes (1953). 
Microscopic traffic models assume that the accelera- 
tion of a driver- vehicle unit a is given by the neighbor- 
ing vehicles. The dominant influence on driving behav- 
ior comes from the next vehicle (a — 1) ahead, called 
the leading vehicle. Therefore, we obtain the following 
model of driver behavior from Eq. (3): 



dVgjt) 

dt 



«(*) 



+ /«,(«-!)(*)■ (15) 



Herein, f a ,( a ~i)(t) < describes the repulsive effect 
of vehicle (a — 1), which is generally a function of 

(i) the relative velocity Av a (t) = [v a (t) — v a -i(t)], 

(ii) the own velocity v a {t) due to the velocity- 
dependent safe distance kept to the vehicle in front, 

(iii) of the headway (brutto distance) d a {t) = 
[x a -i(t) — x a (t)] or the clearance (nctto distance) 
s a (t) = [d a {t) — l a -i], with l a meaning the length 
of vehicle a. 

Consequently, for identically behaving vehicles with 



v o, T a — t, and fa,(a-i) = f i we would have 



f a ,(a-i)(t) = f(s a (t),v a (t),Av a {t)) . 



(16) 



If we neglect fluctuations for the time being and in- 
troduce the traffic-dependent velocity 

v e (s a ,v a ,Av a ) = v + rf(s a , v a ,Av a ) , (17) 

to which driver a tries to adapt, we can considerably 
simplify the "generalized force model" (15): 



dv a 
~dT 



v e (s a ,v a ,Av a ) - v a 



(18) 



1. Non-integer car-following model 

Models of the type (18) are called follow-the-leader 
models (Reuschel 1950a, b; Pipes, 1953; Chandler 



et al, 1958; Chow, 1958; Kometani and Sasaki, 1958, 
1959, 1961; Herman et al, 1959; Gazis et al., 1959; 
Gazis et al, 1961; Newell, 1961; Herman and Gardels, 
1963; Herman and Rothery, 1963; May and Keller, 
1967; Fox and Lchmann, 1967; Hocfs, 1972). One of 
the simplest representatives results from the assump- 
tion that the netto distance is given by the velocity- 
dependent safe distance s*(v a ) = s' + Tv a , where T 
has the meaning of the (effective) safe time clearance. 
This implies s a (t) — s*(v a (t)) or, after differentiation 
with respect to time, dv a (t)/dt = [ds a (t)/dt]/T = 
[dd a (t)/dt]/T = [t> a _i(i) - v a (t)]/T. Unfortunately, 
this model does not explain the empirically observed 
density waves (see Sec. II.E.l). Therefore, one has to 
introduce an additional time delay At sa 1.3 s in adap- 
tation, reflecting the finite reaction time of drivers. 
This yields the following stimulus-response model: 



dv a (t + At) 



i(t) -«„(*)] . 



(19) 



Stimulus 



Herein, 1/T is the sensitivity to the stimulus. This 
equation belongs to the class of delay differential equa- 
tions, which normally have an unstable solution for 
sufficiently large delay times At. For the above case, 
Chandler et al. (1958) could show that, under the in- 
stability condition At/T > 1/2, a variation of indi- 
vidual vehicle velocities will be amplified. The exper- 
imental value is At/T w 0.55. As a consequence, the 
non-linear vehicle dynamics finally gives rise to stop- 
and-go waves, but also to accidents. In order to cure 
this, to explain the empirically observed fundamental 
diagrams, and to unify many other model variants, 
Gazis et al. (1961) have introduced a generalized sen- 
sitivity factor with two parameters mi and m 2 : 



1 

T 



1 



[v a (t + At)Y 



T [X a -l(t) -X a (t)] m i ' 



(20) 



The corresponding non-integer car-following model 
can be rewritten in the form [dv a (t + At)/dt}/[v a (t + 
At)] mi = (l/T Q )[dd a (t)/dt]/[d a (t)] m2 , which is solved 
by f mi (v a (t + At)) = c + cif m2 (d a (t)) with/ fe (z) = 
z 1 ^ k if k 7^ 1 and lnz otherwise (en, c\ being in- 
tegration constants). For mi ^ 1 and m2 / 1 we 
will discuss the stationary case related to identical ve- 
locities and distances. The vehicle density p is then 
given by the inverse brutto distance l/d a , and the 
corresponding equilibrium velocity V e is identical with 
v a . Therefore, we obtain the velocity-density relation 
V e (p) = V [l - (p/pm^)" 12 - 1 ] 1 ^ 1 -" 1 ^ with the free 
velocity Vo and the maximum density p max - Most of 
the velocity-density relations that were under discus- 
sion at this time are special cases of the latter formula 
for different values of the model parameters m\ and 
mi- Realistic fundamental diagrams result, for exam- 
ple, for the non-integer values m\ w 0.8 and m 2 ~ 2.8 
(May and Keller, 1967; Hocfs, 1972) or mi = 0.953 
and m2 — 3.05 (Kiihne and Rodigcr, 1991; Kiihne 
and Kroen, 1992). 



Finally, note that the above non-integer car- 
following model is used in the traffic simulation pack- 
age MITSIM and has recently been made more re- 
alistic with acceleration-dependent parameters (Yang 
and Koutsopoulos, 1996; Yang, 1997). However, other 
car-following models deserve to be mentioned as well 
(Gipps, 1981; Benekohal and Treiterer, 1988; del 
Castillo, 1996b; Mason and Woods, 1997). For a his- 
torical review see Brackstone and McDonald (2000). 



2. Newell and optimal velocity model 

One of the deficiencies of the non-integer car- 
following models is that it cannot describe the driving 
behavior of a single vehicle. Without a leading vehicle, 
i.e. for d a — > co, vehicle a would not accelerate at all. 
Instead, it should approach its desired velocity v a in 
free traffic. Therefore, other car-following models do 
not assume an adaptation to the velocity of the lead- 
ing vehicle, but an adaptation to a distance-dependent 
velocity v' e (d a ) which should reflect the safety require- 
ments and is sometimes called the "optimal velocity". 
While Newell (1961) assumes a delayed adaptation of 
the form 



v a (t + At) = v' e (d a (t)) = v e (s a (t)) , 



(21) 



Bando et al. (1994, 1995a, b) suggest to use the rela- 
tion 



v' e (d) = (vo/2) [tanh(c? — d c ) +t&nhd c 



(22) 



with constants v and d c , together with the optimal 
velocity model 



dv a (t) = v' e (d a (t)) - v a (t) 
dt t 



(23) 



which may be considered as a first-order Taylor ap- 
proximation v a (t + At) w [v a (t) + At dv a (t)/dt] of 
Eq. (21) with r = At. For the optimal velocity model 
one can show that small perturbations are eventually 
amplified to traffic jams, if the instability condition 



dv' e (d a ) dv e (s a ) 



dd n 



ds r 



> 



2t 



(24) 



is satisfied (Bando et al., 1995a), i.e., if we have large 
relaxation times r or big changes of the velocity v e (s a ) 
with the clearance s a (see Fig. 28). Similar investiga- 
tions have been carried out for analogous models with 
an additional explicit delay At (Bando et al, 1998; 
Wang et al, 1998b). 

FIG. 28. Trajectories according to the optimal velocity 
model by Bando et al. (1995a) for each fifth vehicle, if the 
overall density g on the circular road falls into the regime 
of linearly unstable traffic (from Helbing, 1997a). As in 
Fig. 19, small perturbations are amplified and lead to the 
formation of "phantom traffic jams". 



3. Intelligent driver model (IDM) 

As the optimal velocity model does not contain a 
driver response to the relative velocity Av a with re- 
spect to the leading vehicle, it is very sensitive to 
the concrete choice of the function v' e (d a ) and pro- 
duces accidents, when fast cars approach standing 
ones (Helbing and Tilch, 1998). To avoid this, one 
has to assume particular velocity-distance relations 
and choose a very small value of r, which gives un- 
rcalistically large accelerations (Bleile, 1997b, 1999). 
In reality, however, the acceleration times are about 
five to ten times larger than the braking times. More- 
over, drivers keep a larger safe distance and deceler- 
ate earlier, when the relative velocity Av a (t) is high. 
These aspects have, for example, been taken into ac- 
count in models by Gipps (1981), Kraufi et al. (1996, 
1997), Helbing (1997a), Bleile (1997b, 1999), Helbing 
and Tilch (1998), Wolf (1999), or Tomer et al. (2000). 

For illustrative reasons, I will introduce the so- 
called intelligent driver model (Treiber and Helbing, 
1999b; Treiber et al, 2000). It is easy to calibrate, 
robust, accident-free, and numerically efficient, yields 
realistic acceleration and braking behavior and repro- 
duces the empirically observed phenomena. Moreover, 
for a certain specification of the model parameters, 
its fundamental diagram is related to that of the gas- 
kinctic-based, non-local traffic (GKT) model, which is 
relevant for the micro-macro link we have in mind (see 
Sec. III.E). For other specifications, the IDM is related 
with the Newell model or the Nagel-Schrcckcnberg 
model (see Sec. III.B.l). 

The acceleration assumed in the IDM is a contin- 
uous function of the velocity v a , the clearance s a = 
(d a — l a -i), and the velocity difference (approaching 
rate) Av a of vehicle a to the leading vehicle: 



dv a 
~~dt 



= a D 



s*Jv a ,Av a ) 



■ (25) 



This expression is a superposition of the acceleration 
tendency a a [l — (v a /v®) s ] on a free road, and the de- 
celeration tendency f a i a -\) = — a a [s* a (v a , Av a )/s a ] 2 
describing the interactions with other vehicles. The 
parameter 6 allows us to fit the acceleration behavior. 
While 6 = 1 corresponds to an exponential-in-time 
acceleration on a free road, as assumed by most other 
models, in the limit 6 — > oo, we can describe a con- 
stant acceleration with a a , until the desired velocity 
is reached. The deceleration term depends on the 
ratio between the "desired clearance" s* and the ac- 
tual clearance s a , where the desired clearance 



v a Av a 




s* a {v a ,Av a ) 



is dynamically varying with the velocity v a and the 
approaching rate Av a , reflecting an intelligent driver 
behavior. The IDM parameters can be chosen individ- 
ually for each vehicle a, but for the moment, we will 



assume identical vehicle parameters and drop the in- 
dex a for readability. These are the desired velocity vo, 
the safe time clearance T, the maximum acceleration 
a, the comfortable deceleration b, the acceleration ex- 
ponent 6, the jam distances s' and s", and the vehicle 
length I, which has no dynamical influence. To reduce 
the number of parameters, one can assume 6 = 1, 
s" = 0, and I = 0, which still yields good results. 

In equilibrium traffic with dv a /dt = and Av a = 0, 
drivers tend to keep a velocity-dependent equilib- 
rium clearance s e {v a ) to the front vehicle given by 
s e (v) = s*(v,0)[l - (v/vo) 5 }- 1 / 2 . Solving this for the 
equilibrium velocity v = v e leads to simple expres- 
sions only for s" = and 6 = 1, 6 = 2, or 6 — > oo. 
In particular, the equilibium velocity for the special 
case 6 = 1 and s' = s" = is u e (s) = { — 1 + 
v/1 + [4T 2 (v ) 2 ]/s 2 }s 2 /(2v T 2 ). From this equation 
and the micro-macro relation s = (d — I) = (1/p — 
I) = (1/p — 1/pmax) between clearance and den- 
sity follows the corresponding equilibrium traffic flow 
Qe(p) = pVeip) as a function of the traffic density p. 
The acceleration coefficient 6 influences the transition 
region between the free and congested regimes. For 



and , 



0, the fundamental diagram becomes 



triangular-shaped: Q e (p) = min(pw , [1 — p(l + s')]/T). 
For decreasing 6, it becomes smoother and smoother. 



B. Cellular automata (CA) 

Cellular automata are interesting for their speed 
and their complex dynamic behavior (Wolfram, 1984, 

1986, 1994; Stauffer, 1991), including such fascinat- 
ing phenomena as self-organized criticality (Bak et al., 

1987, 1988; Bantay and Janosi, 1992; Olami et al., 
1992), formation of spiral patterns (Markus and Hess, 
1990), or oscillatory and chaotic sequences of states 
(Wolfram, 1984; Markus and Hess, 1990; Nowak and 
May, 1992). Their enormous computation speed and 
efficiency is a consequence of the following proper- 
ties, which are ideal preconditions for parallel com- 
puting: (i) discretization of space into identical cells 
(sites) j of size Ax, (ii) a finite number of possible 
states g(x) (iii) the (parallel) update at times t = i At 
with an elementary time step At, and (iv) globally 
applied update rules, based on (v) short-range inter- 
actions with a finite (small) number of neighbouring 
sites. Despite these simplifications, cellular automata 
and related lattice gas automata have a broad range 
of applications, reaching from realistic simulations of 
granular media (Peng and Herrmann, 1994) or fluids 
(Frisch et al., 1986; Chen et al., 1991) (including inter- 
facial phenomena and magnetohydrodynamics), over 
the computation of chemical reactions (Markus and 
Hess, 1990; Dab et al, 1991), up to the modeling of 
avalanches (Bantay and Janosi, 1992). 

Their application to traffic dynamics has stim- 
ulated an enormous research activity, aiming at 
an understanding and control of traffic instabili- 



ties, which are responsible for stop-and-go traffic 
and congestion, both on freeways and in cities. 
The first cellular automata for freeway traffic go 
back to Cremer and coworkers (Cremer and Lud- 
wig, f986; Schiitt, 1990), and Nagel and Schreck- 
enberg (1992). Since then, there have been over- 
whelmingly many proposals and publications in this 
field. Nevertheless, I will keep my discussion short, 
since it has recently been reviewed by Chowdhury 
et al. (2000b). Complementary, it is worth check- 
ing out the Java applets supplied at the web pages 
http : // r cswww . urz . tu-dresden . de/ ~helbing/Road 
Applet/ and http : //www. traffic.uni-duisbu.rg. 
de/model/index.html, which allow the interested 
reader to compare the dynamics of various cellular 
automata. 

1. The Nagel-Schreckenberg model and its slow-to-start 
variant 

Cellular automata describe the vehicle dynamics 
less detailed than follow-the- leader models, but their 
simplifications are in favour of an extremely fast sim- 
ulation of huge numbers of interacting vehicles. Nagel 
and Schreckenberg (1992) suggest to divide the street 
into cells j of length Ax and the time t into intervals 
i of duration At = 1 s. Each cell is either empty or 
occupied by one vehicle with speed 

Ax 

where Vi £ {0, 1, . . . , w max }. For freeways, one fre- 
quently chooses the cell size Ax = 7.5 m and the 
scaled desired velocity v max = v^At/Ax = 5. The 
car positions are updated in parallel according to the 
following rules: 

(i) Motion: Move the vehicle forward by ii cells. 

(ii) Acceleration: If a vehicle has not yet reached its 
maximum velocity, it would like to drive with the in- 
creased velocity v\ = (ii + 1), corresponding to a con- 
stant acceleration Ax /(At) 2 . 

(iii) Deceleration: If the distance (i.e., the number of 
cells) to the next vehicle ahead is di < v[, the velocity 
is reduced to v" = (di — 1), otherwise v" — v[. 

(iv) Randomization: With probability p, the velocity 
is reduced to v i+ i — (v'J — 1), if this yields a non- 
negative velocity, otherwise Vi + \ = 0. 

According to steps (ii) to (iv), the updated veloc- 
ity as a function of the previous velocity ii can be 
summarized by the formula 

v i+ i = max (0, min(u max , di - 1, t), + 1) - £- p) ) , (28) 

where the Boolean random variable Q = 1 with 
probability p and otherwise. 

Note that, according to rule (iii), drivers move al- 
ways below the speed which advances them by the 



clearance to the car ahead within At. Hence, the pa- 
rameter At plays, at the same time, the role of an up- 
dating step, of the adaptation time r, and of the safe 
time clearance T. Moreover, the discretization length 
Ax agrees with the effective vehicle length /', i.e., 
with the inverse jam density. This makes the Nagel- 
Schreckenberg model so extremely compact and ele- 
gant. 

FIG. 29. Trajectories of each tenth vehicle for unstable 
traffic in the Nagel-Schreckenberg model, if the slowdown 
probability is (a) p — 0.5 and (b) p = 0.001. In the limit 
p — > 0, the jam amplitude goes to zero. 

The model parameter ("slowdown probability") p 
describes individual velocity fluctuations due to de- 
layed acceleration (imperfect driving). For freeway 
traffic, Nagel and Schreckenberg have often set p = 
0.5, which leads to a relatively noisy dynamics (see 
Fig. 29), while p = 0.2 and u max — 2 are suitable val- 
ues for city traffic (Esser and Schreckenberg, 1997). 

Another model variant with a slow-to-start rule, the 
velocity-dependent randomization model (Barlovic 
et al., 1998), assumes p = 0.01 for finite velocities, 
while po = 0.5 for = 0. A similar model was pro- 
posed by Benjamin et al. (1996). The fundamental di- 
agram of the resulting model combines the proper- 
ties of the Nagel-Schreckenberg model with p = 0.01 
at low densities with the properties of the one with 
p = 0.5 at high densities. In between, one finds 
some "crossover region" , in which the solutions of 
both models can exist, but the high- flow states are 
metastable. We will come back to this and other top- 
ics later (see Sec. IV. A). 

2. Some other cellular automaton models 

An interesting variant of the Nagel-Schreckenberg 
model results in the cruise control limit (Nagel and 
Paczuski, 1995), where the fluctuations are turned 
off for di = v max . Then, the life-time distribution 
of traffic jams changes from an exponential (Nagel, 
1994) to a power-law distribution, which implies frac- 
tal self-similarity and points to self-organized critical- 
ity (Nagel, 1994; Nagel and Rasmussen, 1994). 

The cellular automaton by Brilon and Wu (1999) 
modifies the Nagel-Schreckenberg model by addition- 
ally introducing an acceleration probability and re- 
quiring a minimal time clearance for acceleration. 

Takayasu and Takayasu (1993) were the first to in- 
troduce a slow-to-start rule. According to a general- 
ization of their TT or T2 model (Schadschneider and 
Schreckenberg, 1997a), a standing vehicle with veloc- 
ity Vi = will accelerate with probability q = (1 — p), 
if there is exactly one empty cell in front, for di > 2 
it will deterministically accelerate to iii+i = 1. 

The model by Nagel and Herrmann (1993) can 
be viewed as a continuum version of the Nagel- 



Schrcckenberg model. Its slightly generalized version 
by Sauermann and Herrmann (1998) reads 

{max(c?j — 5, 0) for v t > di — a, 
min(wi + a, w mas ) for < rfj - 0, 
Vi otherwise. 

Herein, the acceleration coefficient is determined by 
d. = fimax min(l, di/j), where a, (3, 7, and 5 are the 
parameters of the model. 

Another collision- free continuum-in-space version of 
the Nagel-Schreckenberg model is related to the model 
by Gipps (1981) and given by the equation 

v a (t + l) = max(0,rnd« s (t)-aAi,< s (t))) , 

where rnd(zi, Z2) represents a random number uni- 
formly distributed in the interval [21,22)- The 
desired velocity is determined by v^ cs (t) = 
min (u™ ax , v^ ie (t), v + a At), while the safe veloc- 
ity is calculated via v^ fc (t) = v a -i(t) + b[s a — 
v a -i{t) At]/[v a (t) + b At]. This model was developed 
and investigated in detail by KrauB and coworkers 
(Kraufi et at, 1996, 1997; KrauB, 1998a, b). 

The model developed by Ianigro (1994) is based 
on Petri nets. Yukawa and Kikuchi (1995, 1996) have 
studied traffic models based on coupled maps (see also 
Tadaki et al., 1998). Although there are many more 
cellular automata models of traffic flow (e.g., Wolf, 
1999), I close with pointing to the model v i+1 = 
max (0, min ( Umax, — 1) — with unlimited ac- 

celeration capabilities by Fukui and Ishibashi (1996a; 
see also Braudc, 1996; Wang et at, 1998a). 

Most cellular automata models have a principal dif- 
ference to other models of traffic flow. In contrast 
to deterministic approaches such as the above car- 
following models, they require fluctuations for the ex- 
planation of traffic jams. That is, in the limit p — > 0, 
traffic jams will normally disappear (see Fig. 29). This 
does not disqualify cellular automata models, since 
driver behavior is certainly imperfect and traffic flow 
subject to random influences. However, this raises the 
following questions: (i) Are fluctuations a dominant 
or a subordinate effect? (ii) Which kinds of observa- 
tions require the consideration of noise? (iii) Can cellu- 
lar automata models of traffic flow be mathematically 
connected with other traffic models? 

3. The discrete optimal velocity model 

To construct a mathematical link between the 
Nagel-Schreckenberg model and the optimal veloc- 
ity model, let us discretize the latter: v a (t + At) = 
v a (t) + [v' e (d a (t)) — v a (t)]At/r. Dropping the vehicle 
index a and scaling the time by At = 1 s, distances by 
the cell length Ax (to be specified later on), and ve- 
locities by Ax/ At results in t)j + i = + X[v' e (di) — Vj\, 
where A = At/r. In order to have integer-valued veloc- 
ities v and locations, we need to introduce a tabular 



function h(z). Additionally, we will add some noise 
for comparison with the Nagel-Schreckenberg model. 
The finally resulting discrete and noisy optimal veloc- 
ity model reads (Helbing and Schreckenberg, 1999): 

= max (0,Vi + h{\[v' e {d t ) - v,}) - . (29) 

For the floor function h(z) — [z\ , the argument z is 
rounded down to the largest natural number n < z. 
Then, the above equation implies t)j+i < Xv' e (di) + 
(1 - \)Vi- 

FIG. 30. Trajectories of each tenth vehicle in the con- 
gested regime according to the discrete version of the op- 
timal velocity model for slowdown probability (a) p — 0.5 
and (b) p = 0.001. Jam formation does not disappear for 
small values of p. These just reduce the average distance 
between traffic jams, while their amplitude remains con- 
stant (see Helbing and Schreckenberg, 1999). 

4. Comparison 

The above model is practically as efficient and fast 
as the Nagel-Schreckenberg model, but it allows a fine- 
grained description of vehicle velocities and locations 
by selecting small values of Ax. Moreover, it can be 
specified in a way that vehicles will normally not jump 
into the same cell or over each other. However, like 
in the continuous optimal velocity model itself, the 
discrete version is sensitive to the choice of the pa- 
rameter A = At/r and the integer- valued velocity- 
distance function v' e (d). A value A « 0.77 seems to be 
optimal. For city traffic, the additional specifications 
v' e (d) = min(d — 1,3) and Ax — 6.25 m yield realis- 
tic results. The favourite properties of this model are 
that (i) jam formation is not affected by reducing the 
fluctuation strength (see Fig. 30), (ii) the characteris- 
tic constants of traffic are reproduced (see Sees. II.E.l 
and IV. A. 5), (iii) the critical densities and the char- 
acteristic constants can be calculated, (iv) a relation 
with other cellular automata can be constructed. 

For example, the Nagel-Schreckenberg model corre- 
sponds to the selection A = 1 and v' e (d) = mm(d — 
l,w m ax). Moreover, h(z) = [z\ , if z < 1, otherwise 
h(z) — 1, which limits acceleration, while decelera- 
tion capabilities are unlimited. Nevertheless, these two 
mathematically related models lead to a different dy- 
namics (compare Figs. 29 and 30), probably because 
the assumed velocity changes dv' e (d) / dd with the dis- 
tance assumed in the Nagel-Schreckenberg model are 
not sufficiently large for a linear instability, and the 
ratio between acceleration and braking capabilities is 
too small (see Sec. IV.A.4, Fig. 37). 

C. Master equation 

Some particle hopping models (see Sec. I.E) simplify 
the description of transport processes even further 



than the above cellular automata do. This is for two 
reasons: (i) to identify the minimal requirements for 
certain phenomena in self-driven many-particle sys- 
tems and (ii) to facilitate analytical investigations, al- 
though their treatment is not necessarily easier. How- 
ever, the main difference is that cellular automata are 
discrete in time rather than continuous, allowing us to 
specify deterministic processes with transition proba- 
bility 1. 

For a stochastic description, let us consider a system 
with L states j occupied by rij particles. The distri- 
bution P(n,t) shall denote the probability of finding 
the system configuration, i.e. the occupation vector 



n = (m, . . . ,nj,n j+ i, . . . ,n L ) 



(30) 



at time t. This probability is reduced by transitions 
to other configurations n' , the frequency of which is 
proportional to P(n,t). The proportionality factor is 
the conditional transition probability P(n' ,t + At\n, t) 
of finding the configuration n' at time (t + At), given 
that we have the configuration n at time t. Conversely, 
the probability P{n,t) increases by transitions from 
configurations n 1 to n, which are proportional to the 
occurence probability P(n',t) of n' and to the tran- 
sition probability P(n,t + At\n',t). Considering the 
normalization J2 n > -P( n '> t+At\n, t) — 1, the resulting 
balance equation govering the dynamics of so-called 
Markov chains reads 

P(n, t + At) = J2 p ( n > 1 + A< l n '> ^(V > *) ( 31 ) 

n' 

or, in the continuous limit At — > 0, 
dP(n,t) 



dt 



= W(n\n';t)P(n',t) 

n'(^n) 

W(n'\n;t)P(n,t) , 



E 

n'(^n) 



(32) 



where we have introduced the transition rates 



... ,. P(n,t + At \n',t) , . , , 

W(n\n';t)= lim — — 1 — (n'^n). 

v 1 ; At^O At 

Note that the master equation (32) assumes the con- 
ditional probabilities P(n,t + At\n',t) to depend on 
t and At only, but not on previous time steps. 
For the TASEP (see Sec. I.E.I) we have 



Pin 1 t + At\n t) = l ^(l-%+i) if n' 



otherwise, 



(33) 



where 



rijji = (ni, . . . ,rij - 1, . . . ,rij> + 1, . . . ,n L ) ■ 
In the continuous limit, we assume limAt->o <z/Ai = v. 



1. Solution methods, mapping to spin chains, and matrix 
product ansatz 

Since the master equation is linear, many solution 
methods are available for it (cf. the books and reviews 
by Feller, 1967; van Kampen, 1981; Haken, 1977; Gar- 
diner, 1985; Haus and Kchr, 1987; Weiss, 1994; Hel- 
bing, 1995a; Schiitz, 2000; see also Hclbing and Molini, 
1995). For example, one can map the TASEP to a spin 
system by identifying empty sites with down spins and 
occupied ones with up spins. Following Doi (1976; see 
also Sandow and Trimper 1993b; Kaulke and Trimper, 
1995), the probability distribution P(n,t) can be re- 
lated to a state vector \P{t)) in Fock space, and with 
respect to its orthonormal basis {|«fc)} we have the 
relation \P(t)) — J2k P( n k, t)\nk). The linear master 
equation becomes d\P(t))/dt — C\P(t)), where the 
linear operator L can, similar to a second quantization, 
be expressed in terms of creation and annihilation 
operators usually satisfying Bose commutation rules 
(Doi, 1976; Grassberger and Scheunert, 1980; Peliti, 
1985; Sandow and Trimper, 1993a). When the occupa- 
tion numbers per lattice site are limited to a maximum 
value, one has to take into account the exclusion prin- 
ciple, leading to commutation rules of Pauli operators 
or similar ones (Gwa and Spohn, 1992; Sandow and 
Trimper, 1993b; Rudavets, 1993; Schiitz and Sandow, 
1994; Alcaraz et al, 1994). 

Alternatively, one may rewrite the master equation 
as Schrodinger equation in imaginary time by set- 
ting H = —C (Felderhof and Suzuki, 1971; Siggia, 
1977; Alexander and Holstcin, 1978; Schiitz, 2000). 
The stochastic Hamiltonian H is defined by its ma- 
trix elements (nk\Tt\nk) = J2i(^k) W(ni\rik), and 
(nk\H\ni) = — W(nk\ni). The stationary state corre- 
sponds to the eigenvector of Tt or C with eigenvalue 0, 
i.e. it can be interpreted as ground state. For discrete 
time dynamics, the corresponding equation (31) can 
be analogously rewritten as \P(t + At)) = T\P{t)). 
Here, the stationary state is the eigenvector of the 
so-called transfer matrix T with eigenvalue 1. It can 
be determined by a so-called matrix product ansatz 
(Derrida et al, 1993; Hinrichsen, 1996; Rajewsky and 
Schreckenberg, 1997; Derrida and Evans, 1997; Ra- 
jewsky et al, 1998; Klauck and Schadschneider, 1999; 
Karimipour, 1999c; see also Krebs and Sandow, 1997, 
for an existence theorem). The stationary solution is 



Pst(n) = ^-(W 



H[njV+(l-nj)£] 

3 = 1 



V), (34) 



where \V) and \W) are vectors characterizing the 
boundary conditions, T> and £ are suitable matrices, 
and Zl = (W\(T>+£) L \V) is a normalization constant. 
More general formulas are available for dynamic solu- 
tions, using Bcthe ansatz equations (Stinchcombe and 
Schiitz, 1995a, b; Sasamoto and Wadati, 1997; Schiitz, 
1998). 



2. Mean field approach and Boltzmann equation 



3. TASEP and Nagel-Schreckenberg model 



It is often useful to consider the mean value equa- 
tions for the expected values (rij) = ^2 n rijP(n,t), 
which we obtain by multiplying Eq. (32) with rij, sum- 
ming up over n, and suitably interchanging n and n': 

d ( n j) dP(n,t) 



For the TASEP with the transition rates (33), the 
mean field approach (40) yields 



d ( n j) 
dt 



= Qi-i(t) - Qj{t) 



(41) 



dt ^ J dt 

n 

Here, we have introduced the first jump moments 

m,j(n,t) = J^K' - n 3 )W(n'\n;t) . (36) 

n' 

Let us assume spontaneous transitions with individual 
transition rates w(j'\j) from state j to f and pair 
interactions changing the states from j and k to j' 
and k' with transition rates W2(j', k'\j, k). Defining 



{rrij (n, t)) . (35) witn t ^ e pa ,rticle flow 



n 



kk' . 
jj' 



■ (. . .rij -I . . .rij' + l . . .rife-1 . . . rifc' + l ■ ■ •) , 

the corresponding configurational transition rates are 
given by 



w(j'\j)nj 



if n' 



n 



jj' > 



W(n'\n) = { w 2 (j', k'\j, k)n 3 n k if n> = nf', , (37) 







jj 

otherwise. 



That is, the total rate of spontaneous transitions is 
proportional to the number n 3 of systems which may 
change their state j independently of each other, while 
the total rate of pair interactions is proportional to the 
number n 3 n k of possible interactions between systems 
in states j and k. Inserting Eq. (37) into (36) eventu- 
ally leads to 



mj{n) = w(j\f) + E w 2{j, k\f, k')n k , 



k,k' 



- E w (j'\j) + E w 2(j', k'\j, k)n k 



k.k' 



(38) 



(see, for example, Helbing, 1992a, 1995a). The mean 
field approach assumes 



{m 3 {n)) w m 3 {(n}), 



(39) 



i.e., that the system dynamics is determined by the 
mean value (n), which is true for a sharply peaked, 
unimodal distribution P(n,t). This leads to the gen- 
eralized Boltzmann equation 



dp{j, t) 
dt 



j' k,k' 



E 

j' 



w(j'\j)+^w 2 (f,k'\j,k)p(k,t) 

k,k 



P(f,t) 
p(.i, t) , 

(40) 

where we have introduced the densities p(j,t) = 
(rij) /Ax and W2(j',k'\j,k) — W2W, k'\j, k) Ax. Note 
that this Boltzmann equation neglects the covari- 
ances o- jk (t) = {(n 3 - (n 3 ))(n k - (rife))) = ((n 3 n k ) - 
{rij}{n k }) and the corresponding correlations r 3k (t) = 

Vjk/y/O-jjO-kk- 



Qj(t) = v(rij(l - n 3+1 )) w v(n 3 )(l - (rij)) . (42) 

While in open systems exact calculations are difficult 
clue to correlations, this relation becomes exact for 
the stationary state of periodic systems (Spohn, 1991; 
Schadschneider and Schreckenberg, 1993; Schrecken- 
berg et al, 1995; Liggett, 1999). A spatially continu- 
ous approximation of Eq. (41) results in the Burgers 
equation (53) or, in lowest-order approximation, in the 
Lighthill-Whitham equation (49). 

However, replacing the random sequential update 
of the TASEP by a parallel update, we end up with 
the traffic model by Nagel and Schreckenberg with 
Vmax = 1 and p = (1 — q) (see Sec. III.B.l). For a 
parallel update, the state changes of many particles 
are coupled, leading to significant correlations. Conse- 
quently, we have to take into account higher-order in- 
teractions in the respective transition rates W(n'\n). 
It turns out that the above result for the TASEP corre- 
sponds to the site-oriented mean-field theory (SOMF) 
for the Nagel-Schreckenberg model with w max = 1, 
which underestimates the flows by neglecting the cor- 
relations (Nagel and Schreckenberg, 1992; Schrecken- 
berg et al, 1995). This can be corrected for by elim- 
inating so-called paradisical or garden-of-Eden states 
(Moore, 1962) that cannot be reached dynamically, 
when a parallel update is applied (Schadschneider and 
Schreckenberg, 1998). Alternatively, one can develop 
a car-oriented mean-field theory (COMF), which cal- 
culates the probabilities P n (t) of finding n empty sites 
in front of a vehicle (Schadschneider and Schreck- 
enberg, 1997b). A superior method (Schadschnei- 
der, 1999) is a site-oriented cluster-theoretic approach 
(Schadschneider and Schreckenberg, 1993; Schrecken- 
berg et al, 1995), calculating occurence probabilities 
of states composed of n successive sites (Kikuchi, 1966; 
Gutowitz et al., 1987; ben-Avraham and Kohler, 1992; 
Crisanti et at, 1993). Instead of Eq. (39), a 2-cluster 
approximation gives 

(nj-i{t)nj (t)n 3+1 (t)n 3+2 (i)) 
_ (n 3 -_i (t)n 3 (t)) (n 3 (t)n j+1 (t))(n 3+1 (t)n 3+2 (t)) 
(n 3 (t))(n 3+1 (t)) 

This allows us to calculate the time-dependent aver- 
age flow Q(t) = q((rij(t)(l - n 3+1 (t)))) / At = Q(t)/At 
(Wang and Hui, 1997), if we take into account that 
the occupation numbers change according to 

n 3 (t + 1) = (1 - ef_7 ) (*))n j _ 1 (t)(l - Uj(t)) 

+ nj(t)n j+1 (t) + &\t)n 3 (t)(l - n j+1 (t)) . 



In the above formulas, ((rij)) — gAx indicates an aver- 
age over all lattice sites j. For the stationary case, one 
can finally derive the correct stationary flow-density 
relation 



Q = 1{1 - [1 - 4g«n j »(l - «n,) 



11/2 



} (43) 



(Yaguchi, 1986; Schadschncidcr and Schrcckcnberg, 
1993, 1997b, 1998; Schreckenberg et al, 1995). The 
above results can be gained using (£,j P \t)) = p, 



0, l&\t)¥ = &\t), as well as 



the statistical independence of £j P \t), ^^(t), and 

^ p \t + 1) (Wang and Hui, 1997). The dynamics 
and correlation functions are also different from the 
TASEP (Schadschneider, 1999). 



4. Nucleation and jamming transition 

A master equation model for the jamming transi- 
tion on a circular road of length L has been devel- 
oped by Mahnke and Pieret (1997). It assumes that 
congested regions are eventually merging to form one 
big "megajam". If the spatial dynamics of the N 
cars is neglected, we need only one single occupa- 
tion number n to count the number of cars in the 
traffic jam. Moreover, the model assumes that cars 
leave a standing traffic jam with the constant rate 
w(n — l|n) = 1/T. They join the queue with the rate 
w(n+l\n) = v e (s(n)) / s(n) given by the vehicle speed 
v e (s) = t>ns 2 /[(sn) 2 + s 2 ], divided by the average clear- 
ance s(n) = (L — Nl)/(N — n) of freely moving vehi- 
cles (so being some constant). Mahnke and Kaupuzs 
(1999) derive the stationary distribution P s t(n) of the 
queue length n and show a transition from free to 
congested traffic at a certain critical vehicle density. 
While the maximum of the probability distribution 
P s t(n) is located at n = 1 in free traffic, in congested 
traffic it is located at some finite value n > 1, which 
is growing with the vehicle density. 

Mahnke and Kaupuzs compare the nucleation, 
growth, and condensation of car clusters with the 
formation of liquid droplets in a supersaturated va- 
por (Schweitzer et al, 1988; Ebcling et al, 1990). 
Moreover, they recognize that, for fixed scaled pa- 
rameters and fixed vehicle length I, the fundamental 
flow-density diagram of their model is dependent on 
the road length L. This, however, results from the 
fact that they relate the flow to the overall vehicle 
density g — N/L rather than to the local density 
p=l/[l + s{n)\. 



5. Fokker-Planck equation 

Kiihne and Anstett (1999) have obtained similar re- 
sults like Mahnke and coworkers, and they have man- 
aged to derive the distribution P st (n) and the location 



of its maximum as well. Their calculations are based 
on the related Fokker-Planck equation 

"' =-^[m(y,t)p(y,t)} + ^^[m 2 (y,t)P(y,t)} 

(44) 



dt 



for the scaled, quasi-continuous variable y = n/N. 
This Fokker-Planck equation can be viewed as 
a second-order Taylor approximation of the mas- 
ter equation, after this has been rewritten in the 

form dP(y,t)/dt = E y >[W(y', y - y')P(y - y', t) - 

W(y',y)P(y,t)} with W(y',y) = W((y + y')N\yN) = 
W(n + y'N\n). Correspondingly, the first jump mo- 
ment, which has the meaning of a drift coefficient, 
is rh{y,t) — ^ , y'W(y' ,y), while the second jump 
moment, which can be interpreted as diffusion coeffi- 
cient, is given by rai(y, t) = ^2 y , (y') 2 W(y', y). Details 
regarding the Fokker-Planck equation and its solution 
can be found in the books by van Kampcn (1981), 
Horsthemke and Lefever (1984), Gardiner (1985), and 
Risken (1989). 



D. Macroscopic traffic models 

As we will see in Sec. IV. A. 4, the above findings 
of Mahnke, Kiihne and coworkers regarding the fun- 
damental diagram and the nucleation effect are well 
compatible with those for macroscopic traffic mod- 
els. In contrast to microscopic traffic models, macro- 
scopic ones restrict to the description of the collective 
vehicle dynamics in terms of the spatial vehicle den- 
sity p(x,t) per lane and the average velocity V(x,t) 
as a function of the freeway location x and time t. 
Macroscopic traffic models have been often preferred 
to car-following models for numerical efficiency, but in 
terms of computation speed they cannot compete with 
cellular automata. However, other advantages are (i) 
their better agreement with empirical data, (ii) their 
suitability for analytical investigations, (iii) the sim- 
ple treatment of inflows from ramps (see Sec. IV. B), 
(iv) the possibility to simulate the traffic dynamics in 
several lanes by effective one-lane models considering 
a certain probability of overtaking (see Sec. III.E). 



1. The Lighthill-Whitham model (LW model) 

The oldest and still most popular macroscopic traf- 
fic model goes back to Lighthill and Whitham (1955). 
It appears that Richards (1956) developed the same 
model independently of them. Their fluid-dynamic 
model is based on the fact that, away from on- or off- 
ramps, no vehicles are entering or leaving the freeway 
(at least if we neglect accidents). This conservation of 
the vehicle number leads to the continuity equation 



dp(x, t) dQ(x, t) 



dt 



dx 



0. 



(45) 



Herein, 



Q{x,t) = p(x,t)V{x,t) 



(46) 



is the traffic flow per lane, which is the product of the 
density and the average velocity (see Sec. III.E). 

We may apply the so-called total or substantial 
derivative 

dt dt dx ' 
describing temporal changes in a coordinate system 
moving with velocity V(x, t), i.e. with "the substance" 
(the vehicles). With this, we can rewrite (45) in the 
form d v p(x,t)/dt = —p(x 1 t)dV(x 1 t)/dx, from which 
we conclude that the vehicle density increases in time 
{d v p/dt > 0), where the velocity decreases in the 
course of the road (dV/dx < 0), and vice versa. More- 
over, the density can never become negative, since 
p{x,t) = implies d v p(x,t)/dt = 0. 

Eq. (45) is naturally part of any macroscopic traf- 
fic model. The difficulty is to specify the traffic flow 
Q(x,t). Lighthill and Whitham assume that the flow 
is simply a function of the density: 

Q(x,t) = Q e (p(x,t)) = pV e (p(x,t)) > 0. (47) 

Herein, the so-called fundamental (flow-density) dia- 
gram Q e {p) and the equilibrium velocity-density re- 
lation V e (p) are thought to be suitable fit functions 
of empirical data, for which there are many proposals 
(see, for example, Sec. III. A). The first measurements 
by Greenshields (1935) suggest a linear relation of the 
form 

V e (p) = V (l-p/p iam ), (48) 

which is still sometimes used for analytical investiga- 
tions. 

Inserting Eq. (47) into the continuity equation (45), 
we obtain 



dp , nln , dp 



(49) 



This is a non-linear wave equation (Whitham, 1974, 
1979), which describes the propagation of so-called 
kinematic waves with the velocity 



n( \ d ® e \r ( \ , dVe 
Cip) = ^p- = V ^ )+P ^ p -- 



(50) 



Because of dV e (p)/dp < 0, we have C(p) < V e (p). 
Hence, the kinematic waves always propagate back- 
wards with respect to the average velocity V e (p) of 
motion, namely with the speed c(p) = [C(p) — V e (p)] = 
pdV e {p)/dp<0. 

FIG. 31. (a) Trajectories of each tenth vehicle and (b) 
corresponding spatio-temporal density plot illustrating the 
formation of a shock wave on a circular road according 
to the Lighthill- Whitham model (from Helbing, 1997a). 
Although the initial condition is a smooth sinusoidal wave, 
the upstream and downstream fronts are getting steeper 
and steeper, eventually producing discontinuous jumps in 
the density profile, which propagate with constant speeds. 
The wave amplitude remains approximately unchanged. 



Note that C(p) is the speed of the characteristic 
lines (i.e. of local information propagation), which is 
density dependent. In contrast to linear waves, the 
characteristic lines intersect, because their speed in 
congested areas is lower. This gives rise to changes 
of the wave profile, namely to the formation of shock 
fronts, while the amplitude of kinematic waves does 
not change significantly (see Fig. 31). The property of 
a shock front is characterized by the upstream value 
P- and the downstream value p + of the density (see 
Fig. 32). The difference in the flow Q_ upstream and 
the flow Q + downstream causes the movement of the 
density jump (p + — p_) with speed S. Because of ve- 
hicle conservation, we have (Q+ — Q_) = (p+ — p~)S. 
Consequently, with Eq. (47), the speed of shock prop- 
agation is 

S(p+, p_) = [Q e (p+) - Qe(p-)]/(p+ - P-) (51) 
(see Fig. 32). 

FIG. 32. Illustration of the geometrical relations between 
the fundamental diagram Q e (p) (solid line) with other 
quantities characterizing traffic flow. The vehicle velocity 
V e (p) corresponds to the slope of the line connecting the 
origin (0,0) and the point (p,Q e (p)) (■■■), the propaga- 
tion speed C(p) of density waves to the slope of the tan- 
gent dQ e (p)/dp ( ), and the speed S of shock waves to 

the slope of the connecting line between Q e (p-)) and 

(p+,Qe( P+ )) ("-). 

Based on this information, we can understand the 
phase diagram for the TASEP (see Fig. 33 and 
Sees. I.E.I, III.C, III.C.3). Stimulated by a study by 
Krug (1991), the corresponding results were obtained 
by Schiitz and Domany (1993) and Derrida et al. 
(1993). A detailled overview is given by Schiitz (2000). 
For a system with open boundary conditions with 
average occupation numbers (no) at the upstream 
end and (ul+i) at the downstream end, they found 
the lower density to dominate within the system, if 
(n ) = q /q < 1/2 and (n > = q /q < qL+i/q = 
(1 — (n^+i)), while the higher density dominated in the 
system, if (n ) = q /q > qL+i/l = (1 - and 
(1 — (n^+i)) = qL+i/q < 1/2. This is a consequence of 
downstream shock propagation in the first case, while 
shocks propagate upstream in the second case. When 
the line (n ) = q /q = QL+i/q = (1 - is 
crossed, a phase transition from a low-density state 
to a high-density state occurs. Finally, under the con- 
dition (n ) = qo/Q > 1/2 and qL+i/q > 1/2, i.e., 
( n L+i) < 1/2, the system takes on a maximum flow 
state. For a more detailled analysis see the review by 
Schiitz (2000). Related studies have been also carried 
out by Nagatani (1995c). 



FIG. 33. Phase diagram of the TASEP with open bound- 
aries and other models that can be macroscopically de- 
scribed by the Lighthill-Whitham or Burgers equation, 
with representative density profiles (see insets). The di- 
agonal line separates the areas A\ and An of downstream 
shock propagation from the areas B\ and Bu of upstream 
shock propagation with corresponding changes in the den- 
sity profile. Consequently, crossing this line is related to a 
transition between a low-density phase controlled by the 
upstream boundary and a high-density phase controlled 
by the downstream boundary. Area C is characterized by 
the maximum possible flow. (After Schiitz and Domany, 
1993; Schiitz, 2000.) 



d*(x,t) _ n d 2 *(x,t) 

~~dt~^ U dx 2 (55j 

by means of the Cole-Hopf transformation C(x, t) = 
-[2D/V(x,t)]dV(x,t)/dx (Whitham, 1974). 

The Burgers equation overcomes the problem 
of shock waves, but it cannot explain the self- 
organization of phantom traffic jams or stop-and-go 
waves (cf. Fig. 34). Recently, the main research activ- 
ity focusses on the noisy Burgers equation containing 
an additional fluctuation term (Forster et al., 1977; 
Musha and Higuchi, 1978; Nagel, 1995, 1996). 



2. The Burgers equation 



3. Payne's model and its variants 



The Lighthill-Whitham model is very instructive 
and the basis of a very elaborate shock wave the- 
ory (Haight, 1963; Whitham, 1974; Gazis, 1974; Hida, 
1979; Newell, 1993; Daganzo, 1999b), but the de- 
velopment of shock waves results in serious difficul- 
ties to solve the Lighthill-Whitham model numeri- 
cally. A suitable integration method is the Godunov 
scheme (Godunov, 1959; Ansorge, 1990; Bui et al, 
1992; LeVeque, 1992; Lebaque, 1995). As an alterna- 
tive, Newell (1993), Daganzo (1994, 1995a, b), and 
Lebacque (1997) have recently developed variants of 
the Lighthill-Whitham model. 

To avoid the development of shock fronts, one of- 
ten adds a small diffusion term to the LW model to 
smoothcn the wave fronts. Whitham himself (1974) 
has suggested to generalize Eq. (47) according to 
Q = [Q e {p) - Ddp/dx] or 

V(x,t) = V e (p(x,t))- ^J^ . (52) 



p{x,t) dx 



The resulting equation reads 



dp 
dt 
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The last term is a diffusion term. For a diffusion con- 
stant D > 0, it yields a significant contribution only 
when the curvature d 2 p/dx 2 of the density is great. In 
case of a density maximum, the curvature is negative, 
which results in a reduction of the density. A density 
minimum is smoothed out for analogous reasons. 

For the following analytical considerations, we as- 
sume the linear velocity-density relation (48). With 
the resulting propagation speed C(x,t) = Vo[l — 
2p(x,t)/pj am ] of kinematic waves, we can transform 
(53) into 



dC ^ + C(x,t) dC ^ 



dt 



dx 



D 



d 2 C(x,t) 
dx 2 



(54) 



This so-called Burgers equation is the simplest equa- 
tion containing non-linear propagation and diffusion. 
Surprisingly, this equation can be solved exactly, be- 
cause it is related to the linear heat equation 



Payne (1971, 1979a) complemented the continuity 
equation (45) by a dynamic velocity equation, which 
he derived from Newell's car-following model (21) by 
means of a Taylor approximation. Recently, related 
approaches have been pursued by Helbing (1998a), 
Nelson (2000), as well as Berg et al. (2000). 

Payne identified microscopic and macroscopic ve- 
locities according to v a (t + At) = V(x + V At, t + 
At) w [V (x, t) + V At dV(x, t) /dx + At dV{x, t) /dt] . 
Moreover, he replaced the inverse of the headway d a 
to the car in front by the density p at the place 
x + d a (t)/2 in the middle between the leading and 
the following vehicle: l/d a (t) = p(x + d a (t)/2,t) = 
p(x + l/(2p),f) w [p{x,t) + \/{2p)dp{x,t)/dx\. This 
yields v' e (d a (t)) = V e (l/d a (t)) « {V e (p{x,t)) + 
l/[2p(x,t)][dV e (p)/dp]dp(x,t)/dx}. From the previ- 
ous equations we can obtain Payne's velocity equation 
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p dx 



, (56) 



where we have introduced the density-dependent dif- 
fusion D(p) = -0.5dV e (p)/dp = 0.5\dV e (p)/dp\ > 0. 
The single terms of Eq. (56) have the following inter- 
pretation: 

(i) Like in hydrodynamics, the term VdV/dx is called 
the transport or convection term. It describes a mo- 
tion of the velocity profile with the vehicles. 

(ii) The term —[D(p)/(p At)]dp/dx is called antici- 
pation term, since it reflects the reaction of identical 
drivers to the traffic situation in their surrounding, in 
particular in front of them. 

(iii) The relaxation term [V e (p) — V]/At describes 
the adaptation of the average velocity V(x, t) to 
the density-dependent equilibrium velocity V e (p) . This 
adaptation is exponential in time with relaxation time 
At. 

The instability condition of Payne's model agrees 
exactly with the one of the optimal velocity model, if 
we set r = At (see Sees. III.A.2 and IV.A.l). In the 
limit At — ► 0, one can derive the stable model (52) 
as adiabatic approximation of Eq. (56). Inserting this 
into the continuity equation, we obtain Eq. (53) with 
Payne's density-dependent diffusion function. 



For reasons of numerical robustness, Payne (1979a, 
b) used a special discretization modifying his original 
model by numerical viscosity terms. His model was 
taken up by many people, for example, Papageorgiou 
(1983). Cremer and coworkers multiplied the antici- 
pation term by a factor p/(p + p ) with some con- 
stant po > to get a more realistic behavior at small 
vehicle densities p w (Cremer, 1979; Cremer and 
Papageorgiou, 1981; Cremer and May, 1986; Cremer 
and Mcifincr, 1993). A discrete version of this model is 
used in the simulation package SIMONE (Meifiner and 
Boker, 1996). Smulders (1986, 1987, 1989) introduced 
additional fluctuation terms in the continuity and ve- 
locity equations to reflect the observed variations of 
the macroscopic traffic quantities. Further modifica- 
tions of Payne's freeway simulation package FREFLO 
(1979b) were suggested by Rathi et al. (1987). 

4. The models by Prigogine and Phillips 

An alternative model was proposed by Phillips 
(1979a, b), who derived it from a modified ver- 
sion of Prigogine's Boltzmann-like traffic model (see 
Sec. III.E.l). Similar to Prigogine and Herman (1971), 
he obtained the continuity equation (45) and the ve- 
locity equation 

dV ir 8V ldP 1 , T71 

7H +v ^ = --p^ + ^ v ^- v] > (57) 

but with another density-dependent relaxation time 
r(p). The quantity P(x,t) = p(x,t)9(x,t) denotes the 
so-called "traffic pressure", where 9(x,t) is the veloc- 
ity variance of differently driving vehicles. For this 
variance, one can either derive an approximate dy- 
namic equation or assume a variance-density relation 
such as 9(x,t) — 9 [1 — p(x , t) / Pj am ] > (Phillips, 
1979a, b). According to this, the variance decreases 
with increasing density and vanishes together with the 
equilibrium velocity V e (p) at the jam density Pj a m, as 
expected. 

Like the Payne model, the model by Phillips is un- 
stable in a certain density range. Hence, it could ex- 
plain emergent stop-and-go waves, but it is not numer- 
ically robust. In particular, at high densities p the traf- 
fic pressure decreases with p, so that vehicles would 
accelerate into congested regions, which is unrealistic. 

FIG. 34. Illustration of the typical traffic dynamics ac- 
cording to macroscopic models with a dynamic velocity 
equation. The simulations were carried out with the non- 
local GKT model introduced in Sec. III.E.4 and show the 
spatio-temporal evolution of the traffic density p(x,t) on 
a circular road of circumference 10 km, starting with ho- 
mogeneous traffic, to which a localized initial perturbation 
(108) of amplitude Ap = 10 vehicles/km/lane was added, 
(a) Free and stable traffic at low vehicle density, (b) For- 
mation of stop-and-go waves in the range of linearly unsta- 
ble traffic at medium vehicle density. (From Treiber et al, 
1999; Helbing et al, 2001a, b.) 



5. The models by Whitham, Kiihne, Kerner, Konhauser, and 
Lee 

With relation (47), the Lighthill- Whitham model 
assumes that the traffic flow is always in equilibrium 
Q e {p)- This is at least very questionable for medium 
densities (Kerner and Rchborn, 1996b; Kerner et a/., 
1997), where there is no unique empirical relation be- 
tween flow and density (see Sees. II. B and IV. A. 3). 
Although more conservative researchers still believe 
that the Lighthill- Whitham theory is correct, imply- 
ing that there would be no unstable traffic in which 
disturbances are amplified, Daganzo (1999a) recently 
notes that "large oscillations in flow, speed and cumu- 
lative count increase in amplitude across the detectors 
spanning a long freeway queue and its intervening on- 
ramps", cf. Figs. 23 and 42. Moreover, Cassidy and 
Bcrtini (1999) state that "the discharge flows in ac- 
tive bottlenecks exhibit near-stationary patterns that 
(slowly) alternate about a constant rate .... The onset 
of upstream queueing was always accompanied by an 
especially low discharge flow followed by a recovery 
rate and these are the effects of driver behavior we do 
not yet understand" , cf. Figs. 3b and 39d. 

Whitham himself (1974) suggested a generalization 
of the LW theory, which is related to Phillip's one, 
but without the mentioned problem of decreasing traf- 
fic pressure. He obtained it from the continuity equa- 
tion together with the assumption that the adaptation 
to the velocity (52) can be delayed by some relax- 
ation time r: d v V/dt = [V e (p) - (D/p)dp/dx - V]/t. 
Therefore, Whithams velocity equation reads dV/dt + 
VdV/dx = -(9 /p)dp/dx + [V e {p) - V]/t with suit- 
able constants # = D > and r > 0. 

Unfortunately, this model shows the same problems 
with numerical robustness as the models by Payne and 
Phillips, since it eventually produces shock-like waves. 
Therefore, Kiihne (1984a, b, 1987) added a viscosity 
term vd 2 V '/ 'dx 2 , which has a similar effect as the dif- 
fusion term in the Burgers equation and resolves the 
problem (cf. Fig. 34). In analogy with the Navier- 
Stokcs equation for compressible fluids, Kerner and 
Konhauser (1993) preferred a viscosity term of the 
form (i]/p)d 2 V/dx 2 , leading to the velocity equation 

dV | ydV = 9 dp | V d 2 V | V e (p)-V (5g) 
dt dx p dx p dx 2 t 

A stochastic version of this KK model has been de- 
veloped by Kerner et al. (1995b) as well to consider 
fluctuations in vehicle acceleration (Kiihne, 1987). 

It turns out that the above KK model is rather sen- 
sitive to the choice of the parameters and the velocity- 
density relation. This is probably for similar reasons 
as in the optimal velocity model, which had also no 
dependence of the equilibrium velocity on the relative 
velocity between cars. The most suitable specification 
of parameters was probably made by Lee et al. (1998). 



6. The Weidlich-Hilliges model 

The basic version of the model by Weidlich and 
Hilliges (Weidlich, 1992; Hilliges, Reiner and Weidlich, 
1993; Hilliges and Weidlich, 1995) is discrete: The free- 
way is divided into cells j of equal length Ax, which 
reminds of cellular automata models. Moreover, the 
density p(j, t) = p(j Ax, t) is governed by the spatially 
discretized continuity equation 

dp(j,t) Q(j,t)-Q(J-l,t) = 
dt Ax ' 

However, Weidlich and Hilliges assume for the flow 
the phenomenologically motivated relation 

Q(j,t)= P (j,t)V(j + l,t), (60) 

according to which the drivers in cell j adapt their 
speed to the velocity in the next cell (j + 1). This 
looks like a generalization of the mean value equa- 
tions of the TASEP (see Sec. III.C.3) and reflects 
an anticipatory behavior of the drivers. The choice 
Ax Ri 100 m of the cell length is, therefore, deter- 
mined by the driver behavior. Under the assump- 
tion V(j, t) = V e (p(j,t)), the model cannot de- 
scribe emergent stop-and-go waves, since the con- 
tinuous version of the Weidlich-Hilliges model corre- 
sponds to Eq. (53) with the diffusion function D{p) = 
\V e (p) — pdV e /dp]Ax/2 > 0, as a Taylor expansion of 
p((j - \)Ax,t) and V((j + l)Ax,t) shows. For the 
linear velocity-density relation (48), we have simply 
a diffusion constant D = Vq Ax/2, which implies an 
equivalence to Eq. (53) and the Burgers equation. 

To model self-organized stop-and-go waves, Hilliges 
and Weidlich (1995) have supplemented Eqs. (59), 
(60) by the dynamic velocity equation 

9VU,t) , in , ., V(j + l,t)-V(j-l,t) 
dt +V V> 1) 2 Ax 

= Ve( P (j,t))-V(3,t) 
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with r w 5 s. In the continuous limit Ax — > 0, the 
terms on the left-hand side correspond to the sub- 
stantial time derivative d v V/dt = dV/dt + VdV/dx. 

Moreover, Hilliges and Weidlich (1995) have shown 
that dp(j, t)/dt > if p(j, t) = 0, and dV(j, t)/dt > 
if V(j, t) = 0. This guarantees the required non- 
negativity of the density and the average velocity. Ad- 
ditionally, they have developed simple rules for the 
treatment of junctions, which allow the simulation of 
large freeway networks. 

7. Common structure of macroscopic traffic models 

The previous models are closely related with each 
other (Helbing, 1997a, e), as they can all be viewed as 
special cases of the density equation 
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and the velocity equation 
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The differences are in the specification of the diffu- 
sion D{p), the fluctuations £i(x,i), & t), the traffic 
pressure P(p), the viscosity-like quantity f(p), the re- 
laxation time r(p), and the equilibrium velocity V e (p). 
The Lighthill-Whitham model, for example, results in 
the limit r — > 0. Payne's model is gained with r = At 
and P(p) = [Vq - V e (p)]/(2At), while P(p) = P 6 in 
the KK model (58). Some parameters must be set to 
zero in an obvious way. 

Section III.E will show how to derive macroscopic 
from microscopic traffic models. We will see that D 
and £i should be non- vanishing only in "first-order 
models" restricting themselves to a density equation. 
In "second-order models" with a dynamic velocity 
equation, the pressure P consists of a contribution 
p9, if there is a variation of individual velocities 
due to fluctuations or different driver- vehicle types, 
and other contributions which are due to the non- 
locality of vehicle interactions (Helbing 1996b, 1997a, 
1998a, d; see also Sees. III.E.5 and III.D.3). Daganzo 
(1995c) has fundamentally criticized second-order 
macroscopic traffic models, but his arguments have 
been invalidated (Whitham, 1974; Paveri-Fontana, 
1975; Helbing, 1995b, 1996a, b, 1997a; Helbing and 
Treiber, 1999; Treiber et at., 1999; Aw and Rascle, 
2000; Zhang, 2000). 

E. Gas-kinetic traffic models and micro-macro link 

In the following paragraphs, I will show how the car- 
following models of Sec. III. A can be connected with 
macroscopic traffic models via a mesoscopic level of 
description. 

1. Prigogine's Boltzmann-like model 

Starting in 1960, Prigogine and coworkers proposed, 
improved, and investigated a simple gas-kinetic model 
(Prigogine and Andrews, 1960; Prigogine, 1961). The 
results are summarized in the book by Prigogine and 
Herman (1971), to which I refer in the following. 

Gas-kinetic theories are based on an equation for 
the phase-space density 

p(x,v,t)= p(x,t)P(v;x,t), (64) 

which is the product of the vehicle density p(x, t) and 
the distribution P(v;x,t) of vehicle speeds v at loca- 
tion x and time t. Because of vehicle conservation, we 



find again kind of a continuity equation for the phase- 
space density: 



dt 



dp 

at 



dp 
dx 



(65) 



int 



However, compared to Eq. (45), the right-hand side 
of this conservation equation is not zero this time, be- 
cause of velocity changes. In a coordinate system mov- 
ing with speed v, the temporal change of the phase- 
space density is given by the acceleration behavior and 
interactions of the vehicles. 

Prigogine suggested that the acceleration behavior 
can be described by a relaxation of the velocity distri- 
bution P(v;x, t) to some desired distribution Pq(v): 



PpM) 
r(p(x,t)) 



[P (v)-P(v;x,t)}. (66) 



The quantity r(p) denotes again a density-dependent 
relaxation time. By Po(v), the model intends to reflect 
the variation of the desired velocities among drivers, 
causing the scattering of the actual vehicle velocities 
v. The distribution Po(v) can be obtained by mea- 
suring the velocity distribution of vehicles with large 
clearances (Tilch, 2001). 

Interactions among the vehicles are modelled by the 
Boltzmann-likc term 



dw [1 — p(p)] \w — v\p(x, w, t)p(x, v, t) 

J dw [1 — p(p)] \v — w\p(x, v, t)p(x, to, t) . 

(67) 



This assumes that fast vehicles with velocity w in- 
teract with slower ones with velocity v < w at a 
rate \w — v\p(x,w,t)p(x,v,t), which is proportional 
to the relative velocity \w — v\ and to the product 
of the phase-space densities of the interacting vehi- 
cles, describing how often vehicles with velocities w 
and v meet at place x. Given that the faster vehicle 
can overtake with some density-dependent probability 
p(p), it will have to decelerate to the velocity v of the 
slower vehicle with probability [1 — p{p)\, increasing 
the phase-space density p(x, v, t) accordingly. How- 
ever, the phase-space density is decreased when ve- 
hicles of velocity v meet slower vehicles with velocity 
w < v. This is reflected by the last term of Eq. (67). 
The interaction term (67) can be simplified to 



= [l-p(p(x,t))]p{x,t)[V(x,t) - v]p(x,v,t) 



int 



(68) 



by introducing the average velocity 



V{x,t) = [ dvvP{v;x,t) = f*»^M. (69) 



Additionally, we define the velocity variance 

9(x,t)= J dv[v-V(x,t)] 2 P(v;x,t) (70) 

and the average desired velocity Vq — J dvP (v). To 
derive macroscopic density and velocity equations, we 
multiply the kinetic equation given by Eqs. (65), (66), 
and (68) by 1 or v and integrate over v, taking into 
account that the vehicle density is given by 

(71) 



p(x,t) = J dvp(x,v,t). 

We obtain equations (62) and (63) with D(p) = = 
i/(p), £i(t) = = £2(^)1 the traffic pressure 

P(p,e)=p9, (72) 

and the equilibrium velocity 

V e (p,6) = V -T(p)[l-p(p)}p6. (73) 

Therefore, the gas-kinetic approach allowed Prigogine 
and his coworkers to derive mathematical relations for 
the model functions P and V e , which other researchers 
had to guess by phenomenological considerations be- 
fore. Unfortunately, the above relations are only valid 
for small densities, but theories for higher densities are 
meanwhile available (see Sees. III.E.3 and III.E.4). 

Analyzing the Boltzmann-like traffic model, Pri- 
gogine (1961) found a transition from free to con- 
gested traffic above a certain critical density, which 
he compared to the phase transition from a gaseous to 
a fluid phase (see also Prigogine and Herman, 1971). 
This congested state is characterized by the appear- 
ance of a second maximum of the velocity distribution 
at v = 0. That is, part of the vehicles move, while 
others are completely at rest. This bimodal equilib- 
rium distribution is also found in modern variants of 
Prigogine's model (Nelson, 1995; Nelson et al., 1997), 
which is a consequence of the special acceleration term 
(66). 

To obtain a better agreement with empirical data, 
Phillips (1977, 1979a, b) modified Prigogine's rela- 
tions for the overtaking probability p(p) and the re- 
laxation time r{p). Moreover, he assumed that the 
distribution of desired velocities Pq(v) depends on 
the density p and the average velocity V. Finally, he 
made several proposals for the velocity variance 9. Pri- 
gogine's model has been also improved by Andrews 
(1970a, b, 1973a, b), Phillips (1977, 1979a), Paveri- 
Fontana (1975), Nelson (1995) and others (Helbing, 
1995c, d, 1996a, b, 1997a, d, 1998a; Wagner et al, 
1996; Wagner, 1997a, b, c, 1998b; Nagatani, 1996a, 
b, c, 1997a, b; Wegener and Klar, 1996; Klar et al, 
1996; Klar and Wegener, 1997, 1999a, b; Shvetsov 
and Helbing, 1999). They modified the acceleration 
term, introduced velocity correlations among succes- 
sive vehicles, investigated interactions among neigh- 
boring lanes, and/or considered space requirements. 



2. Paveri-Fontana's model 



3. Construction of a micro-macro link 



Pavcri-Fontana (1975) carried out a very detaillcd 
investigation of Prigogine's gas-kinetic traffic model 
and recognized some strange features. He critized that 
the acceleration term (66) would describe discontin- 
uous velocity jumps occuring with a rate related to 
t(p). Moreover, the desired velocities of vehicles were 
a property of the road and not the drivers, while 
there are actually different personalities of drivers 
(Daganzo, 1995c): aggressive ones, driving fast, and 
timid ones, driving slowly. 

Paveri-Fontana resolved these problems by distin- 
guishing different driver- vehicle types. These were 
characterized by individual desired velocities vq- As a 
consequence, he introduced an extended phase-space 
density p(x, v, t>n, i) and the corresponding gas-kinetic 
equation. By integration over vq, one can regain equa- 
tion (65). While Paveri-Fontata specified the interac- 
tion term analogously to Prigogine's one, a different 
approach was chosen for the acceleration term. It cor- 
responds to the microscopic acceleration law 



dv/dt = (vq — v)/r(p) 
finally resulting in the formula 
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dv 



p(x,v,t) 



V (v;x,t) - v 
r{p(x,t)) 



(74) 



(75) 



(Helbing, 1995d, 1996a; Wagner et al, 1996). Herein, 
the 

quantity Vo(v;x,t) — J dv v p(x, v, Vq, t)/p(x, v, t) w 
V (x,t) + [v- V(x,t)]C'(x,t)/6(x,t) denotes the av- 
erage desired velocity of vehicles moving with velocity 
v, which is greater for fast vehicles. While Vo(x, t) = 
J dvo J dvvop(x,v,vo,t)/p(x,t) represents the aver- 
age desired velocity at place x and time t, C'(x,t) = 
J dv J dv(v - V)(v - V )p(x,v,v ,t)/p(x,t) repre- 
sents the covariance between the actual and desired 
velocities. In spite of these differences, the macro- 
scopic equations for the density and the average veloc- 
ity agree exactly with the ones by Prigogine. However, 
the equations for the velocity variance, the covariance, 
the average desired velocity and so on are different. 

The existence, construction, and properties of solu- 
tions of the Paveri-Fontana equation have been care- 
fully studied (Barone, 1981; Semenzato, 1981a, b). In 
addition, it should be noted that there are further al- 
ternatives to the specification (75) of the acceleration 
term. For example, Alberti and Belli (1978) assume 
a density-dependent driver behavior. In contrast, Nel- 
son (1995) and others (Wegener and Klar, 1996; Klar 
and Wegener, 1997) have proposed an acceleration in- 
teraction leading to mathematical expressions similar 
to Prigogine's interaction term. 



For the derivation of macroscopic equations, one 
may start with the master equation (32) together with 
the transition rates (37), if the state j is replaced by 
the state vector (x,v), and the sums are replaced by 
integrals due to the continuity of the phase space. The 
spontaneous state changes are described by the tran- 
sition rate 



w(x',v'\x, v) — £ lm Q ~^.^( x ' ~ ( x + v ^)) 



x / d£ 



1 



-(£At/T) 2 /(2DAt) 



xs(v'-(v+ V ^±^Atfj, (76) 

reflecting the motion and acceleration of identical 
vehicles with Gaussian-distributed velocity fluctua- 
tions. Inserting this into Eq. (40) and reformulat- 
ing the respective terms as a Fokker-Planck equa- 
tion (see Sec. III.C.5) gives, together with p(x,v,t) = 
( n (x,v))/Ax, the following kinetic equation: 



dp d(pv) d / v - v 
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(77) 



Herein, the quantity D represents a velocity diffusion 
coefficient, and 

W,„rJ d *'J dv 'I dv J dw J dy 'J du ' 

x [w 2 (x, v, y, w\x' , v 1 , y' , w r )p 2 (x', v', y' , w' , t) 
-w 2 (x', v', y', w'\x, v, y, w)p 2 (x, v, y, w, t)] (78) 

is the interaction term, where we have in- 
troduced the so-called pair-distribution function 
p 2 (x,v,y,w,t) = ((i( 1; „)jj (9 ^))/(Ai) 2 . The transi- 
tion rates w 2 {x' , v', y', w'\x, v, y, w) are proportional 
to the relative velocity \v — w\ and proportional to the 
probability of transitions from the states (x, v) and 
(y,w) of two interacting vehicles to the states (x',v') 
and (y 1 , w'). During the interactions, the car locations 
do not change significantly, only the velocities. There- 
fore, we have 

w 2 (x',v',y',w'\x,v,y,w) 
= lim -r-fl — p(p)] \v — io I 

x6(v' -[v + f(x, v, y, w) At])S(x' - x) 

xS(w' -[w + f(y, w, x, v) At])S(y' - y) (79) 

(Helbing, 1997a). Herein, the force / is specified in 
agreement with Eq. (16). The location and velocity of 
the interaction partner are represented by y and w, 
respectively. Quantities without and with a prime (') 
represent the respective values before and after the 



interaction. To simplify the equations, we will assume 
hard-core interactions corresponding to abrupt brak- 
ing maneuvers, when neighboring particles approach 
each other and have a distance d* , i.e. y = (x ± d*). 
Moreover, let us assume the relations (v' — w') = 
— e(l — p)(v — w) and (v' + w') = (1 — p)v + (1 + p)w 
corresponding to 



, l-n-e{l-n) 1 + n + e{l- fi) 
v = v + w, 

w , = l-fJ. + e(l-fi) ^ 1 + n - g(l - m) 



w . (80) 



By variation of the single parameter p reflecting the 
asymmetry of interactions, this allows us to switch 
between vehicle interactions for p, = 1 (velocity adap- 
tation) and granular interactions for /i = (conserva- 
tion of momentum, see Foerster ei aL, 1994). In addi- 
tion, the parameter e allows us to control the degree 
of energy dissipation in granular interactions, where 
e(l — /i) = 1 corresponds to the elastic case with en- 
ergy conservation relevant for idealized fluids. 

Prigoginc's Boltzmann-likc interaction term re- 
sults for p = 1, d* = (point-like vehicles), 
and the so-called vehicular chaos approximation 
p 2 (x,v,y,w,t) s=s p(x,v,t)p(y,w,t), which assumes 
statistical independence of the vehicle velocities v and 
w. This is only justified in free traffic. At higher den- 
sities, we must take into account the increase of the 
interaction rate by a density-dependent factor \ > 1 
(Chapman and Cowling, 1939; Lun et al., 1984; Jenk- 
ins and Richman, 1985), since the proportion 1/x > 
of free space is reduced by the vehicular space require- 
ments (Helbing, 1995b, d, 1996a, b; Wagner et al., 
1996; Klar and Wegener, 1997). For simplicity, we will 
assume that the overtaking probability is just given by 
this proportion of free space, i.e. 



V = i/x 



(81) 



(Trciber et al, 1999). At high vehicle densities, the 
velocities of successive cars will additionally be cor- 
related. This is particularly clear for vehicle pla- 
toons (Islam and Consul, 1991), in which the varia- 
tion of vehicle velocities is considerably reduced (see 
Sec. II. D). For the treatment of vehicle platoons there 
have been several suggestions (Andrews, 1970a, b, 
1973a, b; Beylich, 1979, 1981). Related to an ap- 
proach by Lampis (1978) and others (Poethke, 1982, 
see Leutzbach, 1988; Hoogendoorn and Bovy, 1998a, 
b; Hoogendoorn, 1999), one possibility is to distin- 
guish a fraction c(p) of free vehicles and a fraction 
[1 — c(p)] of bound (i.e. obstructed) vehicles (Helbing, 
1997a, d). The most realistic description of vehicle 
platoons, however, is to set up equations for vehicle 
clusters of different sizes (Beylich, 1979, 1981). Un- 
fortunately, this ends up with a very complex hierar- 
chy of spatio-temporal equations, from which useful 
macroscopic equations have not yet been successfully 
derived. Instead, we may simply assume a correlation 



r of successive vehicle velocities, which increases from 
zero to positive values when the density changes from 
low to high densities (see Fig. 18). Taking also into 
account the fact that vehicle velocities are more or 
less Gaussian distributed (cf. Figs. 13 and 35; Pam- 
pel, 1955; Gerlough and Huber, 1975; May, 1990; see 
also Alvarez et al., 1990), one may approximate the 
two-vehicle velocity-distribution function by a bivari- 
ate Gaussian function 



P 2 {v,w) = 



Vdet R 
2tt 



-R(v,w)/2 



(82) 



with the quadratic form R(v, w) = [(v— V) 2 /9— 2r(v — 
V)(w- V + )/ x /99\+{w-V + ) 2 /9 + ]/{l-r 2 ) and its 
determinant det R = l/[99 {I - r 2 )] (Shvetsov and 
Helbing, 1999). Herein, V denotes the average velocity 
and 9 the velocity variance of vehicles at location x 
and time t, whereas the variables with the subscript 
"+" refer to the location x + — [x + d*(V)] of the 
interaction point. Often, one chooses the interaction 
distance 



d*{V)= 1 



1 



TV 



(83) 



where 7 > 1 is an anticipation factor. I should mention 
that an analytical solution of the gas-kinetic equa- 
tions is not available, even for the stationary case, 
but the above formula seems to be a good approxima- 
tion of numerical results (Wagner, 1997a). The bivari- 
ate Gaussian distribution is favourable for the evalu- 
ation of the interaction integrals, and it makes sense 
for 3\/9 < V, so that negative velocities are negligible. 

FIG. 35. Stationary velocity distributions at different ve- 
hicle densities obtained by numerical solution of a gas- 
kinetic traffic model (after Wagner, 1997a). 

Summarizing the above considerations, a good ap- 
proximation for the pair- distribution function p~2 at 
high vehicle densities is given by 

p 2 (x,v,x + ,w,t) = x + p(x,t)p(x + ,t)P 2 (v,w;x,x + ,t) , 

(84) 

and the resulting Enskog-like gas-kinetic traffic equa- 
tion reads 



dp d(pv) d (JVq—v 



dt dx dv v 

+ (1 -p+)x + pp+ 



d 2 ~ 
6#W 



J "dw(w — v)P 2 (w, v) 



— J dw (v — w)P 2 (v, w) 



(85) 



4. The non-local, gas-kinetic-based traffic model (GKT 
model) 

The probably most advanced macroscopic traffic 
models presently available is obtained by multiplying 



equation (77) with v k (k £ {0,1,2}) and subsequent 
integration over v. The corresponding calculations are 
more or less straightforward, but very cumbersome, 
so the results were only eventually gained (Hclbing, 
1995d, 1996a, b, 1997a, 1998a; Helbing and Treiber, 
1998a; Shvetsov and Helbing, 1999). They can be cast 
into the general form of the continuity equation (62) 
and the velocity equation (63), but instead of V e we 
have a non-local expression 

V e = V - t [1 - p( P+ )}x( P+ )p + B(AV, 8 A ) . (86) 

S v ' 

Braking Interaction Term 

This non-locality has some effects which other macro- 
scopic models produce by their pressure and viscosity 
(or diffusion) terms. Related to this is the higher nu- 
merical efficiency of the non-local GKT model used 
in the simulation package MASTER (Hclbing and 
Treiber, 1999). The dependence of the braking in- 
teraction on the effective dimcnsionless velocity dif- 
ference AV = (V — V + )/y/02 between the actual 
position x and the interaction point x + takes into 
account effects of the velocity variances 9, 9, and 
the correlation r among successive cars (cf. Fig. 18): 
9 A = {9- 2r^fMl+9 + ). The "Boltzmann factor" 

B(AV,9 & ) - 9 & {AV^(AV) + [1 + (AV) 2 ]E{AV)} 

(87) 

in the braking term is monotonically increasing with 
AV and contains the normal distribution $(z) = 
qT z I 2 /^/2tt and the Gaussian error function E{z) = 

J dz'^(z'). 

— oo 

To close the system of equations, we specify the 
velocity correlation r as a function of the density ac- 
cording to empirical observations (see Fig. 18). More- 
over, for a description of the presently known proper- 
ties of traffic flows it seems sufficient to describe the 
variance by the equilibrium relation 9 = Dt of the 
dynamical variance equation (Helbing, 1996b, 1997a), 
which also determines the traffic pressure via P = p9 
or, more exactly, Eq. (119). Empirical data (Helbing, 
1996b, 1997a, c; Treiber et al, 1999) suggest that the 
variance is a density-dependent fraction A(p) of the 
squared average velocity: 



8 = A(p)V 2 . 



(88) 



This guarantees that the velocity variance will vanish 
if the average velocity goes to zero, as required, but 
it will be finite otherwise. The variance prefactor A 
is higher in congested traffic than in free traffic (see 
Fig. 17). 

Finally, the average time clearance of vehicles 
should correspond to the safe time clearance T, if the 
vehicle density is high. Hence, in dense and homoge- 
neous traffic situations with V + = V and 9 + = 9, 
we should have V = (1/p— l/p m ax)/7\ This is to be 
identified with the equilibrium velocity 



V e (p)^l-1 + Jl + 'S\ . ,,0, 



where V 2 = V~o/[TpA(p)(l — p)x\. For the "effective 
cross section" we, therefore, obtain 



[i-p(p)]x(p) = 



V pT 2 



-A(p max )(l - p/p max ) 2 



(90) 



which makes also sense in the low-density limit p — > 0, 
since it implies x — > 1 and p — ► 1 . 



5. Navier-Stokes-like traffic equations 

The above derivation of macroscopic traffic equa- 
tions was carried out up to Euler order only, based 
on the (zeroth order) approximation of local equilib- 
rium. This assumes that the form po(x,v,t) of the 
phase-space density in equilibrium remains unchanged 
in dynamic situations, but it is given by the local val- 
ues of the density p(x, t) and possibly other macro- 
scopic variables, see e.g. Eq. (93). For fluids and gran- 
ular media, however, it is known that the phase-space 
density depends also on their gradients, which is a 
consequence of the finite adaptation time To needed 
to reach local equilibrium (Huang, 1987). The related 
correction terms lead to so-called Navier-Stokes equa- 
tions, which have been also derived for gas-kinetic 
traffic models (Helbing, 1996b, 1998d; Wagner, 1997b; 
Klar and Wegener, 1999a, b). For ordinary gases, the 
Navier-Stokes terms (transport terms) are calculated 
from the kinetic equation by means of the Chapman- 
Enskog method (Enskog, 1917; Chapman and Cowl- 
ing, 1939; Liboff, 1990). Here, I will sketch the more in- 
tuitive relaxation-time approximation (Huang, 1987; 
for details see Helbing and Treiber, 1998c). This as- 
sumes that 

(i) the deviation 5p(x, v, t) — [p(x, v, t) — po(x, v, t)] is 
small compared to po(x,v,t), so that 
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(91) 



with^/dt = d/dt+vd/dx = d/dt+V e d/dx+Svd/dx 
and 5v(x, t) = [v — V e (x, t)], 

(ii) the effect of the linear interaction operator C can 
be characterized by its slowest eigenvalue — I /To < 
(which depends on the density and velocity variance): 



C(8p) « - 



Sp(x, v, t) 



(92) 



In order to reach a time-scale separation, the equi- 
librium distribution is usually expressed in terms of 
the macroscopic variables with the slowest dynamics 



related to the conserved quantities. In fluids, these 
are the particle number, momentum, and kinetic en- 
ergy, while in traffic, this is the particle number only. 
Between the dynamics of the velocity moments (v k ) 
and (v k+1 ) with k > 1, there is no clear separation of 
time scales. Therefore, it is consequential to express 
the equilibrium phase-space density as a function of 
the vehicle density p(x, t) only: 



p (x,v,t) 



[2n0 e (p)]V* P l 29 e (p) /' 



(93) 



where e (p) = A(p)[V e (p)] 2 . The same is done with 
the pair distribution function Pz(v, w). This allows us 
to determine 



d v po _ dfodvP _ dpo ( dp + ^ dp ^dp 
dt dp dt dp \dt e dx dx 



(94) 



The relation dpo/dp can be simply obtained by 
differentiation of po, whereas dp/dt + V e dp/dx = 
—pdV e /dx follows from the continuity equation (45) 
in Euler approximation. Inserting (92) to (94) into 
(91) finally yields 



Sp(x,v,t) = -T \p 



Sv dV e 
e dp 
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26, 
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Tj~dp~ 



+ - 
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dpo 
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-Sv 



Integrating this over v gives exactly zero, i.e. there are 
no corrections to the vehicle density. However, multi- 
plying this with Sv and afterwards integrating over Sv 
yields corrections SV(x,t) = [V(x,t) — V e (p(x,t))] to 
the equilibrium velocity according to 



pSV 



-T aP 

V e -v 



9 e dp d9 e dp dV e dV e 
p dx dp dx dp ^ dx 



T 

d( P 8 e ) 
dx 



(1 -P+)X+P+B 



(95) 



where B is again the Boltzmann factor (87) . Interest- 
ingly enough, this Navicr-Stokes correction contains 
all terms and, therefore, captures all effects of the 
dynamic velocity equation (63) in Euler approxima- 
tion, which makes the success of the Euler-like GKT 
model understandable. Considering V + = V e (p + ) = 
V e (p(x + d,t)) w [V e (p(x,t)) +d(dV e /dp)dp/dx] and 
8 + w [8 e (p(x,t))+d(d9 e /dp)dp/dx], we can carry out 
the linear Taylor approximation 



:(v e -v e ) 



pd fdV e dV e dV e dV e d0 e \dp 



dp + dV + dp 
f d{ P 9 e ) dP cS \ dp 
\ dp dp J dx 



d9 + dp J dx 
(96) 



for the purpose of linear stability analysis. This defines 
the effective traffic pressure P c s which, apart from 
the kinematic dispersion effect p0, contains additional 
terms arising from vehicle interactions, as in Payne's 
model and in Sec. III.E.6. It can be shown that, in 
contrast to d(p6 e )/dp, dP e s / dp is non-negative for rea- 
sonably chosen parameter values and thereby resolves 
the problem of vehicle acceleration into congested ar- 
eas (see Helbing, 1996b, 1997a; Hclbing and Greiner, 
1997). Inserting Eqs. (95) and (96) into the continuity 
equation (45) results in an equation of type (53) with 
the diffusion function 



dP cS 
dp 



P -p) 



(97) 



This diffusion function in the Navier-Stokes-like traf- 
fic model with eliminated dynamic velocity equation 
becomes negative exactly when the Euler-like traffic 
model with dynamic velocity equation becomes lin- 
early unstable. Unfortunately, the resulting equation 
is even more difficult so solve numerically than the 
Lighthill-Whitham model. This problem can, however, 
be resolved by taking into account higher-order terms 
suppressing high frequency oscillations. 



6. Simultaneous micro- and macrosimulation 

We will now discuss a way of obtaining macroscopic 
from microscopic traffic models which is different from 
the gas-kinetic approach (see Hennecke et at, 2000; 
Hclbing et at, 2001b). For this, we define an average 
velocity by, for example, linear interpolation between 
the velocities of neighboring vehicles: 



V(x,t) = 



v a {t)[x a -\{t) -x]+ v a -i(t)[x - x a (t)\ 



x a -\{t) - x a {t) 



(98) 



{xa-i > x > x a ). While the derivative with respect 
to x gives us dV/dx = [v a -i(t) — v a (t)]/[x a -i(t) — 
x a (t)], the derivative with respect to t gives us the 
exact equation 



(99) 



where 



A(x,t) 



a' a {t)[x a - 1 {t)-x\+a' cc _ 1 {t)[x-x a {t)} 



x a -i{t) - x a (t) 



(100) 



is the linear interpolation of the vehicle accelerations 
a' a = dv a /dt characterizing the microscopic model. 
For most car-following models, the acceleration func- 
tion can be written in the form a' a — a'(v a ,Av a , s a ), 
where Av a — (v a — v a -\) is the approaching rate, 
and s a = {x a _\ — x a — l a -i) the netto distance to 



the vehicle in front. In order to obtain a macroscopic 
system of partial differential equations for the average 
velocity and density the arguments v a , Av a , and s a of 
the vehicle acceleration have to be expressed in terms 
of macroscopic fields, in particular the local vehicle 
density p(x, t) defined by 

1 _ 1 | s a {t)[xq-i{t)-x]+s a _ 1 (t){x-x a (t)} 
p(x,t) (t)-xjt) 

(101) 

Specifically, we make the following approxima- 
tions: A(x,t) w a'(V(x,t),AV(x,t),S(x,t)), with 
AV(x,t) = [V(x,t)-V + (x,t)l S(x,t) = \[p- 1 {x,t) + 
p^ 1 (x,t)] — p~*x) an( i tne non-locality given by 
9+(x,t) = g(x + l/p(x,t),t) with g e {p,V}. In 
this way, we obtain a non-local macroscopic velocity 
equation, which supplements the continuity equation 
(45) for the vehicle density and defines, for a given 
microscopic traffic model, a complementary macro- 
scopic model. Note that it does not contain a pres- 
sure term (l/p)dP/dx with P = p9, in contrast to 
the GKT model described above, since we have as- 
sumed identical vehicles and a deterministic dynamics 
corresponding to 9 — 0. For heterogeneous traffic, we 
expect the additional term — (1/ p)d(p9 e )/dx, where 
9 e (p) reflects the resulting, density-dependent veloc- 
ity variance of the vehicles. For the intelligent driver 
model (see Sec. III. A. 3), the above procedure leads 
to a remarkable agreement, at least for identical ve- 
hicles. It is even possible to carry out simultaneous 
micro- and macrosimulations of neighboring freeway 
sections without any significant changes in the density 
and velocity profiles or their propagation speed. This 
has been shown for downstream propagating pertur- 
bations, for traffic jams propagating upstream, and for 
complex spatio-temporal traffic patterns (Hcnnecke 
et al, 2000; Helbing et al, 2001b). 

7. Mesoscopic traffic models 

For the sake of completeness, I also mention some 
"mesoscopic" traffic models, which describe the micro- 
scopic vehicle dynamics as a function of macroscopic 
fields. Examples are the model by Wiedemann and 
Schwerdtfeger (1987) applied in the simulation tool 
DYNEMO (Schwerdtfeger, 1987), and the model by 
Kates et al. (1998) used in the simulation tool ANI- 
MAL. 

IV. PROPERTIES OF TRAFFIC MODELS 

In order to see which model ingredients are required 
to reproduce certain observations, we will start with 
simple one-lane models for identical vehicles, and add 
more and more aspects throughout the next sections. 



A. Identical vehicles on homogeneous freeways 

1. "Phantom traffic jams" and stop-and-go traffic 

For a long time, traffic research has aimed at re- 
producing the fundamental diagram and the kind of 
phantom traffic jams observed by Treiterer and oth- 
ers (see Sec. II.E.l). Beforehand, the fundamental dia- 
gram was normally fitted by calibration of the param- 
eters in the respective equilibrium flow-density rela- 
tion of the model under consideration. This procedure, 
however, makes only sense in the density regime(s) of 
stable traffic, as we will see that unstable traffic does 
not just oscillate around the fundamental diagram. 

Most traffic models have the following mechanism 
of traffic instability in common: 

(i) Overreaction of drivers, which is either reflected by 
a positive slow-down probability p or a delayed reac- 
tion (relaxation time r or reaction time At) . 

(ii) Chain reaction of followers: Before a decelerated 
vehicle manages to return to its previous speed, the 
next car approaches and has to brake as well (Nagcl 
and Paczuski, 1995). Overreacting followers are get- 
ting slower and slower, until traffic comes to a stand- 
still, although everyone likes to drive fast. 

The instability condition for the cellular automa- 
ton by Nagel and Schreckenberg with small slowdown 
probability p w follows directly from (ii). Starting 
with equally distributed vehicles on a circular road, a 
vehicle reaches the leader in one time step At, if the 
average clearance (L — N Ax)/N = (L/N — Ax) be- 
comes smaller than the average distance moved with 
the free velocity (v m & x — p)Ax/ At. Hence, the average 
vehicle density g = N / L above which jams are formed, 
is q w l/[(v max + l)Ax] (Nagel and Herrmann, 1993; 
Liibeck et al, 1998; Eisenblatter et al, 1998). A for- 
mula for finite slowdown probabilities p > 0, which 
takes into account that fluctuations may add up and 
initiate jam formation, has recently been derived by 
Gerwinski and Krug (1999). They calculated the reso- 
lution speed Co(p) = Co{p)Ax/ At of the downstream 
front of a "megajam" and obtained 

e= (102) 

[*max -p + C (p)]Ax 

with C (p) < (1 -p). 

For most traffic models, however, the instability of 
traffic flow is investigated by means of a linear stability 
analysis. This will be sketched in the following. 

For microscopic models of the form 

dv a (t + At) = v e (s a (t),v a (t),Av a (t)) -v a (t) 

dt T 

often with the simple specification v e (s,v,Av) = 
v e {s), we start our analysis by assuming a circular 
road of length L with homogeneous traffic flow of den- 
sity g. Homogeneity implies that all N = gL vehicles 
a have identical distances d = (s + 1) = 1/g to the re- 
spective leading vehicle, and the relative velocity Av a 



is zero. Hence, the situation is similar to a lattice of 
coupled particles in solid state physics, with the dif- 
ference that 

(i) the particles move with the identical velocity v = 
v e (d) along the trajectories x (0) — a(s + I) + vt = 
Xo(0) — ad + vt, 

(ii) the interactions are repulsive and anisotropic, and 

(iii) there is a driving force counteracting the repul- 
sion. 

The dynamics of this spatially homogeneous equilib- 
rium state in the presence of disturbances is, therefore, 
not obvious, but as usual we can carry out a stability 
analysis of the linearized model equations 



dSv a (t) t f5v a (t) 



dt 
dv 



dt 2 
dv 



1 f dv c 



[5x a _i{t) - 5x a (t)] 



dAv 



Due to linearization, this equation is only valid for 
small disturbances Sv a (t) = 5x a (t)/St = [v a (t) — v] <C 
v and Sx a (t) = x a (t) — [xn(0) — ad+vt] <C d. Its general 
solution has the form of a Fourier series: 

Sx a (t) = i ^ c fe cxp ( 27ri^ + (Afe - iu k )t) , 
fc=0 ^ ' 

where c k are the amplitudes of oscillations with wave 
length L/k, frequency Uk, and exponential growth pa- 
rameter \k (k £ {1,2, . . . ,N}). Inserting this into 
the linearized model equations gives consistency rela- 
tions for (Afe — iw/j), which have the form of so-called 
characteristic polynomials. These determine the possi- 
ble eigenvalues (A& — \u>k)- The homogeneous solution 
is stable with respect to small perturbations only if 
A fe < for all k £ {1,2,..., N}. Otherwise it is po- 
tentially unstable with respect to fluctuations £ Q (t). 
Specific analyses for different car-following models 
yield the instability conditions presented in Sec. III. A. 
For details see Herman et al. (1959) or Bando et al. 
(1995a). Typically one finds that traffic becomes un- 
stable if the reaction time At or the relaxation time r 
exceed a certain critical value. 

For macroscopic traffic models, we proceed simi- 
larly. Let us assume spatially homogeneous traffic of 
average density g and velocity V e (g), which is slightly 
disturbed according to Sp(x, t) = [p(x, t) — g] <g and 
5V(x,t) = [V(x,t) - V e (g)} < V e (g). The linearized 
density equation (62) reads, apart from fluctuations, 



dSV B 2 S P 
while the linearized velocity equation (63) is 
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(104) 



5p(x,t) = / dk pl ex v ( [kx + i x i(s) - ^(e)]*) . 
I J 

SV(x, t) = J2 I dkV l k cxp (ikx + [X l k (g) - kj l k (g)]t) , 
i J 

where p l k and V k l are amplitudes. The possible values 
of (X l k — iu) l k ) are given by the eigenvalues = [X l k — 
i(o>k — kV e (g))} of the matrix 



-D{g)k 2 -X{ 



ifc dP(g) 



—ikg 
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Stability requires the real values X l k {g) of the solutions 
of the corresponding characteristic polynomial to be 
negative for all wave numbers k. This leads to the 
stability condition 
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where equality yields the critical densities p c2 and 
p c3 determining the range of linear instability. For 
details see Kiihne and Rodiger (1991), and Hel- 
bing (1997a). As a consequence, traffic flow in the 
Lighthill-Whitham model with D(g) = and r — ► 
is marginally stable, while Eq. (53) and the re- 
lated Burgers equation (54) are always stable, if 
only D(g) > 0. Payne's model with D(g) = 
and dP{g)/dp = \dV e (g) / dp\/[2T(g)] becomes unsta- 
ble on the condition g\dV e (g)/dp\ > l/(2gr). This 
agrees with the instability condition (24) of the op- 
timal velocity model, if we consider that dv' e /dd = 
(dV e /dp)/(dd/dp) = -p 2 (dV e /dp). According to this, 
stability can be improved by decreasing the relax- 
ation time r (i.e. by higher acceleration). The insta- 
bility condition for the models by Kiihne, Kerner, and 
Konhauscr with D(g) = and dP(g)/dp = 0q is 
g\dV e (g)/dp\ > \/9q, given the road is long (Kiihne 
and Rodiger, 1991; Kerner and Konhauser, 1993). 
Hence, traffic flow becomes unstable if the velocity- 
density relation V e (p) falls too rapidly with increasing 
density, so that drivers cannot compensate for changes 
in the traffic situation sufficiently fast. 



2. (Auto-)Solitons, Korteweg-de-Vries equation, and 
Ginzburg-Landau equation 

Instead of using the original equations, in non-linear 
science it is common to investigate spatio-temporal 
patterns appearing in the vicinity of the instability 
threshold by means of approximate model equations 
based on scries expansions (Cross and Hohenberg, 
1993; Haken, 1997). In this way, the original equations 
can sometimes be replaced by the Korteweg-de-Vries 
equation 



The general solution of Eq. (104) is given by Fourier 
series of 8p(x,t) and 5V(x,t): 



where the tilde (~) indicates a scaling of the respec- 
tive variables, or by the Ginzburg-Landau equation 
(Cross and Hohenberg, 1993). Kiihne (1984b) was 
the first to apply this approach to the study of traf- 
fic dynamics (see also Sick, 1989). Moreover, Kurtze 
and Hong (1995) could derive a perturbed (modified) 
Korteweg-de-Vries equation and soliton-like solutions 
from the Kcrncr-Konhauser model for densities imme- 
diately above the critical density p CT . A similar result 
was obtained by Komatsu and Sasa (1995) for the mi- 
croscopic optimal velocity model (see also Whitham, 
1990; Igarashi et al., 1999a). 

A valid derivation of such model equations requires 
that the solutions of the original equations in the 
vicinity of the instability threshold are of small am- 
plitude and stable, otherwise the approximate equa- 
tions can be completely misleading (Cross and Hohen- 
berg, 1993). Therefore, the application of approximate 
model equations to traffic is restricted to the practi- 
cally irrelevant case with p c \ p C 2 ~ p C 3 ~ p C i, where 
instability barely exists. In the optimal velocity model, 
this condition would be fulfilled, if the inverse density 
and relaxation time were close to the critical point 
(d CI , l/r cr ) = (d c , 2 dv' e (d c ) / d d) , where d c is the turn- 
ing point of the velocity-distance relation (22) with 
d 2 v' e {d c )/dd 2 = 0. For that case, Nagatani (1998b, c; 
see also Muramatsu and Nagatani, 1999) derived the 
timc-dependent Ginzburg-Landau equation 
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with x' = a + 2tT[dv' e (d c )/dd] 2 , S' = (d- d c ), 
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and the thermodynamic potential 
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Moreover, Nagatani (1998b, c) calculated the uni- 
form and kink solutions, the coexistence curve de- 
fined by dip'/dS' = 0, and the spinodal line given by 
d 2 4>'/dS 12 = (see also Komatsu and Sasa, 1995). In 
his formalism, the inverse relaxation time 1/r plays 
the role of the temperature in a conventional phase 
transition, and the headways d corresponds to the or- 
der parameter. See also the study by Reiss et aL(1986) 
for an earlier application of thermodynamic ideas to 
traffic. 

Approximate equations for non-linear solutions of 
large amplitudes have been developed by Kerner et al. 
(1997), based on the theory of autosolitons (Kerner 
and Osipov, 1994). Moreover, the features of wide 
moving jams have been compared with autosolitons 
in physical systems (Kerner, 1995). They remind 



of properties found in some reaction- diffusion sys- 
tems (Mikhailov, 1991a, b; Meron, 1992; Kapral and 
Showalter, 1995; Woesler et al, 1996). 



3. Jam formation: Local breakdown and cluster effects, 
segregation, and self-organized criticality (SOC) 

Wave formation in traffic has some particular fea- 
tures. It turns out that the traffic equations are so 
highly non-linear that the linear approximation is only 
of very limited use. For example, the wave number 
k* associated with the largest growth rate A^, (g) at 
a given average density g does not at all determine 
the finally resulting wave length via 2ir/k*. More- 
over, the forming waves are typically not periodic. 
Even when starting with a small sinusoidal pertur- 
bation, the wave amplitude will not just grow lin- 
early. Instead, the extended perturbation will dras- 
tically change its shape, until it eventually becomes 
localized. Kerner and Konhauser (1994) call this the 
"local breakdown effect". The resulting perturbation 
has a characteristic form which can be approximated 
by the function 



p(x,t ) = g + Ap 



cosh 



2 / x — Xo 
w , 



w, _ 2 fx - x - (w + w_ ) 
— L cosh 

w_ \ w_ 
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with suitable parameters to, Xo, w, and w_ deter- 
mining its location and width (see Fig. 36a). When 
this shape has been reached, the perturbation grows 
more and more, playing the role of a nucleus for 
jam formation (Kerner and Konhauser, 1994). While 
small-amplitude perturbations flow with the traffic, 
the propagation speed becomes slower with increas- 
ing perturbation amplitude and eventually becomes 
negative. The final result is a wide traffic jam of char- 
acteristic form, which propagates upstream. This is, 
because vehicles are leaving the standing jam at the 
front, while new ones are joining the traffic jam at its 
end. The traffic jam is localized, which is known as "lo- 
cal cluster effect" (Kerner and Konhauser, 1994; Herr- 
mann and Kerner, 1998). Moreover, it is normally sur- 
rounded by free traffic flow. Hence, one could say that 
there is a phase separation (segregation) between free 
and congested traffic (Kerner and Konhauser, 1994; 
Kerner and Rehborn, 1996a). The result of this noise- 
induced ordering process (Helbing and Platkowski, 
2000) reminds of an equilibrium between two differ- 
ent phases (a freely moving, "gaseous" state and a 
jammed, "condensed" state). In other words, traffic 
jams absorb as much cars as necessary to have free 
traffic in the rest of the system. The resulting state 
of the system is not a partial congestion with all vehi- 
cles moving slow, as the velocity- density relation V e (p) 
would suggest. 

After the traffic jam has developed a stationary 
shape, it moves with constant velocity C < up- 



stream the road. Therefore, we have the relations 
p(x,t) = p(x - Ct,0) and Q(x,t) = Q(x - Ct,0). 
Inserting this into the continuity equation (45) gives 
—Cdp(x*,Q)/dx* + dQ(x*, 0)/dx* = with x* = 
(x — Ct) , which is solved by 

Q(z*,0) = Q + Cp(x*,0) = J(p(s„0)) (109) 

with a suitable constant Qo (Kerner and Konhauser, 
1994). Hence, the flow-density relation of a fully devel- 
oped traffic jam is a linear curve with negative slope 
C < 0, the so-called jam line J(p). While first-order 
macroscopic traffic models like the Lighthill-Whitham 
model have to assume this linear relation for congested 
traffic, in second-order models with a dynamic veloc- 
ity equation it is normally a self-organized relation 
which often differs significantly from the fundamental 
diagram Q e {p) (Kerner and Konhauser, 1994; Kerner 
and Rehborn, 1996a). Therefore, the fundamental dia- 
gram should be fitted only in the range of stable traffic 
flow, while congested traffic flow does not necessarily 
relate to the fundamental diagram. 

Finally, note that the above described mechanism 
of jam formation and the existence of the jam line 
J(p) seems to apply to all models with a deterministic 
instability mechanism, i.e. to models with a linearly 
unstable density regime (Bando et al, 1995b; Herr- 
mann and Kerner, 1998; Helbing and Schreckenberg, 
1999; Treiber et al, 1999, 2000). Many cellular au- 
tomata and other probabilistic traffic models have a 
different mechanism of jamming, but show a tendency 
of phase segregation between free and congested traf- 
fic as well (Nagel, 1994, 1996; see also Schreckenberg 
et al, 1995; Chowdhury et al, 1997a; Liibeck et al, 
1998; Roters et al, 1999). In a strict sense, phase seg- 
regation is found in slow-to-start models, while it is 
restricted in the Nagel-Schreckenberg model. Interest- 
ingly enough, the cruise control variant of the latter 
and a few other cellular automata display power-law 
behavior in the density variations (Choi and Lee, 1995; 
Csanyi and Kertesz, 1995; Yukawa and Kikuchi, 1995; 
Zhang and Hu, 1995; Nagatani, 1995b, 1998a). There- 
fore, Nagel and Herrmann (1993) as well as Nagel and 
Paczuski (1995) have interpreted jam formation as 
self-organized critical phenomenon (Bak et al, 1987), 
which is related with fractal self-similarity (Nagel and 
Rasmussen, 1994). 

4. Instability diagram, stop-and-go waves, metastability, and 
hysteresis 

To characterize the parameter-dependence of the 
possible states of a system resulting in the long run, 
phase diagrams are a very powerful method. They 
are of great importance in thermodynamics with var- 
ious applications in metallurgy, chemistry, etc. More- 
over, they allow us to compare very different kinds of 
systems like equilibrium and nonequilibrium ones, or 
microscopic and macroscopic ones, whose equivalence 



cannot simply be shown by transformation to normal 
forms (Kuramoto, 1989; Manneville, 1990; Eckmann 
et al, 1993). Defining universality classes by math- 
ematically equivalent phase diagrams, one can even 
compare so different systems as physical, chemical, 
biological and social ones, which is done in general 
systems theory (Buckley, 1967; von Bcrtalanffy, 1968; 
Rapoport, 1986). 

The phase diagram for the different traffic states 
on a homogeneous, circular one-lane road as a func- 
tion of the density g is usually called the instability 
diagram. Most traffic models predict a stable traf- 
fic flow at small vehicle densities and unstable traffic 
flow above a certain critical density p cr . At very high 
densities, many models with a deterministic instabil- 
ity mechanism behave stable again corresponding to 
creeping traffic, while the Nagel-Schreckenberg model 
and other cellular automata predict unstable traffic. 
Based on the calculation of limiting cases, a more de- 
tained picture has been suggested by Kuhne (1991a). 
He distinguishes supercritical and subcritical regimes, 
solitary waves, as well as shock fronts with either in- 
creasing or decreasing densities in downstream direc- 
tion. 

Simplifying matters a little, based on a numerical 
analysis of their macroscopic traffic model, Kerner and 
Konhauser found the following picture (see Fig. 36): 
Below some density p c \, any kind of disturbance even- 
tually disappears. Between the densities p c \ and p C 2, 
one wide traffic jam builds up, given a large enough 
perturbation. A series of traffic jams appears in a den- 
sity range between p C 2 and some density p c %. One so- 
called "anticluster" or "dipole layer" can be triggered, 
if the density g is between p c3 and p C 4, while any 
disturbance disappears in stable traffic above some 
density p c \. The critical densities p c k depend mainly 
on the choice of the pressure function, the relaxation 
time, and the velocity-density relation. 

The interpretation of these findings is as follows 
(Kerner and Konhauser, 1994): Traffic is linearly un- 
stable between the densities p C 2 and p c %. Therefore, 
the slightest inhomogeneity can cause a traffic jam. 
Such an inhomogeneity can even be the transition re- 
gion between a traffic jam and the surrounding traffic. 
Hence, an existing traffic jam will produce a new traf- 
fic jam (downstream of it), which generates another 
traffic jam, and so on, until there are not enough vehi- 
cles left to form another jam. Such an irregular series 
of traffic jams is called stop-and-go traffic. Kiihne's 
(1991a) interpretation of this irregularity was based 
on chaotic dynamics, motivated by an approximate 
mapping of his traffic model on the Lorenz equation 
(Kuhne and Beckschulte, 1993). However, numerical 
investigations indicate that the largest Lyapunov ex- 
ponent stays close to zero (Koch, 1996). 

Stop-and-go waves have been compared with waves 
in shallow water (Kuhne, 1984b; Hwang and Chang, 
1987), with the clogging of sand falling through a ver- 
tical pipe (Schick and Verveen, 1974; Baxter et al, 



1989; Lee, 1994; Poschel, 1994; Peng and Herrmann, 
1995) and of lead spheres falling through a fluid 
column (Horikawa et al, 1996; Nakahara and Isoda, 
1997), but the propagation direction is opposite. 

Kerner and Konhauser (1994) have shown that, in 
the density ranges [p c i,p C 2] and [p C 3,p C 4\, an existing 
traffic jam does not trigger any further jams, because 
traffic flow is not linearly unstable anymore. In these 
density regimes, they found metastable traffic, which 
is characterized by a critical amplitude Ap c (g) for the 
formation of traffic jams. This amplitude is zero for 
Q = Pc2 and g — p c3 , i.e. at the edge of unstable traf- 
fic, while it grows towards the edge of stable traffic 
and is expected to diverge at g = p c \ and g = p c ^ 
(see Fig. 36). Perturbations with subcritical ampli- 
tudes Ap < Ap c are eventually damped out (anal- 
ogous to the stable density ranges), while perturba- 
tions with supercritical amplitudes Ap > Ap c grow 
and form traffic jams (similar to the linearly unsta- 
ble density ranges). The situation in metastable traffic 
is, therefore, similar to supersaturated vapor (Kerner 
and Konhauser, 1994), where an overcritical nucleus 
is required for condensation ("nucleation effect"). 

FIG. 36. Schematic illustration of (a) the chosen local- 
ized perturbation and the density-dependent perturba- 
tion amplitudes Ap c required for jam formation in the 
metastable density regime, (b) the instability diagram of 
homogeneous (and slightly inhomogeneous) traffic flow, (c) 
the finally resulting traffic patterns as a function of the re- 
spective density regime. (Simplified diagram after Kerner 
et al, 1995a, 1996; see also Kerner and Konhauser, 1994). 



Interestingly, there are also traffic models without 
the phenomenon of metastability. Investigations by 
Kraufi (1998a, b) for a Gipps-like family of traffic mod- 
els (Gipps, 1981; see also McDonald et al, 1998) sug- 
gest the following general conclusions (see Fig. 37): 

(i) Models with a high ratio between maximum ac- 
celeration a and maximum deceleration b never 
show any structure formation, since traffic flow 
is always stable. 

(ii) Models with a relatively high maximum decel- 
eration b display a jamming transition, which 
is not hysteretic. As a consequence, there are 
no metastable high-flow states and the outflow 
from jams is maximal. 

(iii) A hysteretic jamming transition with metastable 
high-flow states and a characteristic, reduced 
outflow Q out from traffic jams (see Sec. IV. A. 5) 
is found for relatively small accelerations and 
medium decelerations. 



FIG. 37. (a) Schematic phase diagram of model classes 
as a function of the acceleration strength a and the de- 
celeration strength b. (b) Models belonging to phase (i) 
have no critical density p cr and display stable traffic in 
accordance with the fundamental diagram (thick dashed 
and solid lines) , while models belonging to phase (ii) show 
a non-hysteretic jamming transition with a resulting flow- 
density relation represented by solid lines, (c) A hysteretic 
phase transition with metastable high-flow states is found 
for models belonging to phase (iii). The model behavior is 
determined by the relative magnitude of the outflow Q ou t 
from congested traffic compared to the flow Q cr = Q e (p cr ) 
at which traffic becomes linearly unstable with respect to 
small perturbations. (After Kraufi, 1998a). 

The existence of high-flow states requires the flow 
Qout = Qe(pout) to be smaller than the flow Q CI = 
Qe(Pcr) at which traffic becomes linearly unstable, so 
that there is a density range of metastable traffic. By 
variation of model parameters, it is also possible to 
have the case Q cv < Q ou t, where traffic breaks down 
before it can reach high-flow states. This corresponds 
to case (ii) of the above classification by KrauB (1998a, 
b), where traffic flow does never exceed the jam line 
J(p) (see Fig. 37). The Nagel-Schreckenberg model be- 
longs to class (ii). Based on an analysis of the correla- 
tion function (Eisenblatter et al, 1998; Schadschnci- 
der, 1999; Neubert et al, 1999a) and other methods, 
it has been shown that the jamming transition is not 
of hysteretic nature, but continuous for p = 0. For 
< p < 1, there is a strong evidence for a crossover 
(Sasvari and Kertesz, 1997; Eisenblatter et al, 1998; 
Schadschneider, 1999; Neubert et al, 1999a; Chowd- 
hury et al, 2000a), although some people believe in 
critical behavior (Roters et al, 1999, 2000). 

Note that cellular automata with metastable high- 
flow states do exist, basically all slow-to-start models 
(Takayasu and Takayasu, 1993; Benjamin et al, 1996; 
Barlovic et al, 1998). However, it remains to be an- 
swered whether this metastability is of the same or 
a different type as described above, since in proba- 
bilistic models, the transition from high-flow states to 
jamming is due to the unfortunate adding up of fluc- 
tuations. One would have to check whether perturba- 
tions of the form (108) tend to fade away, when they 
are smaller than a certain critical amplitude Ap c , but 
always grow, when they are overcritical. Alternatively, 
one could determine the phase diagram of traffic states 
in the presence of bottlenecks (see Sec. IV.B.l). 

5. Self-organized, "natural" constants of traffic flow 

If the fundamental diagram Q e (p) is not linear in 
the congested regime, the jam line J(p) could, in prin- 
ciple, be a function of the average density g through- 
out the system. The same applies to the propagation 
velocity C of the jam, which is given by the slope of the 
jam line. However, in some traffic models, the jam line 
J(p) and the propagation velocity C are independent 



of the average density g, the initial conditions, and 
other factors (Kerner, 1999b), which is very surpris- 
ing. This was first shown by analytical investigations 
for wide jams in the limit of small viscosity (Kerner 
and Konhauscr, 1994; Kerner et al., 1997), based on 
the theory of autosolitons (Kerner and Osipov, 1994). 
Meanwhile, there is also empirical evidence for this 
(Kerner and Rchborn, 1996a). 

Note that the characteristic jam line J(p) defines 
some other constants as well: The intersection point 
(Pout,Qout) with the free branch of the fundamental 
diagram Q e (p) determines, because of the phase seg- 
regation between jammed and free traffic, the density 
Pout downstream of traffic jams. The quantity p ou t is 
a characteristic density between p c \ and p c2 , which 
depends on the choice of the traffic model and the 
specification of the model parameters. Moreover, the 
value Qout characterizes the outflow from traffic jams. 
Finally, the density pj am at which the jam line J(p) 
becomes zero, agrees with the density inside of stand- 
ing traffic jams. 

In summary, we have the following characteristic 
constants of wide traffic jams: the density p ou t down- 
stream of jams (if these are not travelling through 
"synchronized" congested traffic, see Kerner, 1998b), 
the outflow Qout = Qe(pout) from jams, the density 
Pj am inside of traffic jams, for which J(pj am ) = 0, and 
the propagation velocity C = <2 e (/9out)/(/W - Pjam) 
of traffic jams. Note that these "natural constants" 
of traffic are not at all a trivial consequence of the 
construction of the corresponding models. They are 
rather a result of self-organization within traffic, i.e. 
dynamical invariants. These are related with an at- 
tractor of the system dynamics, determining the form 
of the wave fronts. In non-linear systems, however, the 
size and shape of an attractor usually changes with the 
control parameter, which is the density, here. The at- 
tractor in the congested traffic regime can be only in- 
dependent of the density due to the hysteretic nature 
of the transition and the associated phase segregation 
between free and jammed traffic. 

Not all traffic models display this independence. 
For example, although the microscopic optimal ve- 
locity model produces a similar dynamics as the KK 
model (58) (Berg et al., 2000), for certain velocity- 
distance functions v' e {d) its propagation velocity of 
traffic jams depends on the average density g (Herr- 
mann and Kerner, 1998). This happens, if the result- 
ing jam density pj am varies with the vehicle velocity 
when approaching a jam and, therefore, with the aver- 
age density g and the initial conditions. The density- 
dependence can, however, be suppressed by suitable 
discretizations of the optimal velocity model (Hclbing 
and Schreckenberg, 1999) or by models in which the 
approaching process depends on the relative veloc- 
ity, for example the GKT model or the IDM (Trcibcr 
et al, 1999, 2000). It is also suppressed by particu- 
lar velocity-distance functions such as the simplified 
relation v e {d) = v Q'(d — d ), where O'(z) denotes 



the Heaviside function, which is 1 for z > and oth- 
erwise 0. For this choice, the optimal velocity model 
can be exactly solved. Sugiyama and Yamada (1997) 
find that the dynamics of fully developed jams can 
be represented by a hysteresis loop in the speed-over- 
distancc plane rather than the above Heaviside func- 
tion (see Fig. 38). Moreover, the characteristic con- 
stants of traffic flow are related through the equation 
1/Pout = 1/Pjam + v T' , where V = (t a - i a _i) is 
the characteristic time interval between the moments 
t a , when successive cars a have a distance do to the 
respective leading car (a — 1) and start accelerating 
in order to leave the downstream jam front. The den- 
sities p out and pj am of free and jammed traffic, respec- 
tively, follow from l/p ou t = d + vqT'/2. Moreover, 
the resolution velocity of jam fronts can be deter- 
mined as C = [x a (t a ) - x a -i(t a -i)]/[t a - t a -i] = 
— 1/pjamT' = — l/pj am T. Finally, the time interval 
T" = T between accelerating vehicles is given by 
T'/t = 2(1 -e~ T 'l T ). 

FIG. 38. Hysteresis loop of fully developed traffic jams 
(solid lines with arrows) , emerging in the optimal velocity 
model (23) with the equilibrium speed-distance relation 
v e (d) — voQ'(d-do), see the dashed line. (After Sugiyama 
and Yamada, 1997.) 

Summarizing the above results, one finds character- 
istic traffic constants, if the jam density pj am develop- 
ing at the upstream end of traffic jams is independent 
of the surrounding traffic. In contrast to pj am , the time 
interval T' between accelerating vehicles is determined 
by the dynamics at the downstream jam front. It re- 
sults from the fact that the initial conditions of accel- 
erating vehicles at jam fronts are more or less identi- 
cal: The cars have equal headways rfj am s=s l/pj am to 
the respective leading vehicle and start with the same 
velocity v w 0. The calculation of the characteristic 
constants is nevertheless a difficult task. However, for 
different models, analytical results have been obtained 
by Kerner et al. (1997), Helbing and Schreckenberg 
(1999) as well as Gerwinski and Krug (1999). 

6. Kerner's hypotheses 

It appears that the above theoretical results cover 
many of the observations summarized in Sec. II. E. 
However, Kerner (1998a, b, 1999a, b, 2000a, b, c) was 
surprised by the following aspects and has, therefore, 
questioned previously developed traffic models: 

(i) Synchronized flows of kind (i) and (ii) (see 
Sec. II. E. 2) can exist for a long time (Kerner and Rch- 
born, 1996b), i.e., they can be stable, at least with 
respect to fluctuations of small amplitude (see also 
the discussion of stable congested flow by Westland, 
1998). 

(ii) Wide jams and stop-and-go traffic are rarely 
formed spontaneously (Kerner and Rchborn, 1997; 
Kerner, 1999c, 2000c), but occur reproducible at the 



same freeway sections during rush hours. Instead of a 
transition from free traffic to stop-and-go traffic, one 
normally observes two successive transitions: a hys- 
teretic one from free flow to synchronized flow (Kerner 
and Rehborn, 1997) and another one from synchro- 
nized flow to stop-and-go waves (Kerner, 1998a; see 
Sec. II.E.3). 

In conclusion, the transition to synchronized flow 
appears more frequently than the transition to stop- 
and-go traffic. For an explanation, Kerner (1998a, b, 
1999a, c, 2000a, b) has, therefore, developed some hy- 
potheses in the spirit of his three-phase traffic theory 
(see Sees. II. B and II. E), which I try to sketch in the 
following: 

(i) The whole multitude of hypothetical homogeneous 
and stationary (steady) states of synchronized flow is 
related to a two-dimensional region in the flow-density 
plane. This region is the same for a multi-lane road 
and for a one-lane road. 

(ii) All hypothetical homogeneous and stationary 
(steady) states of free flow and synchronized flow are 
stable with respect to infinitesimal perturbations. 

(iii) There are two qualitatively different kinds of nu- 
cleation effects and the related first-order phase tran- 
sitions in the hypothetical homogeneous states of traf- 
fic flow: 

- the nucleation effect responsible for jam formation, 
the nucleation effect responsible for the phase tran- 
sition from free to synchronized flow. 

(iv) At each given density in homogeneous states of 
free flow, the critical amplitude of a local perturbation 
of the density and/or velocity, which is needed for a 
phase transition from free flow to synchronized flow, 
is considerably lower than the critical amplitude of 
a local perturbation needed for jam formation in free 
flow. In other words, the probability of the phase tran- 
sition from free flow to synchronized flow is consider- 
ably higher than the probability of jam emergence in 
free flow. 

(v) The jam line J(p) (see Sec. IV. A. 3) determines the 
threshold of jam formation. All (i.e. an infinite num- 
ber of) states of traffic flow related to the line J are 
threshold states with respect to jam formation: Traffic 
states below the jam line J(p) would always be stable, 
while states above it would be metastable. At the same 
distance above the line J(p), the critical amplitude of 
local perturbations triggering jam formation is higher 
in free flow than in synchronized flow. Therefore, a 
transition to a wide moving jam occurs much more 
frequently from synchronized than from free flow. 

(vi) The transition from free to synchronized flow is 
due to an avalanche self-decrease in the mean proba- 
bility of overtaking, since it is associated with a syn- 
chronization among lanes. 

(vii) The hysteretic transition to wide jams is related 
to an avalanche-like growth of a critical perturbation 
of the density, average velocity, or traffic flow, see the 
nucleation effect described in Sec. IV. A. 4. 

Simulation models and further empirical evidence 



for these hypotheses are still to be found. 

B. Transition to congested traffic at bottlenecks and 
ramps 

Effects of changes in the number of lanes and ramps 
have been simulated since the early days of applied 
traffic simulations (see, for example, Munjal et ai, 
1971; Phillips, 1977; Cremer, 1979; Makigami et al, 
1983), but a systematic study of the observed phe- 
nomena, their properties, mechanisms, and precondi- 
tions has been missing. Later, there have been various 
theoretically motivated studies by physicists as well. 
These assume, for example, local "defects" , "impu- 
rities", or bottlenecks with reduced "permeability", 
which are mostly associated with a locally decreased 
velocity. Typical results for one single defect are 

(i) a drop of the freeway capacity to a value Qbot at 
the impurity, 

(ii) a segregation into one area of free and and one 
area of congested traffic, 

(iii) a localization of the downstream front of con- 
gested traffic at the bottleneck, 

(iv) a shock-like propagation of the upstream conges- 
tion front in agreement with the Lighthill-Whitham 
theory, particularly formula (51) with Q + = Qbot, 

(v) a stationary length of the congested area on cir- 
cular roads, since the inflow Q- is determined by the 
outflow Q out from the bottleneck. 

Investigations of this kind have been carried out 
for particle hopping models like the TASEP with 
random sequential update (Janowky and Lebowitz, 
1992, 1994) and for traffic models with parallel up- 
date (Chung and Hui, 1994; Csahok and Vicsek, 1994; 
Emmerich and Rank, 1995). One is certainly tempted 
to compare the congested traffic appearing in these 
models with the observed "synchronized" congested 
flow discussed before, but several elements seem to be 
missing: 

(i) the hysteretic nature of the transition, character- 
ized by the nucleation effect, i.e. the requirement of 
an overcritical perturbation, 

(ii) the typical dynamics leading to the formation of 
"synchronized" congested flow (see Sec. II. E. 2), and 

(iii) the wide scattering of flow-density data. 

In other simulations with inhomogencitics, re- 
searchers have recognized wave-like forms of congested 
traffic behind bottlenecks (Csahok and Vicsek, 1994; 
Hilliges, 1995; Klar et al, 1996; Nagatani, 1997c; Lee 
et al, 1998). One can even find period doubling sce- 
narios upstream of stationary bottlenecks (Nagatani, 
1997c) or due to time-dependent local perturbations 
(Koch, 1996; Nicolas et al, 1994, 1996), which are 
analytically understandable. Moreover, Kerner et al. 
(1995a) have simulated slightly inhomogeneous traf- 
fic on a circular freeway with on- and off-ramps, for 
which in- and outflows were chosen identically. They 
found the interesting phenomenon of a localized sta- 



tionary cluster along the on-ramp, which triggered a 
traffic jam, when the on-ramp flow was reduced. How- 
ever, these and other simulations with the KK model 
(58) were not fully consistent with the above outlined 
phenomena, so that Kerner developed the hypothe- 
ses sketched in Sec. IV. A. 6. I believe, the reason for 
the discrepancy between observations and simulations 
may have been the assumptions of periodic boundary 
conditions and uniform traffic, i.e. identical driver- 
vehicle units (see below and Sec. IV. C. 2). This is, how- 
ever, still a controversial topic which certainly requires 
more detailed investigations in the future. 

In 1998, Lee et al. have tried to reproduce the hys- 
teretic phase transition to "synchronized" congested 
flow by simulation of a circular road with the KK 
model including on- and off-ramps, but with different 
parameters and a modified velocity-density function 
V e (p). By a temporary peak in the on-ramp flow, they 
managed to trigger a form of stop-and-go traffic that 
was propagating upstream, but its downstream front 
was pinned at the location of the ramp. They called it 
the "recurring hump" state (RH) and compared it to 
autocatalytic oscillators. Free traffic would correspond 
to a point attractor and the oscillating traffic state to 
a stable limit cycle. In terms of non-linear dynamics, 
the transition corresponds to a Hopf bifurcation, but 
a subcritical one, since the critical ramp flow depends 
on the size of the perturbation. (After Hclbing and 
Treiber, 1998b). 

Lee et al. have pointed out that free traffic survives 
a pulse-type perturbation of finite amplitude, if the 
ramp flow is below a certain critical value. However, 
once a RH state has formed, it is self-maintained until 
the ramp flow falls below another critical value which 
is smaller than the one for the transition from FT to 
the RH state. This proves the hysteretic nature of the 
transition. Moreover, Lee et al. have shown a gradual 
spatial transition from the RH state to free flow down- 
stream of the ramp. They have also demonstrated a 
synchronization among neighboring freeway lanes as 
a result of lane changes. Therefore, they have sug- 
gested that their model can describe the empirically 
observed first-order phase transition to synchronized 
flow, where the two-dimensional scattering of "syn- 
chronized" congested traffic is understood as a result 
of the fact that the amplitude of the oscillating traffic 
state depends on the ramp flow. The similarity with 
empirical data, however, was rather loose. Also the 
central question regarding an explanation of homoge- 
nous and stationary congested traffic (see Sees. IV. A. 6 
and II. E. 2) remained unsolved. 

Independently of Lee et al., Helbing and Treiber 
(1998a) had submitted another study two weeks 
later. Based on the non-local, gas-kinetic-based traf- 
fic model, they simulated a freeway section with open 
boundary conditions and one on-ramp, since the tran- 
sition to "synchronized" congested flow had been 
mostly observed close to ramps. Ramps were modelled 
by means of a source term in the continuity equation: 
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Herein, Q vmp (t) represents the ramp flow per freeway 
lane entering the /-lane freeway over a ramp section 
of effective length L. The simulation scenario assumed 
that the sum 

Qdown = Qup + Qrmp(t)/I (111) 

of the flow Q up upstream of the ramp and the ramp 
flow Q rmp (t)/I per lane had reached the regime of 
metastable traffic flow with Q e (p c i) < Qdown < 
Qe{pc2) downstream of the ramp. Hence, without any 
perturbation, there was free traffic flow. However, a 
short overcritical peak in the ramp flow Q rmp (i) trig- 
gered a growing perturbation, which travelled down- 
stream in the beginning, but changed its propagation 
speed and direction as it grew larger. Consequently, it 
returned to the ramp ("boomerang effect") and initi- 
ated a continuously growing region of extended con- 
gested traffic when it arrived there (see Fig. 39). The 
reason for this is basically the dynamical reduction 
of the freeway capacity to the outflow Q out from con- 
gested traffic (cf. Sees. II. E and IV. A. 5). This allows us 
to understand what Daganzo et al. (1999) call the "ac- 
tivation" of a bottleneck. Note that the outflow Q ou t 
depends on the ramp length L (Hclbing and Treiber, 
1998a; Treiber et al, 2000), but in some traffic mod- 
els and in reality it is also a function of the ramp 
flow Q rmp , the average speed Vbot in the congested 
area, and possibly other variables (such as the delay 
time by congestion, since drivers may change their 
behavior). We expect the relations dQ out /dL > 0, 
<9<3out/<9<3rm P < 0, dQ out /dV hot > 0, and 

lim lim lim Q out (Q lmp /I,V hot ,L) = Q out . 

(112) 

Note that both, Q out and Q ou t are also functions of 
the model parameters Vq, T, etc. (and of their spatial 
changes along the road). 

FIG. 39. Spatio-temporal evolution of the density after a 
small peak in the inflow from the on-ramp. The on-ramp 
merges with the main road at x = km. Traffic flows from 
left to right. In (a), the parabolically shaped region of high 
density corresponds to the resulting "synchronized" con- 
gested flow. Plot (b) illustrates in more detail that the 
congestion starts to build up downstream of the bottle- 
neck (cf. Sec. II. E. 2). The time-dependent inflows Q up at 
the upstream boundary and Qrmp// at the on-ramp are 
displayed in (c). (d) The simulated traffic flows at two cross 
sections show the temporary drop below the finally result- 
ing bottleneck flow ( — ) and discharge flow ( — ) observed 
in real traffic, cf. Fig. 3(b). (After Helbing and Treiber, 
1998a; Helbing et at, 2001a, b.) 

The (congested) bottleneck flow 

Qbot = Qout — Qrmp/I = Qout _ AQ (H3) 



immediately upstream of the ramp is given by the out- 
flow Q ou t minus the ramp flow Q rmp /I per lane defin- 
ing the bottleneck strength AQ. Let /?bot > Pc2 be the 
associated density according to the congested branch 
of the fundamental diagram, i.e. Qbot = Qe(Pbot)- 
Then, the congested traffic flow upstream of the bot- 
tleneck is homogeneous, if pbot falls in the stable 
or metastable range, otherwise one observes differ- 
ent forms of congestion (see Sec. IV.B.l). Hence, the 
above hysteretic transition to homogeneous congested 
traffic (HCT) can explain the occurence of homoge- 
neous and stationary (steady) congested traffic states. 
The scattering of "synchronized" congested flow can 
be accounted for by a mixture of different vehicles 
types such as cars and trucks (Treiber and Helbing, 
1999a, see Sec. IV. C. 2). However, apart from a het- 
erogeneity in the time headways (Banks, 1999), there 
are also other plausible proposals for this scattering, 
which are more difficult to verify or falsify empirically: 
complex dynamics with forward and backward prop- 
agating shock waves (Kcrncr and Rehborn, 1996a; 
Kcrner, 1998b, 2000a, b), changes in the behavior of 
"frustrated" drivers (KrauB, 1998a), anticipation ef- 
fects (Wagner, 1998a; Knospe et al, 2000a, b), non- 
unique equilibrium solutions (Nelson and Sopasakis, 
1998; Nelson, 2000), or multiple metastable oscillating 
states (Tomer et al, 2000). Non-unique solutions are 
also expected for the non-integer car-following model 
(19), where the vehicle acceleration dv/dt is always 
zero when the relative velocity Av becomes zero, no 
matter how small the distance is. However, most of 
these states are not stable with respect to fluctuations. 

Let us now explain why the downstream front of 
extended congested traffic is located at the bottleneck 
and stationary: If it would be located downstream of 
the ramp for some reason, its downstream front would 
travel upstream as in traffic jams. So the question is 
why it does not continue to move upstream, when the 
downstream front has reached the bottleneck. When 
it would pass the on-ramp, it would continue to emit 
the flow Q ou t, leading to an increased flow (Qout + 
Qrmp/I) > (Qup + Qrmp/I) at the ramp, which would 
again queue up. The spatial extension of the congested 
area is only reduced, when the traffic flow Qdown falls 
below the outflow Q out of congested traffic. This is 
usually the case after the rush hour. 

Finally, I should mention that congested traffic can 
be also triggered by a perturbation in the traffic flow 
Q U p(i) entering the upstream freeway section. It can 
be even caused by a supercritical "negative" pertur- 
bation with a reduction of the vehicle density and flow 
as, for example, at an off-ramp (Helbing, 2001). This 
is, because vehicles will accelerate into the region of 
reduced density, followed by a successive deceleration. 
Therefore, a negative perturbation of the density is 
rapidly transformed into the characteristic perturba- 
tion which triggers congested traffic (see Sec. IV. A. 3). 



1. Phase diagram of traffic states at inhomogeneities 

Theoretical investigations of traffic dynamics have, 
for a long time, mainly been carried out for circular 
roads. The corresponding results were, in fact, mis- 
leading, since it is usually not possible to enter the 
linear unstable region g > p c2 of congested traffic on 
a homogeneous road. High densities were produced 
by means of the initial condition, but they cannot 
be reached by any natural inflow into a homogeneous 
freeway with open boundaries, no matter how high the 
inflow is. Consequently, spontaneously forming traffic 
jams are rather unlikely as well, apart from tempo- 
rary ones caused by sufficiently large perturbations. 
However, the situation changes drastically for inho- 
mogeneous roads with ramps or other bottlenecks. 

Helbing et al. (1999) have systematically explored 
the resulting traffic states at inhomogeneities trig- 
gered by a fully developed perturbation (jam) as a 
function of the ramp flow Q rmp // per freeway lane and 
the upstream traffic flow Q up . In this way, they found 
a variety of congested traffic states and their relation 
to each other (see Fig. 40). Free traffic (FT) was, of 
course, observed, if the total downstream traffic flow 
Qdown = (Qup + Qrmp/I) was stable with respect to 
large perturbations, i.e. if Qdown < Qe(Pci)- As ex- 
pected, extended congested traffic states were found 
when the traffic flow Qdown exceeded the maximum 
flow Q m ax = maxp Q e (p), i.e. the theoretical freeway 
capacity for homogeneous traffic. However, this classi- 
cal kind and many other kinds of congestion occurred 
even for smaller traffic flows Qdown < Qmax, triggered 
by perturbations. Consequently, all these kinds of con- 
gestions could be avoided by suppressing perturba- 
tions with suitable technical measures (see Sec. IV. G). 

FIG. 40. Numerically determined phase diagram of the 
traffic states forming in the vicinity of an on-ramp as a 
function of the inflows Q up and Q rm p/ 1 per freeway lane on 
the main road and the on-ramp. Displayed are the phases 
of homogeneous congested traffic (HCT), oscillatory con- 
gested traffic (OCT), triggered stop-and-go traffic (TSG), 
moving localized clusters (MLC), pinned localized clus- 
ters (PLC), and free traffic (FT). The states are triggered 
by a fully developed localized cluster travelling upstream 
and passing the ramp (cf. Fig. 41). Dashed lines indicate 
the theoretical phase boundaries. The lowest diagonal line 
Qdown = Qe(pd) separates free traffic from wide jams, 
while the central diagonal line Qdown = Qout ~ Qe{pc2) 
separates localized from extended forms of congested traf- 
fic. The localized cluster states between these two di- 
agonals require metastable downstream traffic. Finally, 
the highest (solid) diagonal line represents the condition 
Qdown = Qmax limiting the region, in which a breakdown 
of traffic flow is caused by exceeding the theoretically pos- 
sible freeway capacity Q m ax (upper right corner). All con- 
gested traffic states below this line are caused by pertur- 
bations and, therefore, avoidable by technical control mea- 
sures. (After Helbing et al, 1999, 2001a.) 



The "classical" form of extended congested traffic is 



HCT (see Fig. 41a), which is formed if the density /?bot 
in the congested region falls into the (meta-)stable 
regime. Decreasing the ramp flow Q rmp will increase 
the flow Qbot = (Qout - Qrmp/I) in the congested 
region and decrease the related density pbot- When 
it falls into the unstable regime, i.e. when we have 
Pbot < Pc3 or Qbot > Qe(pc3), we can expect oscillat- 
ing congested traffic (OCT, see Fig. 41b). However, in 
deterministic models and without any additional per- 
turbations, homogeneous congested traffic extends up 
to Qbot > Qe(Pcv) with a smaller density p cv < p c3 . 
The reason is that, above a density of p cv , traffic flow is 
convectively stable (Hclbing et ai, 1999a), i.e. pertur- 
bations are convected away (Manneville, 1990; Cross 
and Hohenberg, 1993). While the transition from free 
traffic to congested traffic is hysteretic (see Sec. IV. B), 
the transition from the HCT to the OCT state is con- 
tinuous. A suitable order parameter to characterize 
this transition is the oscillation amplitude. 

FIG. 41. Spatio-temporal dynamics of typical represen- 
tatives of congested states which are triggered at an on- 
ramp by a fully developed perturbation travelling up- 
stream. The middle of the on-ramp is located at x = 8.0 
km. The displayed states are (a) homogeneous congested 
traffic (HCT), (b) oscillatory congested traffic (OCT), (c) 
triggered stop-and-go traffic (TSG), (d) a pinned localized 
cluster (PLC). (From Helbing et ai, 2001a, b; after Hel- 
bing et ai, 1999a.) 

Let us assume we further reduce the ramp flow, so 
that congestion is expected to become less serious. 
Then, at some characteristic ramp flow Q® mp , the os- 
cillation amplitude becomes so large that the time- 
dependent flow-density curve reaches the free branch 
of the fundamental diagram Q e (p). For Q rmp < Q® mp 
we, therefore, speak of a triggered stop-and-go state 
(TSG). A suitable order parameter would be the av- 
erage distance between locations of maximum con- 
gestion. Drivers far upstream of a bottleneck will be 
very much puzzled by the alternation of jams and 
free traffic, since these jams seem to appear with- 
out any plausible reason, analogous to phantom traffic 
jams. The triggering mechanism of repeated jam for- 
mation is as follows: When the upstream travelling 
perturbation passes the inhomogeneity, it triggers an- 
other small perturbation, which travels downstream. 
As described in Sees. II. E and IV. A. 3, this pertuba- 
tion changes its propagation speed and direction while 
growing. As a consequence, the emerging jam returns 
like a boomerang and triggers another small pertur- 
bation while passing the inhomogeneity. This process 
repeats time and again, thereby causing the triggered 
stop-and-go state. 

If the traffic flow Qdown downstream of the inho- 
mogeneity is not linearly unstable, the triggered small 
perturbation cannot grow, and we have only a single 
localized cluster (LC), which normally passes the in- 
homogeneity and moves upstream. The condition for 
this moving localized cluster (MLC, see Fig. 41c) or 



wide moving jam is obviously Qdown < Qe{p C 2) ■ How- 
ever, if Q up < Q e (pci), i-e. when the flow upstream of 
the bottleneck is stable, a congested state cannot sur- 
vive upstream of the inhomogeneity. In that case, we 
find a standing or pinned localized cluster (SLC/PLC) 
at the location of the bottleneck (see Fig. 41d). Os- 
cillatory forms of pinned localized clusters (OPLC) 
are possible as well (Lee et al, 1999; Treiber et ai, 
2000). The transition between localized clusters and 
extended forms of congested traffic is given by the 
condition Qdown > Qout, requiring that the traffic flow 
Qdown exceeds the outflow Qout from congested traffic. 

Finally, I would like to repeat that the transition 
between free and congested traffic is of first order (i.e. 
hysteretic), while the transitions between extended 
forms of congested traffic (HCT, OCT, and TSG) 
seem to be continuous. As a consequence of the for- 
mer, one can have direct transitions from free traffic to 
homogeneous or oscillating congested traffic (HCT or 
OCT) without crossing localized cluster (LC) states 
or triggered stop-and-go waves (TSG). Normally, the 
traffic flows Q U p and Qrmp during the rush hour in- 
crease so drastically that a phase point belonging to 
the area of extended congested traffic (CT) is reached 
before free traffic flow breaks down. The theory pre- 
dicts this to happen at weekdays characterized by high 
traffic flows, while at the same freeway section LC or 
TSG states should rather appear at weekdays with 
less pronounced rush hours. 

It should be stressed that the above analytical re- 
lations for the transitions between free traffic and the 
five different forms of congested traffic (PLC, MLC, 
TSG, OCT, and HCT) are in good agreement with nu- 
merical results (see Fig. 40). Moreover, according to 
the quantities appearing in these analytical relations, 
the transitions are basically a consequence of the in- 
stability diagram, combined with the fundamental di- 
agram and the jam line. Therefore, we expect that 
the above phase diagram is universal for all micro- 
scopic and macroscopic traffic models having the same 
instability diagram. This prediction has been mean- 
while confirmed for several models. Lee et al. (1999) 
have found analogous transitions for the KK model 
(58), but recognized a region of tristability, where it 
is a matter of history whether one ends up with free 
traffic, a pinned localized cluster, or oscillating con- 
gested traffic. This tristability is probably a result 
of a dynamical variation of Qout, as pointed out by 
Treiber et al. (2000), who have also reproduced the 
phase diagram for the IDM (see Sec. III. A. 3). More- 
over, they have presented empirical data confirming 
the existence of all five congested traffic states. These 
have been successfully reproduced in simulation sce- 
narios with the IDM, using the measured boundary 
conditions. Only one single parameter, the safe time 
clearance T, has been varied to calibrate both, the ca- 
pacity of the different freeway sections and their bot- 
tlenecks. It has also been demonstrated that different 
kinds of bottlenecks have qualitatively the same ef- 



fects. This requires a generalized definition of bottle- 
neck strengths AQ according to 

AQ = Qrmp/7 + Qout - Qout . ("4) 

where Qout denotes the outflow from congested traffic 
at the freeway section with the largest capacity and 
Q' out the outflow at the respective bottleneck. Spatial 
changes in the number I(x) of lanes can be reflected 
by virtual ramp flows (Shvetsov and Hclbing, 1999): 

Qrmp _ pV dl(x) 

IT- ~i dx~ m (Ub) 

This follows from the continuity equation d(Ip) jdt + 
d(IpV)/dx = 0. 

Recently, Lee et al. (2000) have presented an em- 
pirical phase diagram, but the corresponding freeway 
stretch contains several capacity changes, which com- 
plicates matters. Such freeway sections show often 
congestion with a stationary (pinned) upstream end. 
It can also happen that a jam travels through a short 
section of stable traffic, because there is a finite pen- 
etration depth (Treiber et al, 2000). Finally, I would 
like to mention a number of other simulation stud- 
ies with open boundaries, which have recently been 
carried out (Bengrinc et al., 1999; Cheybani et al., 
2000a, b; Mitarai and Nakanishi, 1999, 2000; Popkov 
and Schiitz, 2000; Popkov et al, 2000). 

Summarizing the results, a comparison between nu- 
merical simulations and empirical data suggest the fol- 
lowing identifications: Homogeneous congested traffic 
(HCT) seems to be the same as "sychronized" traffic 
flow (ST) of type (i) (see Sec. II.E.2), while oscillat- 
ing congested traffic (OCT) seems to be related to 
synchronized flow of type (iii) . Homogcncous-in-speed 
states, which are also called synchronized flow of type 
(ii) , are found where congested traffic relaxes towards 
free traffic (see Figs. 43b and 5; Tilch and Helbing, 
2000, Fig. 3, right). Moreover, the term "recurring 
hump state" (RH state, see Lee et al, 1998, 1999) 
summarizes all oscillating forms of congested traffic 
(OCT, TSG, and OPLC), while "non-homogeneous 
congested traffic" was used for all forms of unstable 
traffic (Tomer et al, 2000). 

Let us now ask, what the above theory predicts 
for models with another instability diagram. If the 
model is stable in the whole density range such as the 
Burgers equation or the TASEP, we just expect two 
states: free traffic (FT) and homogeneous congested 
traffic (HCT), which is in agreement with numerical 
(Janowsky and Lebowitz, 1992, 1994) and exact an- 
alytical results (Schiitz, 1993). If the model displays 
stable traffic at small densities and unstable traffic at 
high densities, one expects free traffic (FT) and os- 
cillating traffic (OCT or TSG). This was, for exam- 
ple, observed by Emmerich and Rank (1995) for the 
Nagel-Schreckenberg model (see also Diedrich et al, 
2000; Cheybani et al, 2000a). For the case of stable 
traffic at small and high densities, but unstable traffic 



at medium densities, one should find free traffic (FT), 
oscillating traffic states (OCT or TSG), and homoge- 
neous congested traffic (HCT). This could possibly be 
tested with the Weidlich-Hilligcs model. The existence 
of moving and pinned localized clusters (MLC/PLC) 
would normally require a metastable state between 
free traffic and linearly unstable traffic. Finally, as- 
sume there were models with stable traffic at low den- 
sities and metastable traffic above a certain critical 
density p c \, but no range of linearly unstable traffic, 
as suggested by Kerner (see Sec. IV. A. 6). Then, os- 
cillatory congested traffic (OCT) and triggered stop- 
and-go waves (TSG) should not exist. Identifying lo- 
calized cluster states (MLC/PLC) with wide moving 
jams, but homogeneous congested traffic (HCT) with 
synchronized flow would give a similar picture to the 
three-phase traffic theory proposed by Kerner and 
Rchborn (see Sec. II. B): The possible traffic states 
would be free traffic, wide jams, and synchronized 
flow. Transitions to synchronized flow would occur, 
when the traffic flow Qdown would exceed the maxi- 
mum capacity <5 max - This would tend to happen at 
bottlenecks, where the metastable density regime is 
entered first. The transition would cause a queueing of 
vehicles and a capacity drop due to the increased time 
headways of accelerating vehicles at downstream con- 
gestion fronts. In contrast, transitions to wide mov- 
ing jams would always have to be triggered by su- 
percritical perturbations. This would most frequently 
happen in congested areas upstream of bottlenecks 
("pinch effect", see Sec. II. E. 3), as metastable traf- 
fic is rather unlikely on homogeneous stretches of a 
freeway. Supplementary, the scattering of flow-density 
data in "synchronized" congested traffic could be ex- 
plained by some source of fluctuations, e.g. heteroge- 
neous driver- vehicle behavior (see Sec. IV. C. 2). 

The above considerations demonstrate how power- 
ful the method of the phase diagram is. Based on a few 
characteristic quantities (the critical flows), one can 
predict the dynamical states of the respective model 
as well as the phase boundaries. Moreover, one can 
classify traffic models into a few universality classes. 
It does not matter whether the model is of microscopic 
or macroscopic, of deterministic or of probabilistic na- 
ture. 

2. Coexistence of states and "pinch effect" 

In simulations with the IDM, Treiber and Hclbing 
(1999b) have simulated a scenario which reminds of 
the pinch effect reported by Kerner (1998a) and sim- 
ilar observations by Koshi et al. (1983). This sce- 
nario is characterized by the spatial coexistence of 
different forms of congested traffic, which can occur 
upstream of bottlenecks even on rampless road sec- 
tions, as the transitions between the concerned con- 
gested states are continuous (see Sec. IV.B.l). Start- 
ing with homogeneous congested traffic (HCT) at a 



bottleneck, Treiber and Hclbing found oscillating con- 
gested traffic (OCT) upstream of it, but stop-and- 
go traffic even further upstream (see also Helbing 
et al, 2001b). For comparison, upstream of a bottle- 
neck Kcrncr has observed transitions from synchro- 
nized flow to a "pinch region" characterized by nar- 
row clusters, which eventually merged to form wide 
moving jams (see Fig. 42 and Sec. II. E. 3). Personally, 
I valuate this phenomenon as strong evidence for the 
existence of linearly unstable traffic. 

FIG. 42. Illustration of the spatial coexistence of HCT, 
OCT and TSG states in a simulation with the IDM for 
an inhomogeneous freeway without ramps. Traffic flows 
in positive ^-direction. The spatio-temporal density plot 
shows the breakdown to homogeneous congested traffic 
near the inhomogeneity and stop-and-go waves emanating 
from this region. The inhomogeneity corresponds to an in- 
creased safe time clearance T' > T between x — km and 
x — 0.3 km, reflecting more careful driving. Downstream 
of the inhomogeneity, vehicles accelerate into free traffic. 
(From Treiber and Helbing, 1999b; Helbing et al, 2001b.) 

According to the model by Treiber and Helbing 
(1999b), the conditions for this spatial coexistence of 
congested traffic states are the following. The density 
Pbot in the congested region immediately upstream of 
the bottleneck should be in the linearly unstable, but 
convectively stable range [p cv ,Pc3], where perturba- 
tions are convected away in upstream direction (Man- 
neville, 1990; Cross and Hohenbcrg, 1993). In this 
case, traffic flow will appear stationary and homoge- 
neous close to the bottleneck, but small perturbations 
will grow as they propagate upstream in the congested 
region starting at the bottleneck. If the perturbations 
propagate faster than the congested region expands, 
they will reach the area of free traffic upstream of 
the bottleneck. During rush hours, it is quite likely 
that this free flow is in the metastable range between 
Pd and p C 2- Consequently, sufficiently large perturba- 
tions can trigger the formation of jams, which con- 
tinue travelling upstream, while small perturbations 
are absorbed. Note that this phenomenon should be 
widespread, since the range [p C v,Pc3] of convectively 
stable traffic can be quite large, so it is likely to ap- 
pear when the freeway flow Q up exceeds the critical 

flow Q e (p c l). 

C. Heterogeneous traffic 

In reality, driver-vehicle units are not identical, but 
behave differently, which may be reflected by individ- 
ual parameter sets. This is easy to do in microscopic 
traffic models, while gas- kinetic and macroscopic mod- 
els require generalizations. Such generalized models 
have been developed for a continuous distribution 
of desired velocities (Paveri-Fontana, 1975; Helbing, 
1995c, d, 1996a, 1997a; Wagner et al., 1996; Wagner, 
1997a, b, c). However, if several driver- vehicle param- 



eters are varied, it is more practical to distinguish 
several vehicle classes only, such as aggressive drivers 
("rabbits") and timid ones ("slugs") (Daganzo, 1995c, 
1999a) or cars and trucks (long vehicles; see, for ex- 
ample, Hclbing, 1996d, 1997a, d; Hoogcndoorn and 
Bovy, 1998a, b, 2000b; Hoogendoorn, 1999; Shvetsov 
and Helbing, 1999). New phenomena arising from the 
heterogeneity of driver-vehicle units are platoon for- 
mation and scattering of flow-density data, as is out- 
lined in the following. 

1. Power-law platoon formation and quenched disorder 

The simplest models for heterogeneity are particle 
hopping models with quenched disorder. For example, 
Evans (1996), Krug and Ferrari (1996), Karimipour 
(1999a, b), as well as Seppalalnen and Krug (1999) 
study a simplified version of a model by Benjamini 
et al. (1996). It corresponds to the one-dimensional 
driven lattice gas known as TASEP, but with particle- 
specific, constant jump rates q a . Since overtaking was 
not allowed, Krug and Ferrari have found a sharp 
phase transition between a low-density regime, where 
all particles are queueing behind the slowest parti- 
cle, and a high-density regime, where the particles are 
equally distributed. While, at low densities, the slow 
particles "feel free traffic" until the critical density 
is reached (cf. the truck curve in Figs. 15b, c, 44a), 
the growth of particle clusters ( "platoons" ) is charac- 
terized by a power-law coarsening. If particles move 
ballistically with individual velocities v a and form a 
platoon when a faster particle reaches a slower one, 
the platoon size n v \{t) grows according to 

n pl (t)~i (m ' +1)/(m ' +2) , (H6) 

where the exponent ml characterizes the distribution 
Po(v) ~ (v — Wmin)" 1 of free velocities in the neighbor- 
hood of the minimal desired velocity u m i n (Ben-Nairn 
et al, 1994; Ben-Nairn and Krapivsky, 1997, 1998, 
1999; Nagatani, 1996c; cf. also Nagatani, 1995a; Na- 
gatani et al, 1998; for spatial inhomogencities see the 
results by Krug, 2000). Above the critical density, the 
differences among fast and slow particles become irrel- 
evant, because there is so little space that all particles 
have to move slower than preferred. The formation of 
platoons has been compared with Bose-Einstein con- 
densation, where the steady-state velocity of the par- 
ticles is analogous to the fugacity of the ideal Bose 
gas (Evans, 1996; Krug and Ferrari, 1996; Chowdhury 
et al, 2000b). 

Platoon formation and power-law coarsening has 
also been found in microscopic models with paral- 
lel update (Ben-Nairn et al, 1994; Nagatani, 1996c; 
Ben-Nairn and Krapivsky, 1997, 1998, 1999; Nagatani 
et al., 1998; see also Fukui and Ishibashi, 1996a; Na- 
gatani, 2000). An example is the Nagel-Schreckenberg 
model with quenched disorder, i.e. vehicle-specific 



slow-down probabilities p a (Ktitarcv et al, 1997; 
Knospe et al, 1999). 

In real traffic, platoons remain limited in size (see 
Sec. II. D). This is probably because of occasional pos- 
sibilities for lane changes on multi-lane roads. A model 
for platoon size distributions has been developed by 
Islam and Consul (1991). 



2. Scattering 

In order to model the observed scattering of con- 
gested flow-density data, Helbing et al. have simu- 
lated traffic with different driver- vehicle units. In con- 
trast to the models discussed in Sec. IV.C.l, they 
used traffic models with a linearly unstable density 
regime, namely the GKT model (Treiber and Hel- 
bing, 1999a), the discrete optimal velocity model (Hel- 
bing and Schreckenberg, 1999), and the IDM (Treiber 
et al, 2000). Since they could determine the time- 
dependent variation of the percentage of long vehi- 
cles in their empirical data, they distinguished two 
vehicle classes, cars and trucks (long vehicles). These 
were characterized by two different sets of parame- 
ters, i.e. two different fundamental diagrams. Con- 
sequently, the effective fundamental diagram was a 
function of the measured proportion of long vehi- 
cles. Performing simulations with the empirically ob- 
tained boundary conditions and truck proportions (see 
Fig. 24), they managed to reproduce the observed 
transition from free to "synchronized" congested flow 
semi-quantitatively, see Fig. 43 (Treiber and Helbing, 
1999a). 

FIG. 43. The displayed points in flow-density space cor- 
respond to empirical 1-minute data (dark crosses) and 
related simulation results (grey boxes) at two different 
cross sections of the freeway. At the section upstream 
of a bottleneck (left), we find "synchronized" congested 
flow, while the close downstream section (right) displays 
homogeneous-in-speed states. The simulations manage to 
reproduce both, the quasi-linear flow-density relation at 
small densities and the scattering over a two-dimensional 
region in the congested regime. For comparison, we have 
displayed the equilibrium flow-density relations for traffic 

consisting of 100% cars (— ), and 100% trucks ( ). The 

respective model parameters were weighted with the mea- 
sured, time-dependent truck fraction, see Fig. 24. (After 
Treiber and Helbing, 1999a; Helbing et al, 2001a, b.) 

In particular, the spatial dependence of the flow- 
density data on the respective freeway cross section 
was well reproduced, including the relaxation to free 
traffic downstream of the bottleneck, which leads to 
homogeneous-in-speed states, i.e. synchronized flow of 
type (ii). Moreover, the flow-density data had a lin- 
ear structure in the regime of free traffic, while they 
were widely scattered in the congested regime. The 
most important requirement for this was a consider- 
able difference in the safe time clearance T a of cars 



and trucks (cf. Fig. 25). This is certainly in agree- 
ment with facts, and it is also compatible with the 
explanation of randomly sloped flow-density changes 
suggested by Banks (see Sec. II. B). If this interpreta- 
tion of the scattering of flow-density data is correct, 
one should observe a smaller variation of flow-density 
data on days when truck traffic is reduced or prohib- 
ited. However, some scattering is always expected due 
to the wide distribution of time headways (see Figs. 6 
and 25). 



D. Multi-lane traffic and synchronization 

1. Gas-kinetic and macroscopic models 

Traffic models have been also developed for multi- 
lane traffic (Gazis et al, 1962; Munjal et al., 1971; 
Munjal and Pipes, 1971; Makigami et al., 1983; Hol- 
land and Woods, 1997; Daganzo, 1997b; Daganzo 
et al, 1997). In gas-kinetic and macroscopic traffic 
models, this is reflected by additional lane-changing 
terms (R0rbech, 1976; Michalopoulos et al, 1984; Hel- 
bing, 1996d, 1997a, d; Helbing and Greiner, 1997; 
Hoogendoorn and Bovy, 1998a, b; Hoogendoorn, 1999; 
Shvctsov and Helbing, 1999; Klar and Wegener, 
1999a, b). Representing the vehicle density in lane i 
by pi(x,t) and the lane-specific average velocity by 
Vi(x, t), the lane-specific density equations of the form 
(110) are complemented by the additional contribu- 
tions 



Pi-i 



Pi 



'i+1 
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(117) 



Herein, l/r^ is the density-dependent lane-changing 
rate from lane i to lane (i + 1), while 1/t[~ de- 
scribes analogous changes to lane (i — 1). Empirical 
measurements of lane-specific data and lane-changing 
rates are rare (Sparman, 1978; Hall and Lam, 1988; 
Chang and Kao, 1991; McDonald et al, 1994; Brack- 
stone and McDonald, 1996; Brackstone et al, 1998). 
Based on his data, Sparman (1978) proposed the re- 
lation l/rf = cfpi(p maK — Pi±i) with suitable lane- 
dependent constants cf. 

The additional terms in the velocity equations in- 
clude 



+^ + ± {V i -i-V i ) + ^{V i+1 -V i ), (118) 

T i-1 T i+1 

which lead to a velocity adaptation in neighboring 
lanes. This is responsible for the synchronization 
among lanes in congested traffic (Lee et al, 1998) and 
allows the treatment of multi-lane traffic by effective 
one-lane models (Shvetsov and Helbing, 1999). Sim- 
ilar adaptation terms are obtained for the variances. 
More detailed calculations yield further terms reflect- 
ing that overtaking vehicles transfer their higher speed 
to the neighboring lane (Helbing, 1997a). 

The lane-specific traffic equations allow us to de- 
rive simplified equations for the whole freeway cross 



Pi+i 



section. It turns out that, compared to the effective 
one-lane model used before, one has to correct the 
expression for the traffic pressure at small densities, 
where the average velocities in the neighboring lanes 
are not synchronized (Hclbing, 1997a, d): 

P(x,t) « p(x,t)[6(x,t) + W - V) 2 )] . (119) 

The density-dependent contribution ((Vi — V) 2 ) re- 
flects the difference in the lane-specific speeds Vi, 
which also contributes to the overall velocity variance 
of vehicles (see Fig. 15a). In contrast, the interaction 
frequency among vehicles, which enters the braking 
interaction term in (86) and V e , depends on the Zone- 
specific variance 6. 

2. Microscopic models and cellular automata 

A microscopic modelling of lane changes (Fox and 
Lchmann, 1967; Levin, 1976; Daganzo, 1981; Mah- 
massani and ShefE, 1981; Ahmed et al, 1996) turns 
out to be an intricate problem for several reasons: 

(i) Because of legal regulations, lane-changing is not 
symmetric, at least in Europe. 

(ii) Drivers accept smaller gaps than normal in cases 
of lane closures, i.e. they behave differently. 

(iii) The minimum gaps required for lane-changing de- 
pend on the own velocity and the speed differences to 
the neighboring vehicles in the adjacent lane. 

(iv) Without assuming fluctuations or heterogeneous 
vehicles, it is hard to reproduce the high rates of lane 
changes, which drivers can reach by means of tricky 
strategies. 

As a consequence, classical gap acceptance models do 
not work very well (Gipps, 1986; Yang and Kout- 
sopoulos, 1996). Benz (1993) and Hclbing (1997a) 
have, therefore, suggested a flexible lane-changing 
criterion. Lane-changing rules for cellular automata 
have been developed as well (Nagatani, 1993, 1994a, 
b, 1996a; Zhang and Hu, 1995; Rickert et al, 1996; 
Chowdhury et al, 1997a; Wagner et al, 1997; Awazu, 
1998; Nagel et al, 1998; Brilon and Wu, 1999). 

In modern lane-changing models, the lane is 
changed if both, an incentive criterion and a safety 
criterion are fulfilled (Nagel et al, 1998). The incen- 
tive criterion checks whether the considered driver a 
could go faster on (one of) the neighboring lane(s). 
The lane, for which the absolute value \f a p\ of the 
repulsive interaction force with the respective leading 
vehicle [3 is smallest, is the preferred one. It is, how- 
ever, reasonable to introduce a certain threshold value 
for lane changing. The safety criterion is fulfilled, if 
the anticipated repulsive interaction forces to neigh- 
boring vehicles on the adjacent lane stay below a cer- 
tain critical value which depends on the friendliness 
of the driver a. The IDM and other models consid- 
ering velocity differences take into account that the 
gap required for lane changing must be greater if the 



new following vehicle approaches fast, i.e. if its rela- 
tive velocity with respect to the lane-changing vehicle 
is large. 

To obtain reasonable lane changing rates, one has 
to either work with a probabilistic model or distin- 
guish different vehicle types like cars and trucks. In 
a small density range around 25 vehicles per kilome- 
ter and lane, simulations by Helbing and Huberman 
(1998) indicated a transition to a coherent movement 
of cars and trucks, which they could confirm empir- 
ically (cf. Figs. 15b, c and 44a). In this state, fast 
and slow vehicles drive with the same speed, i.e. they 
are moving like a solid block. This is caused by a 
breakdown of the lane-changing rate (see Fig. 44b), 
when the freeway becomes crowded, so that drivers 
do not anymore find sufficiently large gaps for over- 
taking (see Fig. 44c). However, the coherent state is 
destroyed when traffic flow becomes unstable. Then, 
vehicle gaps are again widely scattered (cf. Fig. 44d), 
and lane changing is sometimes possible. 

FIG. 44. The transition to coherent vehicle traffic, see 
(a), is related to the breakdown of the lane-changing rate 
at densities above 20 vehicles per kilometer, see (b). This 
is related to the rapid decay of the proportion of successful 
lane changes (i.e. the quotient between actual and desired 
lane changes), see (c). (d) Opportunities for lane-changing 
are rapidly diminished when gaps corresponding to about 
twice the safe headway required for lane-changing cease to 
exist, i.e. when the quotient of the standard deviation of 
vehicle headways to their mean value is small. (After Hel- 
bing and Huberman, 1998; see also Helbing, 2001). The 
breakdown of the lane-changing rate seems to imply a de- 
coupling of the lanes, i.e. an effective one-lane behavior. 
However, this is already the result of a self-organization 
process based on two-lane interactions, since any signif- 
icant perturbation of the solid-block state (like different 
velocities in the neighbouring lanes) will cause frequent 
lane changes (Shvetsov and Helbing, 1999). By filling large 
gaps, the gap distribution is considerably modified. This 
eventually reduces possibilities for lane changes, so that 
the solid-block state is restored. 

Properties of multi-lane traffic have recently found 
a considerable interest, as is documented by fur- 
ther studies (Chowdhury et al, 1997b; Bclitsky et al, 
2001; see also Goldbach et al, 2000; Lubashevsky and 
Mahnke, 2000). 

E. Bi-directional and city traffic 

Traffic models have also been developed for bi- 
directional traffic (Dutkiewicz et al, 1995; Lee et al, 
1997; Simon and Gutowitz, 1998; Fouladvand and 
Lee, 1999) and city traffic (see, for example, TRAN- 
SYT by Robertson, 1969a, b; Schmidt-Schmiedebach, 
1973; Herman and Prigoginc, 1979; SCOOT by Hunt 
et al, 1982; Cremer and Ludwig, 1986; Williams 
et al, 1987; Mahmassani et al, 1990; MAKSIMOS by 
Putensen, 1994; MIXSIM by Hoque and McDonald, 



1995; PADSIM by Peytchev and Bargiela, 1995). In- 
terested readers should consult the reviews by Papa- 
georgiou (1999), Chowdhury et al. (2000b), and the 
corresponding chapters in Traffic and Granular Flow 
'97, '99 (Schreckenberg and Wolf, 1998; Helbing et al, 
2000c). Here, I will be very short. 

Traffic dynamics in cities is quite different from free- 
ways, since it is to a large extent determined by the 
intersections. The first model for city traffic in the 
physics literature, the "BML model", goes back to 
Biham et al. (1992). It is based on a two-dimensional 
square lattice, in which each site represents a crossing 
of north-bound and east-bound traffic. The sites can 
be either empty, or occupied by a north-bound or east- 
bound vehicle. East-bound vehicles are synchronously 
updated every odd time step, while north-bound vehi- 
cles are updated every even time step, reflecting syn- 
chronized and periodic traffic signals at the crossings. 
During the parallel updates, a vehicle moves ahead 
by one site, if this is empty, otherwise it has to wait. 
Biham et al. (1992) observed a first-order transition 
from free traffic to jammed traffic with zero veloc- 
ity throughout the system, which appears at a finite 
density due to a gridlock of north-bound and east- 
bound traffic along diagonal lines (see also Fukui and 
Ishibashi, 1996b). This first-order transition is also 
characteristic for most other models of city traffic. For 
a detailled discussion of the related literature see the 
review by Chowdhury et al. (2000b). 

F. Effects beyond physics 

There are some effects which are not at all or only 
partly represented by the traffic models described in 
Sec. Ill: 

(i) the anticipation behavior of drivers (Ozaki, 1993; 
Wagner, 1998a; Lenz et al., 1999; Knospe et al, 
2000a, b), 

(ii) the reaction to winkers, 

(iii) the avoidance of driving side by side, 

(iv) the velocity adaptation to surrounding traffic 
(Herman et al., 1973), 

(v) the tolerance of small time clearances for a limited 
time (Treiterer and Myers, 1974), 

(vi) the reduction of gaps to avoid vehicles cutting in 
(Daganzo, 1999a), 

(vii) the avoidance of the truck lane, even if it is un- 
derutilized (Brackstonc and McDonald, 1996; Helbing, 
1997a, b; "Los Gatos effect" : Daganzo, 1997c, 1999a; 
Helbing and Huberman, 1998), 

(viii) motivations for lane changing (Rcdclmeicr and 
Tibshirani, 1999), 

(ix) effects of attention (Michaels and Cozan, 1962; 
Todosicv, 1963; Todosiev and Barbosa, 1963/64; 
Michaels, 1965), 

(x) changes of driver strategies (Migowsky et al, 
1994), 

(xi) effects of road conditions and visibility, etc. 



Certainly, such kinds of effects can be also modelled 
in a mathematical way, and this has partly been done 
in so-called sub- microscopic models of vehicle dynam- 
ics (Wiedemann, 1974) and vehicle simulators such 
as VISSIM (Fellendorf, 1996) or PELOPS (Ludmann 
et al., 1997). However, most of the above effects are 
hard to test empirically and difficult to verify or fal- 
sify. Therefore, we are leaving the realm of physics, 
here, and eventually entering the realm of driver psy- 
chology (Daganzo, 1999a). 

G. Traffic control and optimization 

Traffic control and optimization is an old and 
large field (see, e.g., Cremer, 1978; Cremer and 
May, 1985; Smulders, 1989; Dougherty, 1995; Kiihne, 
1993; Kiihne et al, 1995, 1997; Taale and Middel- 
ham, 1995; Schweitzer et al, 1998a), which can be 
certainly not discussed in detail, here. Overviews 
of optimization approaches by physicists can be 
found in the Proceedings of the Workshops on Traf- 
fic and Granular Flow (Wolf et al, 1996; Schreck- 
enberg and Wolf, 1998; Helbing et al, 2000c; see 
also Claus et al, 1999), while the readers inter- 
ested in the engineering literature are referred to 
the books and review articles by Pagageorgiou (1995, 
1999) as well as Lapierre and Steierwald (1987). 
These also discuss the available software packages for 
the simulation and optimization of street networks 
and/or signalling of traffic lights such as TRANSYT 
(Robertson, 1969a, b), SCOOT (Hunt et al, 1982), 
OPAC (Gartner, 1983), PRODYN (Farges et al, 
1984), CRONOS (Boillot et al, 1992), and many 
others. I should also mention the cellular-automata- 
based projects OLSIM (Kaumann et al, 2000; see 
http : //www . traf f ic . uni-duisburg . de/DLSIM/) 
and TRANSIMS (Nagel and Barrett, 1997; Rick- 
ert and Nagel, 1997; Simon and Nagel, 1998; see 
http: //www-trans ims. tsasa.lanl. gov). A com- 
plementary 

list of microsimulation tools can be found at http:// 
www. its . leeds . ac .uk/smartest/links .html. 

While optimization of city traffic is based on a syn- 
chronization of traffic lights, most optimization ap- 
proaches for freeways are based on a homogenization 
of vehicle traffic, as stop-and-go traffic, jams, and con- 
gested traffic are associated with reduced efficiency 
due to the capacity drop (see, e.g. Huberman and Hel- 
bing, 1999). Remember that this capacity drop is often 
triggered by perturbations in the flow (cf. Sec. IV. B), 
which should, therefore, be suppressed by technical 
measures such as intelligent speed limits, adaptive on- 
ramp controls, traffic-dependent re-routing, and driver 
assistance systems (sec, e.g., Treiber and Helbing, 
2001). 



V. PEDESTRIAN TRAFFIC 

Pedestrian crowds have been empirically stud- 
ied for more than four decades now (Hankin and 
Wright, 1958; Oeding, 1963; Hoel, 1968; Older, 1968; 
Navin and Wheeler, 1969; Carstens and Ring, 1970; 
O'Flaherty and Parkinson, 1972; Weidmann, 1993). 
The evaluation methods applied were based on direct 
observation, photographs, and time-lapse films. Apart 
from behavioral investigations (Hill, 1984; Batty, 
1997), the main goal of these studies was to develop 
a level-of-service concept (Fruin, 1971; Polus et al, 
1983; Mori, and Tsukaguchi, 1987), design elements 
of pedestrian facilities (Schubert, 1967; Bocminghaus, 
1982; Pauls, 1984; Whyte, 1988; Helbing, 1997a; Hel- 
bing et al, 2001c), or planning guidelines (Kirsch, 
1964; Predtetschenski and Milinski, 1971; Transporta- 
tion Research Board, 1985; Davis and Braaksma, 
1988; Brilon et al, 1993). The latter have usually the 
form of regression relations, which are, however, not 
very well suited for the prediction of pedestrian flows 
in pedestrian zones and buildings with an exceptional 
architecture, or in extreme conditions such as evacu- 
ation (see Sec. V.D). Therefore, a number of simula- 
tion models have been proposed, e.g. queueing mod- 
els (Yuhaski and Macgregor Smith, 1989; Roy, 1992; 
L0vas, 1994), transition matrix models (Garbrecht, 
1973), and stochastic models (Mayne, 1954; Ashford 
et al., 1976), which are partly related to each other. In 
addition, there are models for the route choice behav- 
ior of pedestrians (Borgers and Timmermans, 1986a, 
1986b; Timmermans et al., 1992; Helbing, 1992a). 

None of these concepts adequately takes into ac- 
count the self-organization effects occuring in pedes- 
trian crowds. These may, however, lead to unexpected 
obstructions due to mutual disturbances of pedestrian 
flows. More promising with regard to this is the ap- 
proach by Henderson. He conjectured that pedestrian 
crowds behave similar to gases or fluids (Henderson, 
1971, 1974; Henderson and Lyons, 1972; Henderson 
and Jenkins, 1973; see also Hughes, 2000, 2001). This 
could be partly confirmed (see Sec. V.A). However, a 
realistic gas- kinetic or fluid-dynamic theory for pedes- 
trians must contain corrections due to their particu- 
lar interactions (i.e. avoidance and deceleration ma- 
neuvers) which, of course, do not conserve momen- 
tum and energy. Although such a theory can be actu- 
ally formulated (Helbing, 1992a, b; Hoogcndoorn and 
Bovy, 2000a), for practical applications a direct sim- 
ulation of individual pedestrian motion is favourable, 
since this is more flexible. As a consequence, current 
research focusses on the micro simulation of pedestrian 
crowds, which also allows one to consider incoordina- 
tion by excluded volume effects related to the discrete, 
"granular" structure of pedestrian flows. In this con- 
nection, a behavioral force model of individual pedes- 
trian dynamics has been developed (Helbing, 1991, 
1996d, 1997a, 1998c; Helbing et al, 1994; Helbing and 
Molnar, 1995, 1997; Molnar, 1996a, b; Helbing and 



Vicsek, 1999; Helbing et al, 2000a, b) (see Sees. I.C.I 
and V.B). A discrete and simple forerunner of this 
model was proposed by Gipps and Marksjo (1985). I 
also like to mention recent cellular automata of pedes- 
trian dynamics (Bolay, 1998; Blue and Adler, 1998, 
1999; Fukui and Ishibashi, 1999a, b; Muramatsu et al, 
1999; Muramatsu and Nagatani, 2000a, b; Klupfcl 
et al, 2000; Burstedde et al, 2001), emergency and 
evacuation models (Drager et al, 1992; Ebihara et al, 
1992; Ketchell et al, 1993; Okazaki and Matsushita, 
1993; Still, 1993, 2000; Thomson and Marchant, 1993; 
L0vas, 1998; Hamacher and Tjandra, 2000; Klupfcl 
et al, 2000), and Al-based models (Gopal and Smith, 
1990; Reynolds, 1994, 1999; Schelhorn et aZ.1999). 

A. Observations 

Together with my collaborators I have investigated 
pedestrian motion for several years and evaluated a 
number of video films. Despite the sometimes more 
or less "chaotic" appearance of individual pedestrian 
behavior, one can find regularities, some of which be- 
come best visible in time-lapse films like the ones pro- 
duced by Arns (1993). While describing these, I also 
summarize results of other pedestrian studies and ob- 
servations (after Helbing, 1997a, 1998c; Helbing et al, 
2001c): 

(i) Pedestrians feel a strong aversion to taking detours 
or moving opposite to their desired walking direction, 
even if the direct way is crowded. However, there is 
also some evidence that pedestrians normally choose 
the fastest route to their next destination, but not the 
shortest one (Ganem, 1998). In general, pedestrians 
take into account detours as well as the comfort of 
walking, thereby minimizing the effort to reach their 
destination (Helbing et al, 1997a, b). Their ways can 
be approximated by polygons. 

(ii) Pedestrians prefer to walk with an individual de- 
sired speed, which corresponds to the most comfort- 
able (i.e. least energy-consuming) walking speed (see 
Weidmann, 1993) as long as it is not necessary to go 
faster in order to reach the destination in time. The 
desired speeds within pedestrian crowds are Gaus- 
sian distributed with a mean value of approximately 
1.34m/s and a standard deviation of about 0.26m/s 
(Henderson, 1971). However, the average speed de- 
pends on the situation (Predtetschenski and Milinski, 
1971), sex and age, the time of the day, the purpose of 
the trip, the surrounding, etc. (see Weidmann, 1993). 

(iii) Pedestrians keep a certain distance to other pedes- 
trians and borders (of streets, walls, and obstacles; see 
Transportation Research Board, 1985; Brilon et al, 
1993). This distance is smaller the more a pedestrian 
is in a hurry, and it decreases with growing pedestrian 
density. Resting individuals (waiting on a railway plat- 
form for a train, sitting in a dining hall, or lying at 
a beach) are uniformly distributed over the available 
area if there are no acquaintances among the individ- 



uals. Pedestrian density increases (i.e. interpersonal 
distances lessen) around particularly attractive places. 
It decreases with growing velocity variance (e.g. on a 
dance floor, sec Hclbing, 1992b, 1997a). Individuals 
knowing each other may form groups which are enti- 
ties behaving similar to single pedestrians. Group sizes 
are Poisson distributed (Coleman, 1964; Coleman and 
James, 1961; Goodman, 1964). 

(iv) In situations of escape panic (Hclbing et at, 
2000b), individuals are often nervous, and some- 
times they develop blind actionism. Moreover people 
move or try to move considerably faster than nor- 
mal (Predtetschenski and Milinski, 1971). Individu- 
als start pushing, and interactions among people be- 
come physical in nature. Moving and, in particular, 
passing of a bottleneck becomes incoordinated (Mintz, 
1951). At exits, arching and clogging are observed 
(Predtetschenski and Milinski, 1971), and jams are 
building up (Mintz, 1951). The physical interactions 
in the jammed crowd add up and cause dangerous 
pressures up to 4,500 Ncwtons per meter (Elliott and 
Smith, 1993; Smith and Dickie, 1993), which can bend 
steel barriers or tear down brick walls. Escape is fur- 
ther slowed down by fallen or injured people turning 
into "obstacles" . Finally, people tend to show herding 
behavior, i.e., to do what other people do (Quaran- 
telli, 1957; Keating, 1982). As a consequence, alterna- 
tive exits are often overlooked or not efficiently used 
in escape situations (Keating, 1982; Elliott and Smith, 
1993). 

(v) At medium and high pedestrian densities, the mo- 
tion of pedestrian crowds shows some striking analo- 
gies with the motion of gases, fluids, and granular 
flows. For example, I found that footprints of pedes- 
trians in snow look similar to streamlines of flu- 
ids (Helbing, 1992a). At borderlines between oppo- 
site directions of walking one can observe "viscous 
fingering" (Kadanoff, 1985; Stanley and Ostrowsky, 
1986). The emergence of pedestrian streams through 
standing crowds (Arns, 1993; Helbing 1997a, Fig. 2.6; 
1998c, Fig. 1; Helbing et ai, 2001c, Fig. 1) appear 
analogous to the formation of river beds (St0lum, 
1996; Rodrfguez-Iturbe and Rinaldo, 1997; Caldarelli, 
2000). Moreover, similar to segregation or stratifi- 
cation phenomena in granular media (Santra et al., 
1996; Makse et al, 1997), pedestrians spontaneously 
organize in lanes of uniform walking direction, if the 
pedestrian density is high enough (Ocding, 1963; Hel- 
bing, 1991, 1997a, Fig. 2.5; 1998c, Fig. 2; Helbing 
et al, 2001c, Fig. 2). At bottlenecks (e.g. corridors, 
staircases, or doors), the passing direction of pedes- 
trians oscillates (Helbing et al, 1994; Helbing and 
Molnar, 1995). This may be compared to the "saline 
oscillator" (Yoshikawa et al., 1991) or the granular 
"ticking hour glass" (Wu et al, 1993; Pennec et al, 
1996). Finally, one can find the propagation of shock 
waves in dense pedestrian crowds pushing forward (see 
also Virklcr and Elayadath, 1994). The arching and 
clogging in panicking crowds (Hclbing et al, 2000b) is 



similar to the outflow of rough granular media through 
small openings (Ristow and Herrmann, 1994; Wolf and 
Grassberger, 1997). 

B. Behavioral force model of pedestrian dynamics 

For reliable simulations of pedestrian crowds we do 
not need to know whether a certain pedestrian, say, 
turns to the right at the next intersection. It is suf- 
ficient to have a good estimate what percentage of 
pedestrians turns to the right. This can be either em- 
pirically measured or calculated by means of route 
choice model like the one by Borgers and Timmer- 
mans (1986a, b). In some sense, the uncertainty about 
the individual behaviors is averaged out at the macro- 
scopic level of description, as in fluid dynamics. 

Nevertheless, I will now focus on a more flexible 
microscopic simulation approach based on the behav- 
ioral force concept outlined in Sec. I.C.I. Particular 
advantages of this approach are the consideration of 
(i) excluded volume effects due to the granular struc- 
ture of a pedestrian crowd and (ii) the flexible usage 
of space by pedestrians, requiring a (quasi-)continous 
treatment of motion. It turns out that these points are 
essential to reproduce the above mentioned phenom- 
ena in a natural way. 

As outlined in Sec. I.C.I, we describe the differ- 
ent motivations of and influences on a pedestrian a 
by separate force terms. First of all, the decisions 
and intentions where to go and how fast are reflected 
by the driving term v^{t)e { ^{t)/T a , in which v^(t) is 
the desired velocity and e Q a (t) the desired direction of 
motion. The interactions among pedestrians are de- 
scribed by forces 

= + /$(*), (120) 

which consist of a socio-psychological contribution 
f S af}{t) and physical interactions f a h p(t), if pedestrians 
come too close to each other. The socio-psychological 
component reflects the tendency of pedestrians to keep 
a certain distance to other pedestrians and may be de- 
scribed by a force of the form 

/a/3 (*) = A a exp[(r Q/3 - d a p) I B a ]n al j 

x(A a + (l-A a ) 1 + C °^ ) ). (121) 

Herein, A a denotes the interactions strength and B a 
the range of the repulsive interaction, which arc indi- 
vidual parameters and partly dependent on cultural 
conventions. d a p(t) = \\x a {f) — xp{t)\\ is the distance 
between the centers of mass of pedestrians a and /3, 
r af3 = {f a + r@) the sum of their radii r a and rp, and 

nafiit) = (n a(3 (t),nl (t)) = ^ "'^ 

d a p{t) 

the normalized vector pointing from pedestrian (3 to 
a. Finally, with the choice X a < 1, we can take into 



account the anisotropic character of pedestrian inter- 
actions, as the situation in front of a pedestrian has a 
larger impact on his or her behavior than things hap- 
pening behind. The angle (fa/sit) denotes the angle 
between the direction e a (t) = v a (t)/\\v a (t)\\ of mo- 
tion and the direction —n a g(t) of the object exerting 
the repulsive force, i.e. cos<p a p(t) = —n a p(t) ■ e a (t). 
One may, of course, take into account other details 
such as a velocity-dependence of the forces and non- 
circular shaped pedestrian bodies, but this does not 
have qualitative effects on the dynamical phenomena 
resulting in the simulations. In fact, most observed 
self-organization phenomena are quite insensitive to 
the specification of the interaction forces, as differ- 
ent studies have shown (Hclbing and Molnar, 1995; 
Molnar, 1996a; Helbing and Vicsek, 1999; Helbing 
et al, 2000a, b). 

The physical interaction plays a role only when 
pedestrians have physical contact with each other, 
i.e. if r a fj > d a p. In this case, we assume two addi- 
tional forces inspired by granular interactions (Ristow 
and Herrmann, 1994; Wolf and Grassberger, 1997): 
First, a "body force" k(r a p — d a /3)n a p counteract- 
ing body compression and a "sliding friction force" 
n(r a /3 — d a g) Avg Q t a/ 3 impeding relative tangential 
motion: 

= fc0 ( r a/3 - d a p)n al3 

+ nQ(r af3 - d a0 )Avl a t a0 , (122) 

where the function 6(2) is equal to its argument z, 
if z > 0, otherwise 0. Moreover, t a p = {—n a g-,n a g) 
means the tangential direction and Avjj a = (vp — v a )- 
t a g the tangential velocity difference, while k and n 
represent large constants. Strictly speaking, friction 
effects already set in before pedestrians touch each 
other, because of the socio-psychological tendency not 
to pass other individuals with a high relative velocity, 
when the distance is small. This is, however, not im- 
portant for the effects we are going to discuss later 
on. 

The interaction with the boundaries of walls and 
other obstacles is treated analogously to pedestrian in- 
teractions, i.e., if d ab (t) means the distance to bound- 
ary b, n ab (t) denotes the direction perpendicular to 
it, and t ab (t) the direction tangential to it, the corre- 
sponding interaction force with the boundary reads 

fab = {A a exp[(r Q - d ah )/B a ] + kG(r a - d ab )} n ab 
- n<d(r a - d ab )(v a -t ab )t ab . (123) 

Moreover, we may also take into account time- 
dependent attractive interactions towards window dis- 
plays, sights, or special attractions i by social forces of 
the type (121). However, in comparison with repulsive 
interactions, the corresponding interaction range B ai 
is usually larger and the strength parameter A ai (t) 
typically smaller, negative, and time-dependent. Ad- 
ditionally, the joining behavior (Dewdney, 1987) of 



families, friends, or tourist groups can be reflected 
by forces of the type /„"(t) = — C a pn a f)(t), which 
guarantee that acquainted individuals join again, af- 
ter they have accidentally been separated by other 
pedestrians. Finally, we can take into account un- 
systematic variations of individual behavior and vol- 
untary deviations from the assumed behavioral laws 
by a fluctuation term £ Q (t). Therefore, the resulting 
pedestrian model corresponds to equation (3), but 
the sum of interaction forces E/3(^a) fap{t) is re- 
placed by the sum of all the above mentioned attrac- 
tive and repulsive interaction and boundary forces. 
The resulting expression looks rather complicated, but 
in the following, we will drop attraction effects and 
assume X a = for simplicity, so that the interac- 
tion forces become isotropic and conform with New- 
ton's third law. Moreover, the physical interactions 
are mostly relevant in panic situations, so that we 
normally just have repulsive social and boundary in- 
teractions: E/3(^a) /a|(*) + Eb fab(t). 

C. Self-organization phenomena and noise-induced 
ordering 

Despite its simplifications, the behavioral force 
model of pedestrian dynamics describes a lot of ob- 
served phenomena quite realistically. Especially, it 
allows one to explain various self-organized spatio- 
temporal patterns that are not externally planned, 
prescribed, or organized, e.g. by traffic signs, laws, or 
behavioral conventions. Instead, the spatio-temporal 
patterns discussed below emerge due to the non- 
linear interactions of pedestrians even without assum- 
ing strategical considerations, communication, or imi- 
tative behavior of pedestrians. All these collective pat- 
terns of motion arc symmetry-breaking phenomena, 
although the model was formulated completely sym- 
metric with respect to the right-hand and the left- 
hand side (Helbing and Molnar, 1997; Helbing, 1997a, 
1998c; Helbing et al., 2001c). 

1. Segregation 

The microsimulations reproduce the empirically ob- 
served formation of lanes consisting of pedestrians 
with the same desired walking direction (Helbing 
et al., 1994; Hclbing and Molnar, 1995, 1997; Hel- 
bing, 1996d, 1997a; Molnar, 1996a, b; Helbing and 
Vicsek, 1999; Helbing et al., 2001c; see Fig. 45). For 
open boundary conditions, these lanes are dynami- 
cally varying. Their number depends on the width of 
the street (Helbing and Molnar, 1995; Helbing, 1997a), 
on pedestrian density, and on the noise level. Interest- 
ingly, one finds a noise-induced ordering (Hclbing and 
Vicsek, 1999; Helbing and Platkowski, 2000): Com- 
pared to small noise amplitudes, medium ones result 
in a more pronounced segregation (i.e., a smaller num- 



ber of lanes), while large noise amplitudes lead to a 
"freezing by heating" effect (see Sec. V.D.I). 

FIG. 45. Formation of lanes in initially disordered 
pedestrian crowds with opposite walking directions 
(after Helbing et ai, 2000a; Helbing, 2001; cf. also 
Helbing et ai, 1994; Helbing and Molnar, 1995). 
White disks represent pedestrians moving from left 
to right, black ones move the other way round. 
Lane formation does not require the periodic bound- 
ary conditions applied above, see the Java applet 
http : //www . helbing . org/Pedestr ians/ Corridor . html. 

The conventional interpretation of lane formation 
is: Pedestrians tend to walk on the side prescribed in 
vehicular traffic. However, the above model can ex- 
plain lane formation even without assuming a prefer- 
ence for any side (Helbing and Vicsek, 1999; Helbing 
et ai, 2000a). The most relevant point is the higher 
relative velocity of pedestrians walking in opposite di- 
rections. As a consequence, they have more frequent 
interactions until they have segregated into separate 
lanes. The resulting collective pattern of motion min- 
imizes the frequency and strength of avoidance ma- 
neuvers, if fluctuations are weak. Assuming identical 
desired velocities u° = v , the most stable configu- 
ration corresponds to a state with a minimization of 
interactions 

~jf E T f"P ■ e « M Jf E( y o - v a ■ e° ) = «o(l - E) . 
It is related with a maximum efficiency 
N ^ v 

a 

of motion (Helbing and Vicsek, 1999). The efficiency 
E with < E < 1 (where N — Y^ Q 1 is the respec- 
tive number of pedestrians a) describes the average 
fraction of the desired speed vo with which pedes- 
trians actually approach their destinations. That is, 
lane formation "globally" maximizes the average ve- 
locity into the respectively desired direction of mo- 
tion, although the model does not even assume that 
pedestrians would try to optimize their behavior lo- 
cally. This is a consequence of the symmetrical in- 
teractions among pedestrians with opposite walking 
directions. One can even show that a large class of 
driven many-particle systems, if they self-organize at 
all, tend to globally optimize their state (Helbing and 
Vicsek, 1999). This optimal self- organization was also 
observed in the formation of coherent motion in a sys- 
tem of cars and trucks, which was related with a min- 
imization of the lane changing rate (see Sec. IV. D. 2). 
Another example is the emergence of optimal way sys- 
tems by bundling of trails (see Fig. 1). 



2. Oscillations 

In simulations of bottlenecks like doors, oscillatory 
changes of the passing direction are observed, if peo- 
ple do not panic, see Fig. 46 (Helbing et ai, 1994; Hel- 
bing and Molnar, 1995, 1997; Helbing, 1996d, 1997a; 
Molnar, 1996a, b; Helbing et ai, 2001c). 

FIG. 46. Oscillations of the passing direction at a bottle- 
neck (after Helbing, 2001; cf. also Helbing et al, 1994; Hel- 
bing and Molnar, 1995). Dynamic simulations are available 
at http : //www . helbing . org/Pedestrians/Door . html. 

The mechanism leading to alternating flows is the 
following: Once a pedestrian is able to pass the nar- 
rowing, pedestrians with the same walking direction 
can easily follow. Hence, the number and "pressure" of 
waiting and pushing pedestrians becomes less than on 
the other side of the narrowing where, consequently, 
the chance to occupy the passage grows. This leads 
to a deadlock situation which is followed by a change 
in the passing direction. Capturing the bottleneck is 
easier if it is broad and short so that the passing direc- 
tion changes more frequently. Therefore, two separate 
doors close to the walls are more efficient than one 
single door with double width. By self-organization, 
each door is, then, used by one walking direction for a 
long time, see Fig. 47 (Molnar, 1996a, b; Helbing and 
Molnar, 1997; Helbing, 1997a, 1998c; Helbing et ai, 
2001c). The reason is that pedestrians, who pass a 
door, clear the way for their successors, similar to lane 
formation. 

FIG. 47. If two alternative passageways are available, 
pedestrians with opposite walking direction use different 
doors as a result of self-organization (after Molnar, 1996a, 
b; Helbing and Molnar, 1997; Helbing, 1997a, 1998c; Hel- 
bing et ai, 2001c). 

3. Rotation 

At intersections one is confronted with various al- 
ternating collective patterns of motion which are very 
short-lived and unstable. For example, phases dur- 
ing which the intersection is crossed in "vertical" or 
"horizontal" direction alternate with phases of tem- 
porary roundabout traffic, see Fig. 48 (Helbing et ai, 
1994; Helbing, 1996d, 1997a; Molnar, 1996a, b; Hel- 
bing and Molnar, 1997; Helbing et at, 2001c). This 
self-organized round-about traffic is similar to the 
emergent rotation found for self-driven particles (Du- 
parcmeur et ai, 1995). It is connected with small de- 
tours, but decreases the frequency of necessary decel- 
eration, stopping, and avoidance maneuvers consid- 
erably, so that pedestrian motion becomes more effi- 
cient on average. The efficiency of pedestrian flow can 
be considerably increased by putting an obstacle in 
the center of the intersection, since this favours the 



smooth roundabout traffic compared with the com- 
peting, inefficient patterns of motion. 

FIG. 48. Self-organized, short-lived roundabout traffic 
in intersecting pedestrian streams (from Helbing, 1996d, 
1997a; Helbing and Molnar, 1997; Helbing et al, 2001c; 
see also Helbing et al, 1994; Molnar, 1996a, b). 

D. Collective phenomena in panic situations 

The behavior of pedestrians in panic situations (e.g. 
in certain cases of emergency evacuation) displays 
some characteristic features: 

(i) People are getting nervous, resulting in a higher 
level of fluctuations. 

(ii) They are trying to escape from the cause of panic, 
which can be reflected by a significantly higher desired 
velocity. 

(iii) Individuals in complex situations, who do not 
know what is the right thing to do, orient at the ac- 
tions of their neighbours, i.e. they tend to do what 
other people do. We will describe this by an addi- 
tional herding interaction, but attractive interactions 
have probably a similar effect. 

We will now discuss the fundamental collective ef- 
fects which fluctuations, increased desired velocities, 
and herding behavior can have. In contrast to other 
approaches, we do not assume or imply that individ- 
uals in panic or emergency situations would behave 
relentless and asocial, although they sometimes do. 

1. "Freezing by heating" 

The effect of getting nervous has been investigated 
by Helbing et al. (2000a). Let us assume that the in- 
dividual level of fluctuations is given by 

Va = (1 - n' a )r) + n' Q ?7 max , (125) 

where n' a with < n' a < 1 measures the nervous- 
ness of pedestrian a. The parameter 770 means the 
normal and 77 max the maximum fluctuation strength. 
It turned out that, at sufficiently high pedestrian den- 
sities, lanes were destroyed by increasing the fluctua- 
tion strength (which is analogous to the temperature). 
However, instead of the expected transition from the 
"fluid" lane state to a disordered, "gaseous" state, 
Helbing et al. found the formation of a solid state. 
It was characterized by a blocked situation with a 
regular (i.e. "crystallized" or "frozen") structure so 
that this paradoxial transition was called "freezing by 
heating" (see Fig. 49). Notably enough, the blocked 
state had a higher degree of order, although the in- 
ternal energy was increased and the resulting state 
was metastable with respect to structural perturba- 
tions such as the exchange of oppositely moving par- 
ticles. Therefore, "freezing by heating" is just oppo- 
site to what one would expect for equilibrium sys- 
tems, and different from fluctuation-driven ordcrding 



phenomena in metallic glasses and some granular sys- 
tems (Gallas et al, 1992; Umbanhowar et al., 1996; 
Rosato et al., 1987), where fluctuations lead from a 
disordered metastable to an ordered stable state. A 
model for related noise-induced ordering processes has 
been developed by Helbing and Platkowski (2000). 

The precondition for the unusual freezing-by- 
heating transition are the additional driving term and 
the dissipative friction. Inhomogeneities in the chan- 
nel diameter or other impurities which temporarily 
slow down pedestrians can further this transition at 
the respective places. Finally note that a transition 
from fluid to blocked pedestrian counter flows is also 
observed, when a critical particle density is exceeded 
(Muramatsu et al., 1999; Helbing et al, 2000a). 

FIG. 49. Noise-induced formation of a crystallized, 
"frozen" state in a periodic corridor used by oppositely 
moving pedestrians (after Helbing et al, 2000a; Helbing, 
2001). 

2. "Faster-is-slower effect" 

The effect of increasing the desired velocity was 
studied by Helbing et al. (2000b). They found that the 
simulated outflow from a room is well-coordinated and 
regular, if the desired velocities w° = vq are normal. 
However, for desired velocities above 1.5 m/s, i.e. for 
people in a rush (as in many panic situations), one ob- 
serves an irregular succession of arch-like blockings of 
the exit and avalanche-like bunches of leaving pedes- 
trians, when the arches break up (see Fig. 50). This 
phenomenon is compatible with the empirical obser- 
vations mentioned above and comparable to intermit- 
tent clogging found in granular flows through funnels 
or hoppers (Rlstow and Herrmann, 1994; Wolf and 
Grassberger, 1997). 

FIG. 50. Panicking pedestrians often come so close to 
each other, that their physical contacts lead to the build 
up of pressure and obstructing friction effects. This results 
in temporary arching and clogging related with inefficient 
and irregular outflows. (After Helbing et al, 2000b; Hel- 
bing, 2001.) 

As clogging is connected with delays, trying to move 
faster (i.e., increasing v°) can cause a smaller average 
speed of leaving, if the friction parameter k is large 
(Helbing et al., 2000b). This "faster-is-slower effect" 
is particularly tragic in the presence of fires, where 
fleeing people sometimes reduce their own chances 
of survival. As a consequence, models for everyday 
pedestrian streams are not very suitable for realistic 
simulations of emergency situations, which require at 
least modified parameter sets corresponding to less ef- 
ficient pedestrian behavior. Related fatalities can be 
estimated by the number of pedestrians reached by 
the fire front (see http://angel.elte.hu/~panic/). 



Since the interpersonal friction was assumed to 
have, on average, no deceleration effect in the crowd, if 
the boundaries are sufficiently remote, the arching un- 
derlying the clogging effect requires a combination of 
several effects: First, slowing down due to a bottleneck 
such as a door, and second, strong inter-personal fric- 
tion, which becomes dominant when pedestrians get 
too close to each other. It is, however, noteworthy that 
the faster-is-slower effect also occurs when the sliding 
friction force changes continuously with the distance 
rather than being "switched on" at a certain distance 
rp as in the model above. 

The danger of clogging can be minimized by avoid- 
ing bottlenecks in the construction of stadia and 
public buildings. Notice, however, that jamming can 
also occur at widenings of escape routes! This comes 
from disturbances due to pedestrians, who expand 
in the wide area because of their repulsive interac- 
tions or try to overtake each other. These squeeze 
into the main stream again at the end of the widen- 
ing, which acts like a bottleneck and leads to jam- 
ming. Significantly improved outflows in panic situ- 
ations can be reached by columns placed asymmet- 
rically in front of the exits, which reduce the pres- 
sure at the door and, thereby, also reduce injuries (see 
http : //angel . elte .hu/ "panic/). 

3. "Phantom panics" 

Sometimes, panics have occured without any com- 
prehensible reasons such as a fire or another threaten- 
ing event. Due to the "faster-is-slower effect", pan- 
ics can be triggered by small pedestrian counter- 
flows (Elliott and Smith, 1993), which cause delays to 
the crowd intending to leave. Consequently, stopped 
pedestrians in the back, who do not sec the reason 
for the temporary slowdown, are getting impatient 
and pushy. In accordance with observations (Helbing, 
1991, 1997a), one may describe this by increasing the 
desired velocity, for example, by the formula 

«S(*) = [i-«a(*K(o)+<(*)«r c - (126) 

Herein, w™ ax is the maximum desired velocity and 
v®(0) the initial one, corresponding to the expected 
velocity of leaving. The time-dependent parameter 

<m = i - §§ (i2T) 

reflects the nervousness, where v a (t) denotes the av- 
erage speed into the desired direction of motion. Alto- 
gether, long waiting times increase the desired veloc- 
ity, which can produce inefficient outflow. This further 
increases the waiting times, and so on, so that this 
tragic feedback can eventually trigger so high pres- 
sures that people are crushed or falling and trampled. 
It is, therefore, imperative, to have sufficiently wide 
exits and to prevent counterflows, when big crowds 
want to leave (Helbing et al., 2000b). 



4. Herding behavior 

Let us now discuss a situation in which pedestri- 
ans are trying to leave a room with heavy smoke, but 
first have to find one of the hidden exits (see Fig. 51). 
Each pedestrian a may either select an individual di- 
rection e* or follow the average direction {ep(t)) a of 
his neighbours (3 in a certain radius R a (Vicsek et al. , 
1995), or try a mixture of both. We assume that both 
options are again weighted with the parameter n' a of 
nervousness: 

e° a (t) = Norm [(1 - <) e* + n' a (e (t)) a ] , (128) 

where Norm(z) = z/||z|| denotes normalization of a 
vector z to unit length. As a consequence, we have 
individualistic behavior if n' a is low, but herding be- 
havior if n' a is high. Therefore, n' a reflects the degree 
of panics of individual a. 

FIG. 51. Herding behavior of panicking pedestrians in a 
smoky room (black), leading to an inefficient use of avail- 
able escape routes (from Helbing et al, 2000b). 

Our model suggests that neither individualistic nor 
herding behavior performs well. Pure individualistic 
behavior means that a few pedestrian will find an exit 
accidentally, while the others will not leave the room 
in time, before being poisoned by smoke. Herding may, 
in some cases, guide the mass to do the right things, 
but this requires that there is someone who knows a 
suitable solution to the problem, e.g. trained person- 
nel. However, in new and complex crises situations, 
herding behavior normally implies that the complete 
crowd is eventually moving into the same and proba- 
bly wrong or jammed direction, so that available ex- 
its are not efficiently used (see Fig. 51), in agreement 
with observations. According to simulation results by 
Helbing et al. (2000b) , optimal chances of survival are 
expected for a certain mixture of individualistic and 
herding behavior, where individualism allows some 
people to detect the exits and herding guarantees that 
successful solutions are imitated by small groups of 
pedestrians. 

VI. FLOCKING AND SPIN SYSTEMS 

Models for herding behavior have been developed 
already before in order to understand the flocking of 
birds (Reynolds, 1987), the formation of fish schools, 
or the herding of sheep based on a non-equilibrium 
analogue of ferromagnetic models (Vicsek et al., 
1995). More or less related models of swarm forma- 
tion are by Stevens (1992), Miramontes et al. (1993), 
Hcmmingsson (1995), Rauch et al. (1995), Toner and 
Tu (1995), Albano (1996), Bussemaker et al. (1997), 
Schweitzer (1997b), Stevens and Schweitzer (1997), 
and Mikhailov and Zanette (1999) (see also Sec. I.F). 

The possibly simplest model for the behavior of 
birds looks as follows. Let us assume all birds a 



are moving with a certain, finite speed ||u Q (£)|| = 
v a = v o- To make swarm formation less trivial, we 
set fa/3 = and do without the assumption of long- 
range attractive and short-range repulsive interactions 
among birds. Equation (4) can, then, be simplified to 
e a (t) = [e°(i) +Xa(t)], which is solved together with 



dx a [x a (t + At) — x a (t)] 



dt 



At 



= v e a (t). (129) 



We are now assuming that a bird a orients itself at 
the mean direction of motion Norm((e ( g(£ — At)) a ) = 
Norm((vp(t — At)) a ) of those birds (3 which are in 
a neighborhood of radius R a at time (t — At). As a 
consequence, the direction of motion e° is adapted to 
the average one in the immediate environment with a 
reaction time of At: e°(i) = Norm((e / g(£ — At)) a ) . In 
order not to have velocity changes, we assume that the 
fluctuations Xa (t) influence only the directions e a (t) , 
but stronger, the greater the fluctuation level r\ is. The 
vector [e°(i) + Xa(t)} is, therefore, also normalized. 
Consequently, the above self-propelled particle (SPP) 
model can be summarized by Eq. (129) and 

e a (t + At) = Norm(Norm( e/3 (i)) Q + X ' a (tj) , (130) 

where the fluctuations x' a (*) are uniformly distributed 
in a sphere of radius 77. 

FIG. 52. Typical simulation results for the two- 
dimensional, self-propelled particle model given by equa- 
tions (129) and (130). (a) Ferromagnetic ordering for 
vo = 0. (b) For vo > one observes a flocking of parti- 
cles moving into the same direction. (After Vicsek et al, 
1995.) 

The computer simulations start with a uniform dis- 
tribution of particles and mostly use periodic bound- 
ary conditions. Figure 52 shows two representative 
snapshots of the two-dimensional model variant for 
v a = and v 7^ 0. An exploration of the model behav- 
ior in three dimensions yields a phase diagram with an 
ordered phase and a disordered one (see Fig. 53) . For 
uo 7^ 0, no matter how small the average particle den- 
sity g is, there is always a critical fluctuation strength 
r)c(g), below which the directions e a (t) of motion are 
aligned and particle clusters are formed, even without 
attractive interactions (Czirok et al., 1999b). In other 
words, birds are flocking, if the fluctuation strength is 
only small enough. Otherwise, they are flying around 
in a disordered way. At a given fluctuation strength 
77, flocking is more likely, if their density g is high, 
which is in agreement with observations. Moreover, 
the model reproduces that birds are keeping a certain 
(i? Q -dependent) distance to each other, although no 
repulsive interactions have been assumed. 



FIG. 53. Phase diagram of the three-dimensional self- 
propelled particle (SPP) model and the corresponding fer- 
romagnetic system. The nonequilibrium SPP system be- 
comes ordered in the whole region below the curve w c (q) 
connecting the diamonds. In the static equilibrium case 
with vo = 0, the ordered region extends only up to a finite 
"percolation" density, see the beginning of the area given 
by the curve connecting the plus signs. (After Czirok et al, 
1999b.) 

While the critical fluctuation strength r) c (g) sepa- 
rating ordered and disordered behavior seems to be 
independent of the concrete velocity vo > 0, the be- 
havior of the related equilibrium system with vq = is 
considerably different. It corresponds to the behavior 
of a Heisenberg spin system with homogeneous particle 
density g, in which the spins e Q align parallel as in a 
fcrromagnet, if the temperature (which is proportional 
to the fluctuation strength 77) is below a critical tem- 
perature. However, this ferromagnetic ordering only 
takes place if the particle density g exceeds a certain, 
so-called percolation value g c w 0.75 (see Fig. 53). Sur- 
prisingly enough, the behavior for vq = and vq 7^ 
is, therefore, fundamentally different. 

In conclusion, non-equilibrium systems can consid- 
erably differ from equilibrium ones. The motion of 
the particles in a self-driven non-equlibrium system 
seems to have a similar effect as long-range interac- 
tions, since it is mainly a matter of time, until the 
particles meet and interact with each other. As a con- 
sequence, one can find a transition to an ordered phase 
even in an analogous one-dimensional non-equilibrium 
system (Czirok et al, 1999a), while the related equi- 
librium system with i> = is always disordered in the 
presence of fluctuations (Mermin and Wagner, 1966). 



A. Stock markets: Herding and oscillations 

As is reflected in bubbles and crashes of stock mar- 
kets, the reinforcement of buying and selling decisions 
has sometimes features of herding behavior (see, e.g., 
Yousscfmir et al, 1998; Farmer, 1998; Lux and March- 
esi, 1999). Remember that herding behavior is fre- 
quently found in complex situations, where individ- 
uals do not know what is the right thing to do. Ev- 
eryone counts on collective intelligence then, believing 
that a crowd is following someone who has identified 
the right action. However, as has been shown for pan- 
icking pedestrians seeking for an exit, such mass be- 
havior can have undesireable results, since alternatives 
are not efficiently exploited. 

The corresponding herding model in Sec. V.D.4 
could be viewed as paradigm for problem-solving be- 
havior in science, economics, and politics, where new 
solutions to complex problems have to be found in a 
way comparable to finding the exit of a smoky room. 
From the simulation results, we may conclude that 
people will not manage to cope with a sequence of 
challenging situations, if everyone sticks to his own 



idea egocentrically, but they will also fail, if everyone 
follows the same idea, so that the possible spectrum 
of solutions is not adequately explored. Therefore, the 
best strategy appears to be pluralism with a reason- 
able degree of readiness to follow good ideas of others, 
while a totalitarian regime would probably not sur- 
vive a series of crises. This fits well into experimental 
data on the efficiency of group problem solving (Ander- 
son, 1961; Kelley and Thibaut, 1969; Laughlin et at, 
1975), according to which groups normally perform 
better than individuals, but masses are inefficient in 
finding new solutions to complex problems. 

Considering the phase diagram in Fig. 53, we may 
also conjecture that, in the presence of a given level rj 
of fluctuations, there is a certain critical density above 
which people tend to show mass behavior, and below 
which they behave individualistically. In fact, the ten- 
dency of mass behavior is higher in dense populations 
than in dilute ones. In summary, conclusions from 
findings for self-driven many-particle systems reach 
far into the realm of the social, economic, and psy- 
chological sciences. 

Another example is the analogy between irregular 
oscillations of prices at the stock market and of pedes- 
trian flows at bottlenecks (see Sec. V.C.2). The idea 
is as follows: Traffic is a prime example of individu- 
als competing for limited resources (namely, space). 
At the stock exchange market, we have a competition 
of two different groups: optimistic traders ("bulls") 
and pessimistic ones ("bears"). The optimists count 
on growing stock prices and increase the price by their 
orders. In contrast, pessimists speculate on a decrease 
in the price and reduce it by their selling of stocks. 
Hence, traders belonging to the same group enforce 
each others' (buying or selling) actions, while opti- 
mists and pessimists push (the price) into opposite di- 
rections. Consequently, there is an analogy with oppo- 
site pedestrian streams pushing at a bottleneck, which 
is reflected in a roughly similar dynamics (Hclbing, 
2001). 

VII. SUMMARY AND OUTLOOK 

In this review, I have shown that many aspects 
of traffic flow and other living systems can be re- 
flected by self-driven many-particle systems, which 
implies that the respective interactions dominate 
socio-psychological effects. Although there are still 
many interesting open questions and controversial de- 
bates, calling for intensified research activities, one 
can already state that traffic theory is a prime exam- 
ple of a mathematically advanced, semi-quantitative 
description of human behavior. Meanwhile, it appears 
feasible to reproduce the observed complex transitions 
between different traffic states (see Sec. II. E) by simu- 
lation models which are adapted to the measurement 
site by variation of one parameter value only and work 
with the empirically measured boundary conditions 



(see Sec. IV.B.l). One could even say there are kind of 
natural laws behind the behavior of drivers and pedes- 
trians. This includes the self-organized, characteristic 
constants of traffic flow (see Sees. II.E.l and IV. A. 5), 
the collective patterns of motion observed in pedes- 
trian crowds (see Sees. V.C and V.D), and the trail 
formation behavior of humans (see Sec. I.F.5). 

Self-driven many-particle systems display a surpris- 
ingly rich spectrum of spatio-temporal pattern forma- 
tion phenomena. They were described by Newton's 
equation of motion (where Newton's third law ac- 
tio = reactio was sometimes relaxed), complemented 
by driving forces and frictional dissipation. Based on 
such small modifications, we often found consider- 
ably different behavior compared to analogous sys- 
tems from classical or equilibrium statistical mechan- 
ics, but there were also various analogies with driven 
fluids or granular media (see Sees. III.E.3, IV. A. 4 and 
V.A). Important factors for the finally resulting pat- 
terns were (i) the specification of the desired veloci- 
ties and directions of motion in the driving term, (ii) 
inhomogeneities and other boundary effects, (iii) the 
heterogeneity among the particles, and (iv) the degree 
of fluctuations. Many of the systems were character- 
ized by non- linear feedback mechanisms, so that small 
perturbations could have strong effects like "phantom 
traffic jams" (see Sees. II.E.l and IV. A. 1) or "phan- 
tom panics" (see Sec. V.D. 3). These effects were of- 
ten counterintuitive, such as "freezing by heating" 
(see Sec. V.D.I) or the "faster-is-slower effect" (see 
Sec. V.D. 2). Summarizing the discoveries, in the con- 
sidered self-driven dissipative many-particle systems 
with repulsive interactions, one can find all sorts of 
counter-intuitive breakdown and structure formation 
phenomena, when the "temperature" (i.e. the fluctu- 
ation strength) is increased. In contrast, from classi- 
cal many-particle systems with attractive interactions, 
we are used to the idea that increasing temperature 
breaks up structures and destroys patterns, i.e. fluid 
structures are replaced by gaseous ones, not by solid 
ones, etc. Nevertheless, the combination of repulsive 
interactions with a driving term and dissipation gives 
also rise to phenomena typical for attractive interac- 
tions, e.g. the segregation effects observed in the for- 
mation of traffic jams (see Sec. IV. A. 3) and pedestrian 
lanes (see Sec. V.C.I), or the emergence of the coher- 
ent moving state in a system of cars and trucks (see 
Sec. IV.D.2). 

The description of the phenomena in self-driven 
many-particle systems uses and generalizes almost 
the complete spectrum of methods developed in non- 
equilibrium statistical physics and non-linear dynam- 
ics. It has close relations with kinetic gas theory 
and fluid dynamics, soft matter, solid state physics, 
and even quantum mechanics (see Sec. III.C.l), while 
the corresponding non-equilibrium thermodynamics is 
still to be developed. However, there are aspects of 
universality found in traffic, as is exemplified by the 
power-law behavior of traffic flows (see Sees. II.E.l 



and IV. A. 3) and the phase diagrams of traffic dynam- 
ics and model characteristics (see Figs. 36, 37, and 40). 
In the spirit of general systems theory (Buckley, 1967; 
von Bcrtalanffy, 1968; Rapoport, 1986), one could 
even go one step further: Many phenomena found in 
traffic systems, where drivers or pedestrians compete 
for limited space, may have implications for other bio- 
logical and socio-economic systems, in which individ- 
uals compete for limited resources. In fact, jamming 
and herding phenomena are found in various living 
systems, from animal colonies, markets, administra- 
tions, or societies, up to science, economics, and pol- 
itics. Therefore, traffic is a good and, in particular, 
concrete example to study certain aspects of socio- 
psychological phenomena in an experimental setup. 
Possibly, this can help to gain a better understanding 
of more complex human behavior in the future. 

Self-driven many-particle systems are also in- 
teresting from a practical point of view. They 
are related to aspects of collective intelligence 
(see, e.g., the work carried out at the Santa Fe 
Institute, http : //www . santaf e . edu/sf i/research/ 
indexResearchAreas.html, by scientists like Botcc 
and Bonabeau, 1998; Kube and Bonabeau, 1998; 
Wolpert and Turner, 1999). Based on suitable interac- 
tions, many-particle and many-agent systems can or- 
ganize themselves and reach an optimal state (Hclbing 
and Vicsek, 1999). Such a distributed strategy based 
on local optimization is more efficient and robust than 
classical, centralized control strategies. In particular, 
the outcome on large scales is relatively insensitive to 
the local failure of control. 
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